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Losses and energy efficiency 
of drive motors and systems 


Replacement of the Sankey diagram of power decrease 
in the direction of power flow by a diagram of power 
increase opposite to the direction of power flow opens a new 
perspective of research of drive motors and systems 


Zygmunt Paszota, Prof. 
Gdansk University of Technology, Poland 


ABSTRACT 


Losses and energy efficiency of every drive motor must be presented as functions of 
physical quantities independent of losses in the motor. Such quantities are speed and load 
required by the machine or device driven by the motor, changing in the drive operating 
field. Speed and load of the motor decide of the instantaneous useful power of the motor 
and also in a differentiated way of kinds and values of losses occurring in the motor. 
However, losses and energy efficiency of the hydrostatic drive motors and systems are 
evaluated by researchers and manufacturers as functions of parameters depending on 
the losses. The basic cause of such situation is the traditional, commonly accepted but 


erroneous, view of the power flow in the drive motors and systems represented by the Sankey diagram 
of power decrease in the direction of power flow. It is necessary to replace the Sankey diagram by the 
proposed diagram of increase of power in the motor and in the drive system in the direction opposite to 
the direction of power flow. The proposed view of losses and energy efficiency should be applied to all 
types of motor and drive systems. The aim of this paper is showing the resulting problems of the above 
postulates, exemplified by operation of a rotational displacement motor in a hydrostatic drive system. 
In order to make possible objective evaluation of the energy behaviour of different motor and system 
solutions and sizes, the losses and energy efficiency should be described and compared as dependent 
on the motor speed coefficient @,, and load coefficient M changing in the drive system operating field 
(0 <@y< Dumao OM u< M ymad. The presented proposals open a new perspective of unavoidable research 
of drive motors and systems, making it possible to compare objectively the energy efficiency of different 
types of motors and drive systems. 


Key words: drive motors; drive systems; energy losses; energy efficiency; 
Sankey diagram; new research field 


INTRODUCTION 


Losses and energy efficiency of every drive motor must 
be presented as functions of physical quantities independent 
of losses in the motor. Such quantities are speed and load 
required by the machine or device driven by the motor, 
changing in the drive operating field. Speed and load of 
the motor decide of the instantaneous useful power of the 
motor and also in a differentiated way of kinds and values 
of losses occurring in the motor. 

However, losses and energy efficiency of the drive motors 
and systems are evaluated by researchers and manufacturers as 
functions of parameters depending on the losses. An example 


of the wrong interpretation may be presentation of energy 
efficiencies of the hydrostatic rotational displacement motors 
as dependent on the flow intensity of the motor feeding liquid 
and on the pressure decrease of the motor. 

Energy efficiency of turbines is presented in a similar way. 

A cause of such situation is the traditional, commonly 
accepted but erroneous, view of the power flow in the drive 
motors and systems represented by the Sankey diagram of 
power decrease in the direction of power flow. 

It is necessary to replace the Sankey diagram by the 
proposed diagram of increase of power in the motor and in 
the drive system in the direction opposite to the direction of 
power flow [1 + 23]. 
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The proposed view of losses and energy efficiency 
should be applied to all types of motor and drive systems, 
for instance: 

— piston internal combustion motors, 
— turbines, 

— electric motors, 

— hydrostatic displacement motors, 
— ship screw propellers. 


The aim of this paper is showing the resulting problems of 
the above postulates, exemplified by operation of a rotational 
displacement motor in a hydrostatic drive system. 

In order to make possible objective evaluation of the 
energy behaviour of different motor and system solutions and 
sizes, the losses and energy efficiency should be described and 
compared as dependent on the motor speed coefficient © q and 
load coefficient My changing in the drive system operating field 
(0 SOMS Ommay 0< Mm z Muimax)- 

The presented proposals open a new perspective of 
unavoidable research of drive motors and systems, making 
it possible to compare objectively the energy efficiency of 
different types of motors and drive systems. 


INDEPENDENT AND DEPENDENT 
PARAMETERS OF A DRIVE MOTOR AND 
SYSTEM OPERATION, EXAMPLIFIED 
BY OPERATION OF A DISPLACEMENT 
MOTOR IN A HYDROSTATIC DRIVE 
SYSTEM 


The work of a rotational or linear hydraulic motor as an 
element of hydrostatic drive and control system, directly 
connected with the driven machine (device) must provide 
parameters required by the driven machine (angular speed @,, 
(rotational n,,) of the shaft or linear speed v,, of the piston rod 
and shaft load (torque) M,, or piston rod load (force) Fy) and 
also ensure the required machine movement direction. 

The mechanical parameters of a motor (speed ,, (Ny) or 
Vy and also load M,, or Fy) change in the range from zero to 
maximum Values Oumas (Ayimax) OF Vaimax Nd Mymax OF Fuma 

The required current speed Oy (Ny) or Vy and required 
current load M,, or Fy of the driven machine are an effect 
of its work cycle and the work task. The current driven 
machine speed and load values are independent of the type 
and structure of that machine driving system control (e.g. 
an electrical or hydrostatic system). 

The current speed and current load of a hydrostatic 
system driven machine have a direct or indirect impact 
on the mechanical, volumetric and pressure losses in the 
hydraulic motor, pump and other system elements, a system 
with determined motor speed control structure. The losses 
are also an effect of the viscosity of the used working liquid 
(hydraulic oil). 

The current speed @,, (ny) or Vy and current load My 
or F,, of the driven machine influence, in consequence, 
the current hydraulic motor absorption capacity Qu and 
pressure decrease Ap,, and also (depending on the used 
motor speed control structure) the current pump capacity 
Q, and discharge pressure pp»). 

Ifin effect of the increasing, required by the driven machine 
(device) hydraulic motor speed @\y (Ny) Or Vy, as Well as in 
effect of the increasing, required by the driven machine motor 
load M,, or Fy, and also in effect of the mechanical, volumetric 
and pressure losses of the hydrostatic drive system elements, the 
maximum drive system capacity (determined by the maximum 
pump capacity Q,,,,, Or Maximum pump discharge conduit 
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pressure Pp» nx Limited to the system nominal pressure p,) is 
fully used, then further increase of @,, (Ny) or Vy as well as 
My, or Fy, will not be possible. 

Maximum pump capacity Qpmax 1s smaller than its theoretical 
capacity Q,,. The pump theoretical capacity Q,, is a product of 
the theoretical capacity qp, per one pump shaft revolution and 
the no-load pump shaft speed npọ. The pump Qpmax Capacity, 
however, results from the loaded pump speed n, , lower than 
the speed np, and from volumetric losses in the pump. 

The system nominal pressure p, is a maximum permissible 
continuous operation pressure Pp,,,, determined in the pump 
discharge. 

The maximum speed values Oyma; Mymax) OF Vyuma as Well 
as the maximum load values Mymax OF F ymas Of the hydraulic 
motor used in a hydrostatic drive system are limited 
by the maximum pump capacity Q,,,,, or by the system 
(pump) nominal pressure p, and also by the corresponding 
mechanical, volumetric and pressure losses in the remaining 
system elements, the losses being also an effect of viscosity of 
the working liquid used. Therefore, the © ymax (ymax) OF Vmax? 
Mymax OF F ymax Values are dependent variables. 

The current mechanical operating parameters of the 
hydraulic motor used in a hydrostatic drive system (current 
motor speed @,, (Ny) OF Vy, and current motor load My or Fy) are 
independent values in the motor, deciding of the motor losses 
and of the hydraulic parameters (the current motor absorbing 
capacity Q,, and current pressure decrease Ap,, also depending 
on the motor mechanical, volumetric and pressure losses). The 
current motor absorbing capacity Q,, and current pressure 
decrease Ap,, are dependent variables in the motor. 

In the hydraulic motor (hydrostatic drive system) operating 
field (0 < Oq (My) < Ommas ymax)» 9S Mu <Moimax) OF (0 < Vy< 
< Vumaw Ô < Fu < Fuma), the pressure and flow intensities in 
the system and also the energy losses in the motor, in the pump 
and in the whole system, power of energy losses and energy 
efficiencies of the system elements should be considered the 
functions of the current speed @,, (Ny) or Vy and the current load 
My, or Fy, required by the system driven machine (device). Also 
the torque M, that the pump loads the driving (electric or 
internal combustion) motor and the speed n, that the motor 
drives the pump with should be considered the functions 
of the current speed and the current load required by the 
system driven machine. 

The decrease of speed n, that the electric or internal 
combustion motor drives the pump with is connected with the 
increase of torque M, that the pump loads the motor with. The 
decrease of speed depends on the operating characteristics of 
the motor, which is not a component of the hydrostatic drive 
system. Therefore, the pump driving speed n, should be 
treated as a parameter independent of the hydrostatic 
system (of the pump). 


NON-DIMENTIONAL COEFFICIENTS 
OF THE DISPLACEMENT MOTOR 
PARAMETERS, COEFFICIENTS OF 

ENERGY LOSSES IN THE HYDROSTATIC 
DRIVE ELEMENTS 


The energy efficiency of the hydrostatic drive system and its 
elements is described by mathematical models as functions of 
the hydraulic motor (hydrostatic drive system) speed coefficient 
®w and load coefficient My. 

The current angular speed @,, (rotational speed ny) 
required of a rotational motor or the linear speed v,, 
required of a linear motor, operating in a hydrostatic drive 
system, are replaced in the energy efficiency mathematical 


models by the motor speed non-dimensional coefficient 
On: 


Oy = Ou _ Om _©m Ime _ MM Mm 
Om Ny MQ Qs, 
or 
— Vm _ Vm Sui 
v= 
Vt Qrt 


The rotational hydraulic motor speed coefficient ©y is 
a ratio of the current angular speed @,, (rotational speed ny), 
required of the motor by driven machine, to 
2T1Q), 


mt 
Qr 


Mt 

which would correspond with the theoretical capacity Qp, of the 
motor driving pump and with the theoretical motor absorbing 
capacity qm per one shaft revolution. The speed Om (Nm) would 
be achievable on the condition, that there are no volumetric 
losses in the hydrostatic drive system (including the pump and 
the hydraulic motor) and the pump is driven by an (electric or 
internal combustion) motor operating with constant rotational 
speed n, = np independent of its load. 

The theoretical angular speed @,,, (rotational speed ny) 
of a rotational motor is treated as a constant reference value 
for the motor current angular speed @,, (rotational speed ny). 

The linear hydraulic motor speed coefficient @y, is a ratio 
of the current linear speed v,,, required of the motor by driven 
machine, to 

Qer 


MI 

which would correspond with the theoretical capacity Qp, of 
the motor driving pump and with effective area S,,, of the 
motor piston in the inlet chamber. The speed v,,, would be 
achievable on the condition, that there are no volumetric losses 
in the hydrostatic drive system (including the pump and the 
hydraulic motor) and the pump is driven by an (electric or 
internal combustion) motor operating with constant rotational 
speed np = Np independent of its load. 

The theoretical linear speed v,,, of a linear motor is 
treated as a constant reference value for the current motor 
linear speed vy. 

The current torque M,, required of a rotational motor 
or current force F,, required of a linear motor, operating in 
a hydrostatic drive system, are replaced by the motor load 
non-dimensional coefficient My: 


theoretical angular speed Oy, = 


(theoretical rotational speed Ny, = 


theoretical linear speed Vy, = 


M, = Mu B 2I Myu 
Mm QmtPn 
or 
— F F 
Mu = —— = M 
Fue SwiPn 


The rotational hydraulic motor load coefficient My, is 
a ratio of the current torque My, required of the motor by 
driven machine, to 


theoretical torque M,,, = Amt Pn n 


2I 
which would correspond with the theoretical absorbing 
capacity qm per one motor shaft revolution and with the 
hydrostatic system nominal pressure p,. The torque Mw would 


be achievable on the condition that there are no mechanical or 
pressure losses in the hydraulic motor and in the remaining 
system elements (except the pump) and the pressure pp, 1n 
the pump discharge conduit is equal to the system nominal 
pressure p,- 

The rotational motor theoretical torque M wis treated as 
a constant reference value for the current motor torque My. 

The linear hydraulic motor load coefficient My is 
a ratio of the current force Fy, required of the motor by driven 
machine, to 


theoretical force Fy = Sy, Py» 


which would correspond with the effective area Sọ, of the 
motor piston in its inlet chamber and with the system nominal 
pressure p,. The force Fw, would be achievable on the condition 
that there are no mechanical or pressure losses in the hydraulic 
motor and in the remaining system elements (except the pump) 
and the pressure Pp» max 1N the pump discharge conduit is equal 

to the system nominal pressure p,. 

The linear motor theoretical force Fw is treated as 

a constant reference value for the current motor force Fy. 
The mechanical, volumetric and pressure losses in 

a hydraulic motor, pump and in the remaining hydrostatic 
drive system elements are described in the mathematical 
models of the losses, power of losses and energy efficiency 
by the coefficients k; — relations to the values connected 
with the values of the hydrostatic drive system characteristic 
parameters: 

— theoretical capacity qp, per one pump shaft revolution, 

— theoretical absorbing capacity qy, per one rotational 
hydraulic motor revolution or effective piston area Sy; in 
the linear motor inlet chamber, 

— theoretical pump capacity Q,,, 

— system nominal pressure p,. 


The basis of energy evaluation of the particular design 
solutions and size of the volumetric machines is a catalogue 
of the coefficients k, of energy losses in various types of 
pumps and hydraulic motors used in the hydrostatic drive 
systems, operating with different levels of pump theoretical 
capacity Q,, and system nominal pressure p,, with the 
working liquid reference viscosity v,. 


THE MOTOR OPERATING FIELD 
IN A DRIVE SYSTEM 


Figure | presents the operating field of a rotational or 
linear displacement motor in a hydrostatic drive system. The 
operating field is determined in the plane of motor mechanical 
parameters, i.e. speed coefficient Oy and load coefficient My, 
independent of the motor and of the system. 

The limit values Ummas = f(My) or Mumax = f©m) of the 
hydraulic motor operating field are determined by the maximum 
motor feed capability in the hydrostatic drive system. The 
values Oyma and Mygnax are dependent on the motor and 
on the system losses. 

The maximum motor absorbing capacity Q mma 
achieved in the system by the applied motor speed control 
structure, should be equal or close to the instantaneous 
maximum pump capacity Q,,,,, (resulting from the theoretical 
capacity Q,, decrease of the pump shaft rotational speed n, and 
the intensity of pump volumetric losses Q,,). 

The maximum possible motor pressure decrease Apyjmax 
should be equal or close to the system nominal pressure p, 
determined in the pump discharge conduit, reduced by the 
pressure losses Ap, in the system conduit. (In the motor series 
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Mumae =f( Qu, ki, V/Un ) 
A. Tuma depends on the current 3 value, 
S coefficients k, of the mechanical 


hydraulic motor (hydrostatic drive system) load coefficient My 


the 


Fo es a the machine working point (GysMy) 


7 


Sar, OL we el tye es er Ser see Mak ay ce pe ge e's 


` the machine working point with instantaneous values © >° 7 7 
* of the speed coefficient Wy and load coefficient My 
` required by the machine from the hydrostatic drive system; 
* fhe current values Gy,My are generally independent of each 
` other and independent of the hydrostatic drive system losses ` ` 


and pressure losses in the hydrostatic 
system elements (excep! the pump) 
and on the ratio v/v, of working liquid 
viscosity v to reference viscosity uy, 


@umax= f (Mus kis V/Un ) 

Gumo depends on the current 
My value, coefficients k; of 
the volumetric losses 

in the hydrostatic system 
elements and on the ratio 
v/v. of working liquid 
viscosity v to the reference 
viscosity vs 


0 the hydraulic motor (hydrostatic drive system) speed coefficient Oy 1 


Fig. 1. The range of motor speed coefficient ©, and load coefficient My, (0 < Dy < Duna OX 


throttling speed control structure, the maximum slot area of the 
throttling valve, proportional directional valve or servo-valve 
should allow to minimize the pressure decrease APproyma With 
the set Quimax ~ Qpmax): 

Therefore the limit values O\ymax Of the hydraulic motor 
speed coefficient are a function of the current motor load 
coefficient My, coefficients k; of the volumetric losses in the 
hydrostatic system elements (including coefficient k, of the 
pump shaft rotational speed decrease An,) and a function of 
the ratio v/v, of working liquid viscosity v to the reference 
viscosity V,,. 

On the other hand, the limit values Myjmax Of the hydraulic 
motor load coefficient are a function of the current motor 
speed coefficient Om, coefficients k; of the mechanical and 
pressure losses in hydrostatic system elements and a function 
of the ratio v/v, of working liquid viscosity v to the reference 
viscosity v,. 


DIRECTION OF POWER INCREASE 
IN A DRIVE SYSTEM 


The Sankey diagram is the fundamental cause of the 
incorrect evaluation of losses in a drive system. 

The Sankey diagram connected with a drive system 
suggests an evaluation defining the output power of the system 
as a difference between the system input power and sum of the 
powers of losses occurring in the system. Therefore the energy 
losses in the system are evaluated as a function of parameters 
describing the input power. 

In accordance with the proposed by the author diagram of 
power increase in a drive system opposite to the direction of 
power flow, the system input power is a sum of the system 
output power and powers of losses occurring in the system. 
The energy losses in the system are evaluated as a function 
of parameters describing the system output power. 

Full picture of the energy losses in a drive system is 
a picture of power of energy losses in the system elements. 
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fu < Mumay) in a hydrostatic drive system 


Shaft power of the hydrostatic system feeding pump is equal 
to the sum of hydraulic motor shaft (or piston rod) power 
and powers of individual losses in the power flowing from the 
pump shaft to the hydraulic motor shaft (or piston rod). 


The quantity of power increases, in order to overcome 
the power of energy losses, in the opposite direction to 
the direction of power flow. Therefore, the image of power 
of energy losses in the system should be constructed in the 
direction from the hydraulic rotational motor shaft or from 
linear motor piston rod towards the system feeding pump 
shaft. 

Power of the hydraulic motor, pump and of other system 
element losses (mechanical, volumetric, pressure losses) should 
be determined as functions of the parameters independent of 
those losses and influencing directly those losses. 

Powers of energy losses in the system elements and also 
power of those elements must be precisely defined. 


DIAGRAM OF POWER INCREASE 
IN A ROTATIONAL DISPLACEMENT 
HYDRAULIC MOTOR, REPLACING THE 
SANKEY DIAGRAM 


Figure 2 illustrates the diagram, proposed by the Author, 
of power increase in a rotational hydraulic motor opposite to 
the direction of power flow, replacing the Sankey diagram of 
power decrease in the direction of power flow. 

The power Pw. consumed by the hydraulic motor is a sum 
of motor shaft useful power Pyu, and powers of three different 
energy losses in the motor. The losses occur in series increasing 
power in the opposite direction to the direction of power flow. 
In effect, the power in the motor increases from the shaft 
useful power Py, to the working liquid power Pu. consumed 
by the motor: 


Sankey diagram 


direction 
of power flow 


inlet channel 


“shaft = working 
chambers assembly 


—) 
direction” 
of power flow 


outlet channel 


direction 
of power flow 


un Fuc a AFup -Afw -A 


diagram according to Z.Paszota 


RP 
irection 
í power flow 


viscosity » of the working fluid 


inlet channel (Ap,,,) Ru R Ra =(Pur Per) UFAN Raa + AR 
py P 


ARF APAS 
Apt Pupu "shaft - working 
chambers” assembly 


of power flow 


\working chambers” shoft speed ny (cy) 


= shaft torque My, 


~ 


outlet channel (Ap,,2) = 
rection 
f power flow 


R= Pu Au 


Fue” Fiu + ARim + Ay + Aip 


a0 


outlet pressure py *0 


Fig. 2. Diagram of power increase in a rotational displacement hydraulic motor, opposite to the direction of power flow, replacing the Sankey diagram of 
power decrease in the direction of power flow (example of motor with theoretical (constant) capacity qn Vu) per one shaft revolution) 


Power increases from the motor useful power P „, required on the motor shaft by the driven machine (device) 
to power P,,, consumed and required by the motor of the working liquid. 
The increase of power is an effect of the powers of losses in the motor: power AP „„ of mechanical losses in the „shaft - working chambers” 
assembly, power AP „ of volumetric losses in the working chambers and power AP „, of pressure losses in the channels. 
Powers AP y» AP x, and AP y of the losses are functions of the output parameters of the motor assembly where the losses occur and diversified functions 


of the working liquid viscosity v: power AP q, 


im 


of mechanical losses is a function of torque M,, and shaft speed n, (@,) required of the motor by the driven 


machine (device) and a function of the working liquid viscosity v, power AP „ of volumetric losses is a function of the pressure decrease Ap, indicated in 
working chambers (of torque M,, indicated in the chambers) and of the shaft rotational speed n,, as well as a function of the working liquid viscosity v, 
power AP,,, of pressure losses is a function of motor capacity Qy, and of the working liquid viscosity v. 
Power P; indicated in the working chambers: Py = Py, + AP im» power P ye; Of the working liquid consumed in the working chambers: 
Pii = Piu + AP ym + AP, power Py, of the working liquid consumed by the motor: 


P 


Me P iu t AP yim t AP t AP yp 


The proposed diagram replaces the Sankey diagram of distribution of power in hydraulic motor causing 
incorrect loss evaluation during the hydraulic motor energy investigations. 


Mechanical losses (and power AP wm of mechanical losses) 
occur in the ,,shaft - working chambers” assembly. 

Volumetric losses (and power AP,,, of volumetric losses) 
occur in the working chambers. 

Pressure losses (and power AP,,, of pressure losses) occur 
in the channels. 


COMPLEX DEPENDENCE OF 
ENERGY LOSSES IN A ROTATIONAL 
DIPLACEMEMENT HYDRAULIC MOTOR 
ON THE SHAFT SPEED, ON THE SHAFT 
TORQUE, ON THE CAPACITY PER ONE 
SHAFT REVOLUTION AND ON THE 
WORKING LIQUID VISCOSITY 


There exists a direct dependence of the torque Mum of 
mechanical losses in the ,,shaft - working chambers” assembly 
on the torque M,,, on the motor shaft rotational speed n,,, on 
the capacity qm (OF quey) per one shaft revolution as well as on 
the working liquid viscosity v. 

There is a complex dependence of the intensity Qu, of 
volumetric losses in the working chambers on the shaft loading 
torque M,, and on the torque Mum of mechanical losses in 
the „shaft - working chambers” assembly (decrease Apy; of 
pressure indicated in the working chambers depends on M,, and 
Mym ON qm: (OF qmeyv) and has direct impact on Quy) and also on 
the shaft speed n,, (influencing in a diversified way the torque 


Mym Of mechanical losses and intensity Qu, of volumetric 
losses). The intensity Qu, of volumetric losses depends on 
a diversified impact of the working liquid viscosity v: indirectly 
by impact of v on the torque Mun of mechanical losses in the 
„Shaft - working chambers” assembly and directly by impact of 
v on the intensity Q,,, of losses in the working chambers. 

It can be said, that it is a complex dependence of pressure 
losses Ap,,, in the channels on the shaft rotational speed ny, 
on the capacity qm (OF dmgv) per one shaft revolution and on 
intensity Q,,, of volumetric losses in the working chambers. The 
intensity Qu, of volumetric losses influences the motor capacity 
Q,, and at the same time Qu, depends in a complex way on the 
shaft loading torque M,, and on the torque Mym of mechanical 
losses in the ,,shaft - working chambers” assembly. Pressure 
losses Ap,,, in the motor channels are also dependent on the 
diversified impact of the working liquid viscosity v: indirectly 
by impact of v on the torque M,,,, of mechanical losses in the 
„Shaft - working chambers” assembly and by impact of v on 
the intensity Q,,, of volumetric losses in the working chambers 
and directly by impact of v on the losses Ap,,, of pressure in 
the channels. 

Contrary to the commonly used, both by manufacturers and 
researchers, methods of evaluation of the rotational hydraulic 
motor losses, it is unacceptable to create a ,,sum” of the torque 
Mym Of mechanical losses in the „shaft - working chambers” 
assembly and the ,,torque” of pressure losses Ap,,, in the motor 
channels, and also such a,,sum” must not be evaluated as 
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directly dependent on the same chosen parameters, because 
those losses are of different character and depend on different 
parameters: 


Mym = £ (My Dms qmi (dma). V) 
APup = f (Qus V) 


The impact of hydraulic oil viscosity v on the energy 

losses in a hydraulic motor, i.e. on: 

— torque Mym of mechanical losses in the „shaft - working 
chambers” assembly, 

— intensity Q,, of volumetric losses in the working 
chambers, 

— pressure losses Ap,,, in the channels, 

is diversified. 


Dependence of the motor losses on the hydraulic oil 
viscosity v should be presented in expressions describing also 


the dependence of those losses on other parameters which 
influence them directly: 


Mym = f (Mms Dms qmi (dmav)» V) 
Qw = f (APmis Nu» V) 
APwp = f (Qus v) 
The motor overall efficiency ny, as a function of 


Mw; Dy; qm (ama) and v, is a product of Nim Ny and Nimp 
efficiencies: 


Nm = f (Mys Dy» Amt (miqv)sV) = 


— Pu | Mum _ 2TIMyny 7 
Pue  APmQm ApmQu 
z Mmm Moy Mp 


where: 
P,,, is the motor useful power, 
Puc is the motor consumed power. 


Each of the three efficiencies, as a factor in the product 
describing the overall efficiency, is evaluated as a function 
of parameters directly influencing the respective losses and 
a function of parameter to which the losses are ,,added”: 

— motor mechanical efficiency Nym: 


(Mu + Mum) Om 

_ STM Bing 
~ QTM yp +My) Day 

— Mu 
Mu +Mun 


= f (Mumm (amv )- Y) 


where: 
P,, is the power indicated in the motor working chambers: 


Pui = AP Mid (Amg) Om = 
= (Mu + Mum) Om = 2M u +M um) Ds 
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— motor volumetric efficiency Ny: 


Pui AP Miu (Ame) Am 


Nav = SS OO 
Pui AP: (Gan (dug ) Nm + Quv) 


_ AP Mime (Amg) Am — m (dmg) Am 


AP 4iQua AmlaAmg) Oa Qw 


2 dm (dmev ) Ny 
Qu 


where Pua is the power consumed in the motor working 
chambers: 


= f (AP i> doe (Guiev ).Ty Vv) 


Puci = APumiCqmlAmg Nm + Quv) = APMiQu 
— motor pressure efficiency nyp: 
Pees AP ui (Ami (Amey) Am + Quv) 
Pue (Api + APup) (Am lAmg) Am + Quv) 


AP mi 


— ÂPmiQm _ 
APmi + AP up 


APy Qu 


AP ui 
=—“‘=f(A a : 
Apu (APmi» Qu» V) 


In order to present the motor volumetric efficiency ny, as 
a factor in the Numm Nmv Nmp Product describing ny, i.e. to present 
Ny as a complex dependence on the (My, nw qm:lAmev) V) 
parameters describing n,, and dependent on the mechanical 
losses, the intensity Quy = f(APis Is Ime(Amev),V) Of volumetric 
losses in the working chambers should be determined with 


211 (Ma, tM um) 
M7. nA 
dt Mev) 


and with torque Mum of mechanical losses in the „shaft 
- working chambers” assembly as an Mym = f (My, 1y, 
dm(Amev)-V) function. 

In order to present the motor pressure efficiency ny, as 
a factor in the Num Ny Mp Product describing nm, i.e. to present 
mp aS a complex dependence on the (My, Ams qmlAmev)Y) 
parameters describing n,, and dependent on the mechanical 
and volumetric losses in the motor, the pressure losses Apy,, = 
= f (Qus V) in the channels must be determined with 


Ap 


Ou= Amlqmev)Am + Quv 


then intensity Qu, = f (APmi Dy, V) of volumetric losses in the 
working chambers must be determined with 


= 211 (My +M um) 


AP ui = 
qdm(q Mev) 
and the torque Mum of mechanical losses in the „shaft - working 
chambers” assembly must be determined as an Mym = f (My, 
ny» AmlAmev)» V) function. 
The characteristic of the hydraulic motor overall 


efficiency Ny = f (Mys Nys AmlAmgv)» V) presents a complex 


picture as a product nmm Ny Nmp Of three efficiencies 


correctly described by: 
- mechanical efficiency Nym = f (Mys Dys Gni(Qmev)> Y), 


- volumetric efficiency Ny, = f (APyis AmlAmev)» nvs Y) 
- and pressure efficiency nyp = f (APyis Qu V)- 


The picture of the hydraulic motor overall efficiency ny = 
= f (Mms Dv: qmi (Amev) V) must be supplemented by assessment of 
the hydraulic motor operating field (0 < y | Omma > 0 < My < 
( Myimax) in the hydrostatic drive system, i.e. assessment of 
the range of ny and My (Op and My) parameters. 

In the motor (and the hydrostatic drive system) (0< © M| OMmax 
0 < Myu | Mumma Operating field, its current speed ny, (Oy) 
and load My (My) are an effect of the demand of the motor 
(system) driven machine (device) and are independent of 
the losses in the hydraulic motor and in the motor driving 
hydrostatic system. 

However, the hydraulic motor (system) operating field 
limit values nNymas (Mma) and Mymas (Mima) depend on 
the maximum capacity of the hydraulic motor driving system. 
The values NMmax (OMmax) and Mmmax (Myimax) determine 
simultaneously the corresponding motor energy efficiency ny 
and the overall system efficiency n. 

The limit parameters of the hydrostatic drive system 
operation result from the pump theoretical capacity Qp 
and the system nominal working pressure p, as well as 
from the actual energy losses in the hydraulic motor, 
conduits and pump and also losses in the motor speed 
throttling control assembly (if it is installed). Therefore, 
the limit values nymas (@ygmax) and Mymax (Mima) are also 
dependent on the working liquid viscosity v changing in the 
Vmin = V < Vmax range. The hydrostatically driven hydraulic 
motor operating field is also influenced by the n,= f(M,) 
characteristic of the (electric or internal combustion) motor 
in the pump driving system. 


NECESSITY OF USE OF THE ENERGY 
LOSS MATHEMATICAL MODELS 
WITH THE LOSS COEFFICIENTS IN 
DISPLACEMENT MOTOR AND INA 
HYDROSTATIC DRIVE SYSTEM 


Evaluation of the hydraulic motor overall efficiency ny = 
= f (Myp Dys qme (Amey), V) as a product Nym Novy Mvp Of three 
motor efficiencies can be performed only by means of the 
mathematical models of losses and efficiencies, where the 
defined coefficients k, of energy losses in the motor and in 


the motor drive system are used. 
Evaluation of the hydraulic motor energy efficiency 


is performed together with evaluation of the efficiency of 
a hydrostatic drive system where the hydraulic motor is 
used (including also the energy efficiency of pump, conduits 
and the hydraulic motor speed throttling control assembly 
(if it is used)). 

In the proposed method, based on the mathematical 
models of losses, each kind of energy losses is a function of 
parameters directly influencing the losses and independent 
of those losses. 

Evaluated are the values of the k; coefficients of energy 
losses, relating the mechanical, volumetric and pressure losses 
in the hydraulic motor, pump and other system elements to the 
reference values of driving system: nominal pressure p, of the 
system, theoretical capacity Qp, of the system driving pump, 
theoretical torque M,, of the pump shaft as well as theoretical 
torque M,,, of the hydraulic motor shaft. The k, coefficients are 
determined at the hydraulic oil reference viscosity v,. At the 


same time the impact is determined of the viscosity ratio v/v, in 
the Vmin < V < Vmax range on each kind of energy losses. 

The method allows to evaluate the values and proportions 
of mechanical, volumetric and pressure losses in the hydraulic 
motor, pump, conduits and in the throttling assembly (if 
installed), as well as the dependence on the hydraulic oil 
viscosity v. 

The energy investigations of a pump and hydraulic 
motor as independent displacement machines are limited 
to determination of the k, coefficients of losses in them. 
The energy efficiency characteristics of those machines 
are determined in parallel with efficiency evaluation of the 
hydrostatic drive system where they are used. 

The knowledge of k, coefficients of the mechanical, 
volumetric and pressure losses in the drive system elements 
allows to obtain, with the numerical method, the 
characteristics of the hydrostatic drive system overall 
efficiency, pump efficiency, hydraulic motor efficiency, 
conduit efficiency and the motor speed throttling control 
assembly (if it is used) structural efficiency in the motor 
(system) (0 < ðm l ymax > 0 < My ( Mymax) Operating field 
at a selected ratio v/v, of the hydraulic oil viscosity v to the 
reference viscosity v,. 

Characteristics of the overall efficiency of elements used 
in a hydrostatic drive system: pump np, hydraulic motor 
Ny conduits yn, and throttling control assembly (if it is 
used) structural efficiency n,, are defined as functions of 
the hydraulic motor (system) speed coefficient ®, and load 
coefficient My and the hydraulic oil viscosity ratio v/y,. 

At the same time the hydraulic motor (0 < y ( Oymax> 
0 < My | Myma,) Operating field in the hydrostatic drive 
system is determined at the selected ratio v/v, of the 
hydraulic oil viscosity to the reference viscosity. 

Characteristics of energy efficiency of the pump and 
hydraulic motor with the determined constant coefficients k, 
of losses and also of efficiency of the conduits are different 
in hydrostatic systems with different motor speed control 
structures. 

The method is precise by definition and simple in use. 
It simplifies the laboratory investigation of pumps, hydraulic 
motors and hydrostatic drive systems. It allows to seek for 
energy saving solutions of pumps and hydraulic motors. It 
allows also to evaluate the overall energy efficiency of the drive 
and to find energy saving hydrostatic drive system structures. 


CONCLUSIONS 


1. Losses and energy efficiency of every drive motor 
and system must be presented as functions of physical 
quantities independent of losses in the motor and 
system. Such quantities are speed and load required 
by the machine or device driven by the motor, 
changing in the (0 < Oy < Ommas > 0 < Mm < Mmmax) 
drive operating field. Speed and load of the motor decide 
of the instantaneous useful power of the motor and also in 
a differentiated way of kinds and values of losses occurring 
in the motor. 

2. In energy considerations of a motor and a drive system, it 
is necessary to replace the Sankey diagram by the proposed 
diagram of power increase in the direction opposite to power 
flow. 

3. Itis necessary to determine the motor operating field in its 
drive system, i.e. to determine the motor speed coefficient 
® wand load coefficient My dependent on and independent 
of losses in the motor and in the system. 

4. The example of operation ofa rotational displacement motor 
in a hydrostatic drive system shows a complex relation of 
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energy losses in the motor and in the system to the motor 
shaft speed and load, to the capacity per one shaft revolution 
and to the working liquid viscosity. Evaluation of the motor 
overall efficiency ny as a product of mechanical efficiency 
Nm > Volumetric efficiency Nym, and pressure efficiency ny, 
can be performed only by means of mathematical models of 
losses and efficiencies with the use of defined coefficients 
of losses in the motor and in the motor driving system. 
The presented proposals open a new perspective of 
unavoidable research of drive motors and systems, 
making it possible to compare objectively the energy 
efficiency of different types of motors and drive 
systems. 
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ABSTRACT 


The basic selecting peculiarities of the optimal project characteristics of the small waterplane area twin 
hull ships compared to conventional ships are considered. The description of the mathematical model and 
the ship operating model is given. The choice of the optimization method is justified. 


Keywords: small waterplane area twin hull ship; mathematical model; constraints; objective function; optimization 


INTRODUCTION 


Small waterplane area twin hull (SWATH) ships have 
excellent seaworthiness and are used as pilot, research, 
passenger, patrol, pleasure yachts due to the peculiarities of 
their hull form. Because of novelty and lack of its quantity, 
the design experience of such type of ships is little. Besides, 
an engineer has to solve a lot of problems that are not inherent 
for the traditional types of ships when choosing the project 
characteristics of SWATH ship. 

It should also be considered that in case of intense 
competition in a very short period there must be designed the 
best ship variant that will be much better than the competing 
projects. Using the traditional design methods by the original 
or the variation method gives an opportunity to get the project 
of SWATH ship that meets the task requirements but doesn’t 
guarantee its high efficiency. The way out of such situation is 
to move to an optimization design. But the application of the 
optimization approach requires special knowledge and skills 
of the engineers, especially in the design of such complex 
objects as small waterplane area twin hull ships. Therefore, 
the problem of improving the decisions quality at the initial 
designing stages of small waterplane area twin hull ships is 
rather important. 

Review of domestic and foreign publications has shown 
that there are few papers devoted to the application of the 
optimization approach for the SWATH ship design, for 
example [1-6]. The analysis of these studies gives grounds for 
the authors to conclude following issues that require further 
study: 


1. In most models the determine problem is considered that 
doesn’t allow take into account the effect of uncertainty of 
the initial information on the project efficiency. 

2. At the initial design stage the comfort requirements of 
passengers are hardly set and the reliability factor is 
ignored. 


The aim of the article is to consider the basic selecting 
peculiarities of the optimal project characteristics of the small 
waterplane area twin hull ships including the uncertainty of 
the initial information. 


BASIC MODEL 
Design problems 


The basis for the project design is a technical task (C vector) 
containing the SWATH ship specifications set by the owner. 
Vector C is as follows: 

— required service speed (kN); 

— number of passengers; 

— endurance (day); 

— hull material (steel, aluminum Alloy, glass-reinforced 
plastic); 

— superstructure material (steel, aluminum Alloy , glass- 
reinforced plastic); 

— number of struts (single, tandem); 

— type of machinery (medium speed diesel, high speed diesel, 
diesel electric, gas turbine). 
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Tab. 1. Design variables and parameters of SWATH ship 


Variable 


Description 


relative length of lower hull L,,/D, 


slenderness coefficient of strut Lyts 


waterplane area strut coefficient 


relative waterplane area Awps/V?3 


ratio of the distance between lower hull center-line to the length of the ship B,/L,, 


ratio of the ship draft to the lower hull diameter d/D,, 


ratio of the lower hull beam to its depth B,,/H, 


lower hull prismatic coefficient 


factor of the lower hull nose shape 


factor of the lower hull tail shape 


factor of the lower hull cross section shape 


strut nose and tail shape factor 


hull nose length to lower hull ratio Ly,/Ly, 


strut nose length to strut length ratio L,,/L, 


Let’s denote the vector of independent variables through 
X = (Xj, X», .-., Xa). The X vector includes the design variables 
and parameters of SWATH ship (Table 1). 

The optimization problem of SWATH ship at the initial 
stage of design is formally stated as follows: 

Minimize or maximize objective function: 


f(X, C) — min(max) (1) 
Subject to the bound constraints: 
xm <x ox j=], n (2) 


and functional constraints: 
g(X)20;1=1,...,m (3) 


where: 

m, — total number of constraints; 

xj, xia — lower and upper bounds on the independent 
variable (see Table 1); 

n — number of independent variables. 


The functional constraints of the g(X) > 0 task include 
inequalities that define the ship performance requirements. The 
following constraints include: 

- requirements for intact stability (High Speed Craft Code); 
- rolling period; 

- equality between weight and displacement; 

- minimum value of lower hull diameter; 

- maximum value of draft; 

- maximum value of breadth; 

- minimum and maximum values of strut tail length; 

- minimum and maximum values of strut tail length; 

- minimum and maximum values of lower hull tail length; 
- motion sickness indexes (MSI) and etc. 


All of these constraints are got on the basis of the analysis 
of technical requirements to the ship characteristics. There is 
a possibility to regulate the feasible search space by enabling/ 
disabling of certain constraints. 
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payload coefficient W,/A 
strut setback S,/L,, 


The objective function or criterion optimization (1) 
represents the expected value of efficiency indexes [7]. This 
criterion seeks the maximization of expected (average) profit 
or the minimization of expected operational cost: 


f(X, C) = M{EI} Prob — max 
f(X, C) = M{EI}(1 — Prob) — min 


(4) 


where: 
EI — efficiency indexes; 
M{...} — average. 


The data of the problem assumes that the payoff (or cost) 
associated with each decision alternative is probabilistic. 

Also there is a choice of another criterion form: aspiration 
level criterion, utility function. 

In this formulation the optimization problem is usually 
nonlinear and conditionally divided into two parts. The first 
part deals with the ship mathematical model development, 
the second one provides the selection of the optimal solution 
search method. 

These parts for each type of ship have their own peculiarities 
that effect the whole process of the problem solution. 
Particularly for SWATH ship it is possible to point further 
features. 

The first SWATH ship feature as an optimization object is 
the technical solutions variety used while creating in the part 
of principle project and constructive layouts and their possible 
combinations. The studies gave an opportunity to reveal the 
significant changes in the relationship nature between the 
main structural elements of SWATH ship and well-known 
monohull ship. For example, the twin hull construction means 
the increasing of the dependence of the design characteristics 
on the size and configuration of the hulls and struts, much 
changing of the external load effect nature (forces and moments 
that act in cross direction become the most important), the value 
changing of the total resistance components, and, as a result, 
the constructive measures for its decrease. 


Possible relationship variants of the hull and strut sizes, 
shape parameters are so variable that while their proving it is 
necessary to make much research work. Even small changes 
of the hull form parameters while displacement increasing 
or decreasing influence the required propulsion power and 
the weight. The same way they influence the SWATH ship 
seakeeping performance. Thus, according to the study [5], 
for the pilot ship with a small waterplane area at the stage 
of preliminary design it was necessary to make 120 steps of 
design, each of them included the calculation of resistance and 
seakeeping control. Such number of variant elaborations can be 
done only with the help of special software. Besides, the hull 
shape, got by computer-aided design, must be adapted for the 
machinery layout and rudders. 

The second SWATH ships feature is that the process of 
their optimal design is much complicated than in the traditional 
monohull ships and catamarans with simple hull configuration. 
It is known that all the performance of monohull ships are 
mainly defined by their principle dimensions (length, beam, 
draft and depth), block coefficient. To the catamarans with 
simple hull configuration the separation of demi-hulls and 
vertical clearance (Table 2) are added. And in case of SWATH 
ship it is important to consider the geometrical characteristics 
and shape parameters not only of each single hull, but also the 
struts, and, in addition, hydrodynamic interference effect of 
hulls and struts (Fig. 1) [2]. 

Besides, when designing SWATH ship it is important to 
consider the fact that external forces are defined not only by 
the geometrical hull, strut and their connections characteristics 
but also their positional relationships. 


Tab. 2. Number of design variables and combinations 


Number of 
combinations 


Ship type 


Monohull 
Catamaran 
SWATH 


35 = 243 
38 = 6561 
32 = 531441 


The next feature is the lack of sufficient design and 
construction experience, and a small amount of built ships 
which can be used as the prototypes. Nowadays approximately 
80 ships are built, and only 20 of them are for passengers. 

While the SWATH ship project design there should be 
a justification of those characteristics that are not included into 
the design task, but at the same time have strong influence on 
the SWATH ship efficiency, notably: 

- the selection of constructive and arranging type: mono, 
catamaran, trimaran, and other variants, number of struts; 

- the selection of constructive materials with an opportunity 
to combine variants that are different for the hull and 
superstructure: steel, aluminum, glass-reinforced plastic; 

- justification of the spectrum of the comfort level with the 
arranging the passengers according to the categories, as 
well as according to the decks and cabins along the length 
of the ship; 

l 


- justification of the spectrum of the cruising range and 
seakeeping levels in combination with the maximal and 
operational speed. 


The development of the general arrangement of the ship 
also needs special attention because great area of the SWATH 
ship decks provides completely new ship space structure and 
presents almost unlimited possibilities in the inhabited area 
organization. 

When solving the optimization problem there should be 
considered the existence and trustworthiness of the initial 
economical information used when estimating the ship 
efficiency. It is very difficult at the initial design stage to set 
the price of fuel, crew expenditures, port charges etc. These 
characteristics change during a season, not mentioning the 
operation time of 15-25 years. That’s why consideration of the 
economic situation instability is rational to carry out by moving 
to the stochastic formulation of the optimization problem. 

According to the information above there was developed 
a program complex (PC) SWATH Ship in order to find the best 
elements of the ship. 

The basis of the PC mathematical support consists of the 
mathematical model ofthe ship as an engineering building and 
operational model. 


Mathematical model of SWATH ship 


The SWATH ship mathematical model contains analytical 
dependences that allow to define (Fig. 2): geometrical ship 
characteristics; lightship weight and deadweight; capacity 
(required areas for passengers and areas of service, public 
and sanitary rooms); intact stability and the stability (GZ) 
curve; geometrical characteristics of the fins; ship seakeeping 
performances; building cost. 

The mathematical model is realized into two units: 
«SWATH_ model» and «Resist». 

The «SWATH model» unit contains the algorithm for 
calculation the basic SWATH ship characteristics. 

The initial data for calculation are the start values of the 
independent variables, the parameters noted in the design task 
and extra data. 

The selection of the main dimensions begins with the 
calculation of the ship payload: 


W, = (P + Nenar i Prw + P 


where: 

— one passenger mass, [kg]; 

Pew — fresh water for a passenger per day, [kg]; 
Pov — provision for a passenger per day, [kg]; 
— number of passengers; 

end ~ endurance, [day]. 


) - Np, /1000; [t] (5) 


Pas Prov. Pas 


Then the ship design displacement is estimated: 
A= Wn; [t] (6) 


where: 
Np = 


payload coefficient. 


Fig. 1. SWATH ship design variables scheme 


POLISH MARITIME RESEARCH, No 1/2013 13 


Initial data 


3 : 
' I 
i Design task i 
' I 
' I 
I t 
' I 
' t 
I I 
i i 


Economical data 


Mathematical model 


Main ship dimensions 


Initial hydrostatics 


Resistance and engine 
power 


Operation 
model 


r 
' I 
[j 
' [j 
1 [j 
' [j 
i Mission performance ' 
! probability : 
1 I 
' [j 
Jei ' i 
Weight ' Operating cost ' 
i i 
' [j 
I j 
i i 
' [j 
1 1 
' [j 
' I 
i i 
i i 
1 ! 
i i 
' ' 


Intact Stability 
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Fig. 2. Block diagram of the SWATH ship mathematical model and 
operational model 


At the next stage the calculation of the basic SWATH ship 
geometrical characteristics is performed. Some of dependences 
are listed in the Table 3. 

The formulas shown in Table 3 are obtained under the 
assumption that the lower hull forebody is elliptical, the stern 
is parabolic, the strut nose and tail are parabolic. 

The values of the geometric SWATH ship characteristics 
are used to generate the hull surface, as well as to calculate the 
propulsive performance of the ship, weight and seakeeping. 

Determination of the lightship weight and deadweight mass 
is performed by solving the weight equation: 


A=W,;+ DW (7) 
where: 
Wis — lightship weight, t; 
DW -— deadweight, t. 


At the early stage of design the calculation of the SWATH 
ship lightship weight is reasonable to perform in the following 
groups: 


Wis = Wau + Wsup t Wm t Wou t Wsm (8) 
where: 
Wau — hull weight, [t]; 
sup ~ Superstructure weight, [t]; 
Wu — machinery weight, [t]; 
Wo, — outfit weight, [t]; 
sm — design margin, [t]. 


The biggest difficulty in calculation is the hull weight 
because of the lack of information about the prototype weight 
and statistic dependences. In PC in order to define the hull 
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weight there was taken a basic approach. According to this 
approach [2] SWATH ship hull weight is estimated through the 
structural part thickness taking to account the operation load 
and requirements of the Ship-Classification Society: 


Whun = A+Co) 0+ C) Wy (9) 


where: 

C, = 0.085 — coefficient that takes into account the weight 

of additional components (painting, welding 

material and margin); 

coefficient that takes into account the weight of 

other than the plate (stiffeners); 

W; — plate weight ofthe main SWATH ship structural 
part: lower hulls, struts, sponsons, box, inside 
decks and platforms, longitudinal and transverse 
bulkheads. 


The plate weight of the SWATH ship hull structural parts 
is defined according to the following dependence: 


i=l 


(10) 


where: 
S; t, qi — area [m?], thickness [m] and density of the material 
[kg/m3] of the i hull part appropriately. 


The surface area S; of the structural parts of the SWATH 
hull is determined using the parametrical model and is directly 
connected with independent variables. 

For the ships with glass-reinforced plastic hulls the similar 
coefficients are difficult to obtain. That’s why according to 
the calculations for several small waterplane area twin hull 
ships made of glass-reinforced plastic there was obtained 
the following dependence between the hull weight and ship 
displacement: 


Won = 0.2168A + 4.6129 (11) 


The superstructure at the first approximation is calculated 
depending on the material by the formula: 


Wee = Bsir V sup? [t] (12) 
where: 
Zsu — the superstructure volume density, [t/m3]; 
Vey» — the superstructure volume with regard for the 


wheelhouse, [m3]. 


The other lightship weight groups with some changing 
and improvements are defined by formulas that are used at the 
designing of the high speed passenger catamarans. 

Providing of the passenger SWATH ship capacity is 
performed by calculating necessary areas: 


n 
As=))A,5 [m] (13) 
i=l 
where: 
A= required room area, [m7]; 
n — number of rooms on the ship. 


When calculating the cost of ship building at the initial 
stage design the following expression was used: 


Cs = (1 + ki)(Cu + Co + Cw) (14) 


where: 

k, — coefficient of the commercial expenditures and planned 
contributions, includes the value-added tax and the profit 
of the shipbuilding plant; 


Tab. 3. Calculation of the SWATH ship geometrical characteristics 


Waterplane area, [m7] 


Formula 


Ayps = V?3kw 


Fore waterplane area coefficient of the strut 


Cy._-=n/(1 +n) 


Aft waterplane area coefficient of the strut 


Cywr_a=n/( +n) 


Strut length, [m] 


L; = [Awesls 
2C wps 


Strut thickness, [m] 


ts=1/Ls 


Midship section coefficient of the lower hull 


aie 
Dh 


Cun “2 A 1) 
nh, A EN ER 
2°) -n, r(3 + a ) 


Block coefficient of the lower hull 


Cpu = CunCen 


Lower hull beam, [m] 


from the equation solution 


2C nulu Bir , Awesla-DBu yo 
by by 


Lower hull depth, [m] 


Hy = By/by 


Lower hull diameter, [m] 


D,=/HyBy 


Lower hull length 


Ly =1yDy 


Ship draft, [m] 


d=1,H, 


Midship area of the lower hull, [m7] 


Aun = CunHyBy 


One lower hull displacement volume, [m3] 


Va = CpyAmnln 


Strut submerged depth, [m] 


Hgs = d- Hy 


Strut submerged volume, [m3] 


Vs = (AwesHss)/2 


Strut setback, [m] 


S,=s _ bLy 


Length overall, [m] 


Loa = max(Lg + Sy; Ly) 


Box length, [m] 


Lyx = min(Ls + Sy Loa) 


Vertical (box) clearance, [m] 


Hpk = max(0.75h,,,; 0.625./ B,,Hy) 


Strut depth, [m] 


hs = Hss + Hpk 


Distance between lower hull center line, [m] 


Box beam, [m] 


Bs = hcLy 
Brox = Bs + By 


Depth of cross structure box, [m] 


Hopp = (Bgox — 2By)/7.5 


Depth up to the [m]ain deck, [m] 


Cu — material cost, [US $]; 
Co — equipment cost, [US $]; 
C,, — the labour cost, [US $]. 


Cost of labour is calculated as follows: 


N 
Cas +k )> cL, 
i=l 


where: 

c,, — unit hourly wage, [US $/man-hours]; 
Lp; — labour man hours [8]; 

N  — number of parts; 


(15) 


kọ — coefficient that takes into account the overhead cost. 


D = Ha + h, + Hp, 


Hull material cost: 


N 
Cu=0+k) > eW (16) 
TA 
where: ! 
k, — coefficient that takes into account the material loses; 


c; — specific is the cost of 1 ton material, [US $/t]. 
Equipment cost: 
K 
Co=> OW, (17) 
j=l 
where: 


W, — weight of unit, [t]; 


J . . 
c; — cost of outfit per unit weight or power. 
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The «Resist» unit implements the resistance calculation 
method R, and the propulsion engine power P. For defining 
the R, there was used the following formula: 


Ry = R; + Ry + Rap t Rar Ray (18) 
where: 
Rp Rap Raa — frictional, appendages and aerodynamic 
resistance; 
Rw Rsp — wave and spray resistance. 


The calculation of frictional, appendages and aerodynamic 
resistance is performed by the known dependences of the ship 
theory considering the SWATH ship construction peculiarities 
[3, 9]. 

Spray resistance calculation is based on the results of the 
model tests, as described in [10]. 

The wave resistance of the small waterplane area twin hull 
ship is defined according to the formula, kN: 


Ry =) Rwi + 22 AR wij 
where: ' _ 


> Rwi — individual wave resistance of every body that is 
i a part of the SWATH ship, [kN]; 


1 AR wij ~ additional wave resistance as a result of the 
ij wave systems interference, [kN]. The wave 
resistance is defined for those bodies such as 
underwater hull, fore and aft struts. There is also 
an opportunity to calculate the resistance for the 

Slice type ships and single-hull SWATH ship. 


(19) 


In order to calculate the components of the SWATH ship 
wave resistance there was used the Michell integral: 


4pg" ny TE o E 
Rw=—= fa +J)w’ sec? @d0 (20) 
0 

where: 

L 

2 0 
J@)=- | Pii O)cos(kyx sec 0)dxdz; 

L 
0 ° — integration variable; 


k, = g/U2 — wave number, [1/m]; 
w= 2cos(kbsec?0sin0); 


z 0 
1(0) =- | Yes Z)exp(k )z sec *9)sin (ky xsec 8)dxdz; 
Ty. = 
2b k — distance between center line of the hull, [m]; 
n — fluid density, [t/m3]; 
U — ship speed, [m/s]; 


Y(x,z) — SWATH ship lower hull or strut offsets, m. 
When calculating the wave integral, the original 
integration procedure according to the Filon rule 
was used. 


Then the main engines power is defined: 


P= Pr (SM +D, [kW] (21) 
NND 
where 
Ns — propulsive coefficient; 
SM — sea margin power; 


P,=R,U — effective power, [kW]. 
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In the program there provided the power curve output 
or the output of the resistance curve, that visually show the 
dependence of total and other types of resistance from the 
ship speed. 


Operational problems 


The operational costs and the SWATH ship economic 
efficiency indexes are defined by the economic analysis in 
operational model. Operational model allows examining the 
ship dynamic operation being effected by chance factors. The 
chance factors under the environmental effect are generally 
defined by the hydrometeorological conditions, that typical 
for the examined operational area, as well as by the initial 
uncertain data, that are used in calculation of the operational 
economic indexes. 

The operational model of the ship that makes regular 
scheduled cruises between two points is realized in the 
«Simulation» unit and contains three blocks: «Meteo», 
«Voyage», «Statistic». Before the ship setting out on a voyage 
there is a check of being ready to do that. In case of storm the 
ship voyage is canceled during the storm. Its value is generated 
in the «Meteo» block. In other cases the ship voyage between 
the departure point and destination is being modeled. While 
voyage there were calculated the coefficient of ship loading, 
the average cruise speed taking into account the wave height 
and wind speed, MSI and other operation indexes. In the port 
the passengers loading and unloading is being modeled. The 
process is repeated till the simulation time ends. Then the 
control is given to the «Statistic» block, where the statistic 
processing of the simulation results is performed and one of 
the economic efficiency indexes is calculated. 

The program provides an opportunity to simulation the 
operation of one, two and three ships on the line according to 
the following shape of service: a ship that returns and doesn’t 
return to the destination point on the same day; a ship that 
makes several trip per day; two or three ships that make series 
trip; two or three ships that make opposite trips. The voyages 
can be made every day or special days. 

The economic analysis provides the calculation of one 
of chosen ship efficiency indexes: net present value (NPV), 
required freight rates (RFR), payback period (PP), net income 
(NI) and profitability index (PI). 

The operation costs are defined as the sum of the following 
components: 

Co = Crix + Cvar (22) 
where: 
Cris — fixed cost; 
C\,, — voyage cost. 


Fixed operation cost depends on the crew number, ship 
building cost and is calculated by: 


Crix = Corew + Cr + Cig + Cp + Cor (23) 
where: 
Corew — Crew costs; 
C, — repair and maintenance; 
Ci, — insurance; 
Cp — depreciation; 
Cor — administration. 

Voyage costs include the following components: 
Cvar = Cp + Cp + Coi (24) 

where: 
C, — port charges; 
C, — fuel; 


Co; — lubrication oil. 


The SWATH ship economic efficiency analysis is performed 
with a glance of the risk factor and the factor of probability of 
the mission performance during the whole life cycle. For the 
passenger ships that perform scheduled voyages the probability 
of the mission performance can be estimated as following: 


Prob = P P,P,P, (25) 


where: 
probability of the voyage performance; 


P, — probability of keeping the given average speed during 
the trip; 

P, — economic risks (probability of the nonnegative profit 
receiving); 

P, — reliability meaning the probability of the constructions 


and equipment accident — free operation. 

The values of the probability of the mission performance 
and the parameters of the optimization criterion distribution 
law are defined via the simulation modeling of the ship trip 
elements. 


The simulation modeling [11] is based on the computer 
reproduction of the extensive ship operation process with 


Initial data 
input 
Independent variable 
values (X) 

Ship mathematical 
model 
Pseudorandom 
numbers 


Ship operational 


Output information 
accumulation: NPV, 
RFR, PP, NI, PI 


Statistical data 
manipulation 


Optimization 
algorithm 
Calculation results 
output 


Pseudorandom 


number 
generators 


a glance of external environment interaction. The SWATH 
ship operation process is presented as serial manual of the ship 
operation process between the departure and destination point, 
load/unload it in ports, etc., taking into account the hydro- and 
meteorological conditions. As a result of such modeling certain 
events and conditions are fixed that allows to define the system 
efficiency characteristics. 
The main stages of simulation modeling: 
— Accumulation and statistical data manipulation in order to 
determine the distribution law; 
— random numbers generation with given distribution laws 
using random numbers generators; 
— construction and realization of the ship operation model; 
— carrying out of the simulation experiment; 
— statistical manipulation of the modeling results. 


The block diagram of the simulation modeling procedure 
is shown in the Fig. 3. 

The second part of the problem is solved with the 
optimization method depending on the design variables 
vector length, availability constraints, non-linear criterion 
and constraints. According to the analysis results the genetic 


External action 


Fig. 3. Block diagram of the simulation modeling 
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algorithm (GA) is suggested to use while solving the problem. 
GA is a simple evolution model in natural world that is realized 
as a computer program [12-14]. In the genetic algorithm, the 
analogues of natural genetics and natural selection are used. In 
general, the GA search strategy is described by the following 
cycle. At the first iteration the initial «population» is formed (the 
whole set of the project solutions). Then for each «individual» 
(problem solutions) the fitness function values are calculated, 
that helps to identify the best «individual». After that GA 
generates anew «population» with the genetic operators of 
selection, crossover, mutation and elitism strategy. For the 
new «population» the estimation of the fitness function value, 
etc. The process is repeated till one of the stopping criteria is 
performed. 

The optimal solution search method with the application of 
genetic algorithm is realized in the «GeneticAlg» units. 

Apart from the given units that realize the SWATH ship 
design methods, computer program contains the interface 
units that provide the operational comfort for the system user. 
The computer program interface allows to input the initial 
data, to choose the objective function type and economic 
efficiency indexes, to set the initial conditions of simulation, 
the parameters of the genetic algorithm and the determined and 
stochastic economic data, to output the calculation results in 
the graphical or tabular style. 

The SWATH Ship program complex is designed in the 
Borland Delphi Professional 7.0 programing system and can 
be used in the Windows 98/XP/Vista operation systems. 

The program complex can be used in the following range 
of characteristics: passenger carrying capacity — 20...450 
persons; service speed — 20...40 knot; endurance — 100...500 
miles; ship length — 20...50 m. Besides, using the SWATH Ship 
it’s possible to carry out different experiments dealing with the 
check of the models validity, sensitivity etc. 


Results and verification 


In order to verify the calculations using the developed 
methodology and based on it programs complex, the series of 
the SWATH model towing tests (Fig. 4) were carried out in the 
towing tank of the National University of Shipbuilding. 

The towing tests of the models were carried out at speeds 
of 0.5 to 3.0 m/s, that correspond to the Froude numbers in 


the length from 0.17 to 1.07. The experiment was carried 
out in two stages. At the first stage the model with two struts 
on each hull was tested. After the first stage ended the space 
between the struts was sewn. And then the SWATH one-strut 
model was tested. 

The results of the model testing were recalculated for the 
32 m long full-size ship. Then with the help of the developed 
program the impedance values were obtained. Recalculation 
from the model to the full-size ship and the results of theoretical 


Fig. 4. SWATH models: a) single strut; b) tandem strut 


Tab. 4. Estimation of total ship resistance from model experiment and SWATH Ship program 


Ry [KN] 


from model 
experiment 


o 
Theory by program Erron ikal 


33.654 


72.071 


107.933 


122.994 


135.932 


151.919 
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Fig. 5. Comparison of measured and calculated total ship resistance: a) Single strut SWATH ship, b) Tandem strut SWATH ship 


Tab. 5. The SWATH optimal performance for the Odessa — Varna route 


Hull/Superstructure material 


Description 


Lower hull length, [m] 


Lower hull beam, [m] 
Lower hull depth, [m] 
Hull nose length, [m] 
Hull tail length, [m] 
Strut length, [m] 
Strut thickness, [m] 
Strut height, [m] 


Aluminum Alloy Steel 


Bier Rive /Aluminum Alloy | /Aluminum Alloy 


Strut nose length, [m] 


Strut tail length, [m] 


Waterplane area strut coefficient 


Box clearance, [m] 


Distance between lower hull center line, [m] 
Ship draft, [m] 
Depth up to the [mJain deck, [m] 
Length overall, [m] 26.378 25.123 25.808 
Box length, [m] 26.378 25.123 25.808 
Box beam, [m] 12.979 11.413 12.623 
Depth of cross structure box, [m] 1.004 0.98 1.004 
Displacement, [t] 250 150 225 
Deadweight, [t] 39.26 34.75 37.95 
Main Engines, [number x kW] 2 x 3460 2 x 2300 2 x 3460 
Generator, [kW] 
Crew 
Cost of ship, thousand [US $] 
Payback period, [year] 
Net Present Value, thousand [US $] 
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calculation of the one-strut SWATH towing resistance are shorn 
in the Table 4. Similar results were obtained for the two-strut 
SWATH. 

According to the tests and calculated data using the 
developed program, the diagrams of the towing resistance and 
Floude number dependence for the one- (Fig. 5a) and two- strut 
(Fig. 5b) SWATH were constructed. 

The results, obtained with the help of theoretical calculation 
using the developed program, quite accurately match the test 
data of the small size models. 

The technique’s working efficiency is shown by the example 
of solving the problem of choosing the best performance of 
the passenger SWATH for the Odessa — Varna route. The 
time schedule of the Krymskaya Strela catamaran is used in 
the calculation. During the calculating of capital investment 
in the ship construction it is supposed that the buyer’s own 
funds are 20%, and the rest 80% of investment is the bank 
loan for 8 years under the 6...10 % interest rate per year. The 
ship operational lifetime is 15 years. The net present value is 
used as an economical efficiency factor. The SWATH optimal 
performance values for the passenger Odessa — Varna shipping, 
obtained after the work of optimization program, are listed in 
the Table 5. 

In order to define the main SWATH characteristics the 
genetic algorithm with the following parameters was used as an 
optimization method: population number — 50 chromosomes, 
gene capacity — 32 bit, crossover probability — 0.9, mutation 
probability — 0.1, inversion probability — 0.05, initial penalty 
— 0.5, extreme achieving accuracy — 0.000001. The elitism 
strategy was used during the optimization. These parameters 
are set experimentally as a result of multiple test runs of the 
program. 


| tL 


The calculation results have shown that the most 
economically efficient SWATH model is the one that is made 
of aluminum because it brings the highest return at the lowest 
expenditures and has less payback period. 


CONCLUSIONS 


1. The SWATH ship is more complicated for the optimization 
research than the conventional monohull ships and 
catamarans with traditional hull shape. First of all it’s 
connected with the variability of the used technical solutions 
in the project and construction arrangement and their 
possible combinations. It’s also connected with the difficult 
optimization process and lack of design and construction 
experience. 

2. The general problem statement of the SWATH ship optimal 
design is characterized such complexity factors as large 
number of independent variables, presence of constraints, 
necessity to account the stochastic and uncertainty external 
agencies. The solving process of such problem provides the 
use of the penalty function approach (for the constraints 
account), genetic algorithm (for the direct optimum search) 
and simulation modeling (for the accounting of the data 
uncertainty). 


Further research work is advisable to direct for improvement 
of the algorithm calculation of the propulsive coefficient, 
seakeeping performance and for enlargement of the model for 
other types of small waterplane area twin hull ship. 


Appendix A 


Fig. A.1. Hull form parameters 
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midship area of the lower hull, [m?] 
required room area, [m?] 

waterplane area, [m?] 

ratio of the lower hull beam to its depth 
breadth of the ship, [m] 

box beam, [m] 

lower hull beam, [m] 

distance between lower hull center line, [m] 
vector of technical task 

block coefficient 

block coefficient of the lower hull 

crew costs, [US $] 

depreciation cost, [US $] 

fuel cost, [US $] 

fixed cost, [US $] 

coefficient that takes into account the weight of 
other than the plate 

insurance cost, [US $] 

material cost, [US $] 

midship section coefficient of the lower hull 
equipment cost, [US $] 

administration cost, [US $] 

lubrication oil cost, [US $] 

cost of port charges, [US $] 

lower hull prismatic coefficient 

repair and maintenance cost, [US $] 

cost of ship, [US $] 

voyage cost, [US $] 

the cost of the work of the shipbuilding plant, 
[US $] 

aft waterplane area coefficient of the strut 
fore waterplane area coefficient of the strut 
waterplane area strut coefficient 

coefficient that takes into account the weight 
of additional components (painting, welding 
material and margin) 

ship draft, [m] 

depth up to the main deck, [m] 

lower hull diameter, [m] 

deadweight, [t] 

efficiency criterion 

Froude number 

ship performance requirements 

the superstructure volume density, [t/m3] 
ratio of the distance between lower hull center- 
line to the length of the ship 

strut depth, [m] 

wave height of 3% probability 

depth of cross structure box, [m] 

vertical (box) clearance, [m] 

lower hull depth, [m] 

strut submerged depth, [m] 

relative waterplane area 

wave number, [1/m] 

coefficient of the commercial expenditures and 
planned contributions, includes the value-added 
tax and the profit of the shipbuilding plant 
ratio of the ship draft to the lower hull diameter 
relative length of lower hull 

hull nose length to lower hull ratio 

strut nose length to strut length ratio 
slenderness coefficient of strut 

length of the ship, [m] 

length overall, [m] 

box length, [m] 

lower hull length, [m] 

strut length, [m] 

motion sickness indexes 

factor of the lower hull tail shape 

factor of the lower hull nose shape 

factor of the lower hull cross section shape 
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strut nose and tail shape factor 

number of passengers 

endurance, [day] 

engine power, [kW] 

effective power, [kW] 

fresh water for a passenger per day, [kg] 
one passenger mass, [kg] 

probability of the mission performance 
provision for a passenger per day, [kg] 
probability of the voyage performance 
probability of keeping the given average speed 
during the trip 

economic risks (probability of the nonnegative 
profit receiving) 

reliability [m]eaning the probability of the 
constructions and equipment accident-free 
operation 

density of the [m]aterial hull part [kg/m3] 
aerodynamic resistance, [kN] 

appendages resistance, [kN] 

frictional resistance, [kN] 

spray resistance, [kN] 

total resistance, [kN] 

wave resistance, [kN] 

strut setback 

strut setback, [m] 

area of the i hull part [m?] 

sea margin power 

thickness of hull part [m] 

strut thickness, [m] 

ship speed, [m/s] 

model speed, [m/s] 

superstructure volume, [m3] 

hull weight, [t] 

lightship weight, [t] 

outfit weight, [t] 

payload, [t] 

machinery weight, [t] 

plate weight of the main SWATH ship structural 
part: lower hulls, struts, sponsons, box, inside 
decks and platforms, longitudinal and transverse 
bulkheads 

design margin, [t] 

superstructure weight, [t] 

minimun value of independent variable 
maximum value of independent variable 
vector of independent variables 

SWATH ship lower hull or strut offsets, [m] 
distance between center line of the hull, [m] 
gamma-function 

displacement, [t] 

payload coefficient 

propulsive coefficient 

fluid density, [t/m3] 


individual wave resistance of every body that is 
a part of the SWATH ship, [kN] 


additional wave resistance as a result of the 
wave systems interference, [kN] 


volume displacement, [m3] 
one lower hull displacement volume, [m3] 
strut submerged volume, [m3] 
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Analysing the potential for application 
of the phase shift method in endoscopic 
examination of marine engines 


Zbigniew Korczewski, Prof. 
Gdansk University of Technology, Poland 


ABSTRACT 


The article presents theoretical principles of image processing in digital endoscopy 
which makes use of a miniaturised spectral scanner “PhaseProbe” designed by General 
Inspection Technologies LP, the producer of the measuring videoendoscope Everest XLG3. 
The technology of optoelectronic 3D mapping of the examined surface, which consists in 
measuring the phase shift between the emitted and reflected light waves, is briefly described. 
The efficiency of the ““PhaseProbe” based method of endoscope measurements is compared 
with the earlier developed methods, such as “ShadowProbe” and “LaserDots”’, which 


were used by the author in real conditions of diagnostic investigations of engines in operation in marine 

power plants. The processed results of the endoscope analysis of the broken gas turbine engine rotor 

blades which were earlier dimensioned using the inspection probes of “StereoProbe”, “ShadowProbe”’ 

and “LaserDots” type and now are examined by the author with the aid of the phase shift method make 
a valuable complement and synthesis of the discussion presented in the article. 


Key words: endoscopic diagnostics; phase shift method; marine engines 


INTRODUCTION 


Despite rapid technological development in measuring 
methods used in industrial endoscopy, precise measurement 
of the dimensions of the detected surface defects is still 
a challenging problem. New opportunities in this area have 
been brought by various 3D scanning techniques which consist 
in converting the real three-dimensional shape of the examined 
surface into a digital form. The theoretical principles of the 
phenomena of light wave diffraction and interference which 
accompany this mapping were worked out over 200 years ago 
by Thomas Young), but only as recently as in 2010 the most 
modern optoelectronic devices for light emission and detection 
made it possible to design the measuring videoendescope 
Everest XLG3 by General Inspection Technologies LP. This 
videoendoscope makes use of the spectral scanner miniaturised 
to the overall dimension of the inspection probe head of 
an order of 6 mm [4, 5, 6]. Moreover, a method has been 
developed which uses the results of the light wave phase shift 


measurement”) to create a digital 3D model of the examined 
surface being a precise mapping of its real geometric shape 
and dimensions?). 


PHASE SHIFT METHOD - 3D 
MEASUREMENT 


A key element in the design of the inspection probe 
of the measuring videoendoscope Everest XLG3 in the 
“PhaseProbe” option is the diffraction lens mounted in the 
inspection probe head — Fig. 1. The design of the lens bases 
on a pair of symmetrically fixed diffraction gratings (2) 
which generate a static interference pattern having the form 
of halftone screens with known parameters, (i.e. gratings of 
dark and bright interference fringes with certain density) on 
the examined surface. The part of the halftone screen which 
is within the visual field of the image processing optics (1) 
is recorded on the matrix of the CCD camera as a cloud of 
points, the X, Y, Z coordinates of which are known. Then this 


D Thomas Young (1773-1829) — an English physicist and physiologist. Initiated the wave theory of light by stating that the light waves are 
transverse waves. In 1801 he discovered diffraction and interference of light waves (Young’s experiment). A detailed description of these 


phenomena was issued by him in 1807. 
2) United States Patent No.: US 7,821,649 B2. 


3) Inthe promotional materials issued by Everest VIT in Polish, the light wave phase shift measurement method is more briefly referred to 


as the “3D phase method”. 
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pattern is analysed and converted, in the triangulation process 
done by the central unit of the videoendoscope, to a grid of 
triangles which transform the virtual image of the examined 
surface into its 3D equivalent. A complementary element of the 
“PhaseProbe” lens is the window of the standard illumination 
system fed with a 75-watt high-intensity discharge arc lamp 
(HID) having the luminous flux equal to 4300 lumen, as the 
measure of its brightness [3]. The resultant colour temperature 
of the light source?) is approximately equal to 5000 K. 
a) b) 


Fig. 1. Diffraction lens for measuring the phase shift with the aid of 
“PhaseProbe” of the videoendoscope Everest XLG3 [4]. a) straight-ahead 
direction of observation (field of vision angle “FOV” — 105°, depth of 
field “DOF” - 8-250 mm, b) side direction of observation (“FOV” - 105°, 
“DOF” - 7-250 mm). 1) image processing optics — CCD camera, 

2) diffraction grating, 3) window of standard illumination system. 


The standard illumination system of the videoendoscope 
works in the continuous mode only during routine observation 
of the examined surfaces. During the measurement the standard 
illumination system is automatically switched off to reach the 
maximal possible resolution of the interference fringe pattern 
generated by the diffraction gratings, which are illuminated in 
this time by electro-luminescence diodes (LED) radiating the 
precisely defined light wave length. 

When the distance of the “PhaseProbe” lens from the 
examined surface changes, the halftone screens of diffraction 
fringes also undergo relevant changes, according to the schematic 
diagram shown in Fig. 2. Therefore based on the records and 
analyses of the deformation pattern of particular diffraction 
fringes we can conclude about qualitative and quantitative 
surface deformations of the examined object. Since the observed 
dimensions of the surface defects are not only the function of 
their real dimensions but also of the distance of the inspection 
head lens from the examined surface, its precise determination 
is a basic metrological problem in 3D measurements. For this 
purpose we can apply the phase shift method, initially described 
by Thomas Young — Fig. 3. Monochromatic, incoherent light 
waves (revealing the same length but different and time- 
dependent phase shifts) are emitted in a sequence by the light 
sources LED1, LED2, LED3. These light waves undergo double 
diffraction: initially on the gaps which equalize the vibration 
phases, and then, as the monochromatic coherent waves, on the 
diffraction grating gaps, thus creating a pattern of bright and 
dark interference fringes of the same intensity on the examined 
body surface. The bright fringes (maxima of interference) are 


gap 1A which equalizes the 


vibration phase 
LED — emission 
group A 


} 


image mapping optics 
— CDD camera 


3D head distance Di interference fringes 


dis Te D2 


Fig. 2. Schematic diagram of the method of 3D measurement with the aid 
of the “PhaseProbe” inspection probe diffraction lens 


situated in places of highest superposition of the two interfering 
light waves of the same length, while the dark fringes (minima 
of interference) — in the place of their highest difference. We 
can talk about so-called constructive interference, for which 
the condition of the equality of the light wave optical path 
length difference (phase difference) and a multiple of the wave 
length is met, or destructive interference — when the optical 
path length difference of the light waves is the odd multiple 
of the half of the light wave length. Using the nomenclature 
for macroscopic quantities shown in Fig. 3, and taking into 
account that the distance D of the inspection head lens from the 
examined surface is much longer than the diffraction grating 
constant d (D>>d) we can assume that the line segment AC is 
almost perpendicular to the optical paths r, and r,. That means 
that the angle ABC is almost equal to the angle PRO (denoted 
as a). Hence the optical path length difference between the 
interfering light waves (being the effect of the phase shift) can 
be given by the following formula: 


BC =Ar=d-: sina (1) 


Thus the condition for the appearance of the bright 
interference fringe on the examined surface is identical with 
the constructive interference condition: 


d-sina=n-A (2) 


where: 
n=1,2,3... - is the order of the successive interference 


fringe. 


For point P in Fig. 3, the phase difference of the emitted 
light waves is equal to their four lengths. 

The condition for the appearance of the static dark fringe in 
interference pattern is identical with the destructive interference 
condition: 


d: sina = (2n + 1)-A/2 (3) 


where: 
n=1,2,3... - is the order of the successive interference 


fringe. 


4 The colour temperature is the method of visible light evaluation used in illumination, photography and filming. The colour temperature 
of the source of light is assessed by comparing its chromaticity to that of the perfectly black body. The colour temperature is usually 
measured in Kelvin degrees, and it is the temperature to which the perfectly black body is to be heated to obtain the given colour. Contrary 
to a common sense, a higher colour temperature (above 5000K), means “cool” (green-blue) colours, while a lower colour temperature 
(2700-3000K) means “warm” (yellow-red) colours. The “cool” light is considered better in case of tasks which need concentration [www. 


ModernHome.pl]. 
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diffraction grating 


gaps which equlize 
vibration phase 


examined surface 


interference fringes 


U3 


N 


Fig. 3. Method of calculating the distance between the inspection head and the examined surface in the phase shift method. S, S, S; — gaps which 
equalize the vibration phase, Sp, Sp>....Sp, — diffraction grating gaps (for clarity of the picture, their structure is presented in part only and is over- 
enlarged), r, r, — optical paths of the light waves, d — diffraction grating constant (distance between two neighbouring gaps, diaphragm width), 

D — distance of the inspection head lens from the examined surface, y — distance of the interference fringe of a given order from the zero-order fringe, 
0, 1, 2, 3, 4 ... n — order of the successive interference fringe 


If the dark interference fringe was expected to appear at 
point P in Fig. 4, the phase difference between the emitted light 
waves should be equal to the odd multiple of the half of the 
wave length, i.e. four and half of the wave length. 

Since the light beam deflection angle a at which the bright 
fringes are created increases for each successive interference 
number, the value of the deflection angle sine function increases 
as well, thus determining the distance between particular fringes 
in the interference pattern. The first-order fringe is always 
situated opposite the diaphragm, between the gaps d and is the 
symmetry axis of the interference pattern. 

If we manage to create a stable and sharp pattern of widely 
distributed interference fringes on the examined surface as 
aresult of light wave diffraction and interference, we can 
measure with high precision their distances from the zero-order 
fringe and then, having known the length à of the emitted 
light wave and the diffraction grating constant d, calculate the 
distance of the videoendoscope inspection head lens from the 
examined surface D using the following formula: 


y =n-aA 


d- 


Hence: 


(5) 


Of high importance, from the technological point of 
view, in the production of the probe lens of the “PhaseProbe” 
videoendoscope is precise incision of two diffraction gratings, 
consisting of parallel and equidistant grooves having minimal 
dimensions of an order of 1-2 mm, on two glass plates using 
a diamond cutting edge. These grooves compose the diffraction 
grating diaphragms which separate its optical gaps. The 
width of the gap is to be approximately equal to the width 
of the groove and, simultaneously, to the light wave length. 


The sine qua non condition for obtaining a stable pattern 
of the interference fringes from the diffraction grating is 
securing a monochromatic and coherent light sources). The 
light waves emitted by the LED diodes are monochromatic, 
which means that they have the same wave length A, but they 
are not coherent, i.e. the sinusoidal waves of the light beam 
do not vibrate with the same amplitude, frequency and phase, 
or time-independent phase difference. To solve this problem, 
the incoherent source of LED light generates the waves on 
the diffraction gratings of the “PhaseProbe” head lens via 
additional gaps which equalize their phases (according to the 
Huygens principle®, the waves which leave the diffraction 
grating gaps are generated by vibrating points situated on the 
same wave front, which means the presence of a constant phase 
difference between them). 

To provide opportunities for complete and accurate 3D 
mapping of the examined surface, the phase shift method 
makes use of a system of two diffraction gratings, each o 
which is fed from three LED emission groups which emit, in 
a sequence, sinusoidal light waves mutually shifted by a phase 
angle equal to 120° (this solution secures high sharpness of the 
static interference pattern) — Fig. 4. 

This way a rapidly changing projection and recording 
of two patterns (half-tone screens) of interference fringes 
with known density take place on the examined surface. 
These fringes are deformed adequately to the scale of surface 
deformation (they deflect left when the material has piled up 
and right in case of material decrement). The pattern of the 
light reflected from the examined surface is recorded on the 
CCD camera matrix. 

The videoendoscope processor system (CPU) identifies 
and localises particular fringes. Then it performs a comparison 
analysis of the phase shift between the known emitted 
light and the reflected light. The results of the calculations, 
performed taking additionally into account parameters of 
camera calibration and the angle between the direction of the 
projection and that of the halftone screen recording, have the 


5 This condition can be only met by laser light sources. 


© Christiaan Huygens (1629 -1695) — Dutch mathematician, physicist and astronomer. In the area of physics, he performed a number of 
works on optics. He formulated the Huygens principle, which says about the wave propagation in the vicinity of obstacles. His experiments 
confirmed the theory on the wave nature of light. He also examined the phenomenon of light refraction in various materials. 
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form of X, Y, Z coordinates for each camera pixel. The cloud 
of points created in the above manner, the number of which 
results from the resolution of the applied camera, makes the 
basis for dimensioning the detected surface defects. 
a) gaps LA,2A3A 
LED - — 
emission 
group A 


diffraction 
grating A 


diffraction 
grating B 


LED - 
emission 
group B 


Fig. 4. Phase shift in the 3D measurement method [6]. a) interference 
pattern generated by the inspection head lens, b) phase shift pattern of the 
sinusoidal light waves emitted by three emitting groups LED 1, 2 and 3 and 

the diffraction gratings ,,A” and „B”. 


ADVANTAGES AND DISADVANTAGES OF 
THE 3D MEASUREMENT METHOD 


The videoendoscope Everest XLG3 is equipped with a 3D 
measurement probe bearing the name of “PhaseProbe” which 
is a „flexible three-dimensional eye” of the operator. These 
borescopes reveal numerous advantages, which are the reason 
why they are in more and more frequent use. The basic advantage 
of the 3D measurement technology with respect to older 
measuring methods, such as “StereoProbe”, “ShadowProbe”, 
and “LaserDots” for instance, is its ergonomics. Switching from 
the standard (qualitative) measurement mode to the measuring 
(quantitative) mode does not require withdrawing the inspection 
probe from the inside of the examined object and changing the 
optic lens. Apart from complicated manual work to be done to 
change the lens, a big problem in those cases was finding again 
the earlier detected surface defect, the more so that the measuring 
lens reveal, as arule, low quality of image processing of the 
examined surface observed from a larger distance. For instance, 
the depth of field for “ShadowProbe” lens is within 7-30 mm, 
while for “PhaseProbe” it ranges between 7-250 mm! 

On the other hand, “PhaseProbe” is very sensitive to 
image movements (vibration) and light reflection from highly 
reflective surfaces. The measurement cannot be performed, in 
practice, when the operator does not manage to keep the tip of 
a flexible and a number of meters long probe still and situated at 
a proper angle for at least two seconds. It is extremely difficult 
to reach this state, even if a special probe rigidiser with a tube 
gripper is used for this purpose [2, 3]. From the practical point 
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of view it is a serious disadvantage of the “PhaseProbe” based 
3D measurement method and it should be eliminated in the 
future by its inventors. 

Taking into account high cost of purchase of the measuring 
set, which, despite the decreasing trend, still remains at the 
level of 40-50 thousand Euro, each time a decision about its 
purchase should have strong rational background. 

When discussing possible areas of application of the 
measuring videoendoscope with “PhaseProbe” in marine 
engine diagnostics, working spaces should be mentioned for 
which surface wear is decisive for the efficiency of the energy 
conversion processes taking place in the engine. Tasks to be 
done in those cases include precise determination of parameters 
characterising the roughness of the surfaces composing the 
borders of the working space. This should be done over 
a relatively large area, which from the point of view of 
endoscope diagnostics requires its full numerical mapping and 
dimensioning. In those cases the quantitative measure of wear 
of the examined surface is the arithmetic mean deviation of its 
profile from the average line determined along a normalised 
elementary line segment, or the roughness height calculated 
using ten points of this profile. 

The research experience gained in the past by the author 
suggests that the phase shift method can be efficiently used for 
diagnosing the following constructional elements: 

1. For working spaces of piston engines: 

a) cylinder bearing surface — honing grooves, 

b) valve seats — valve set faces, wear thresholds on valve 

heads, 

c) air and exhaust gas flow ducts — shape and geometric 
dimensions, active flow sections, condition of inner 
surfaces, 

d) turbo compressor rotor assembly - shape and geometric 
dimensions of blades, condition of blade profile surfaces 
and inter-blade passages; 

2. For flow parts of turbine engines: 

a) rotor assemblies of compressors and turbines — shape 
and geometric dimensions of stator and rotor blades, 
condition of blade profile surfaces and inter-blade 
passages, 

b) combustion chamber - shape and geometric dimensions 
of blade systems in flame tube swirl vanes, condition of 
inner and outer surfaces of flame tubes. 


The measuring methods of “Stereo”, “Shadow” and/or 
“LaserDots” type provide opportunities for digital processing 
of stereoscope effects, which makes it possible to dimension 
the observed objects in such a way that they give an impression 
of quasi three-dimensionality with its depth, massiveness and 
mutual distribution. Unfortunately, their application to full 
(sufficiently detailed), three-dimensional mapping of larger 
surfaces is very limited for the following reasons: 

1. The “Stereo” method: the 3D measurement bases on 
precise correlation of two identical characteristic points of 
the left and right view of the dispersed stereoscope image [4]. 
A problem can arise when no clear characteristic points can 
be found on the examined surface. The second limitation is 
remarkable labour consumption of the measurement which 
results from small angles of the visual field and depths of field 
of the used “Stereo” type lens, which are equal to: ,,FOV” 
— 45+60°, ,,DOF” - 2+80 mm. This prolongs the time of the 
measurements necessary for full three-dimensional mapping 
of the entire examined surface. For these reasons the “Stereo” 
method is more efficient in cases when only a limited number 
of surface points is required for measuring basic dimensions 
of the detected defects. 


2. The “LaserDots” method: the 3D measurement consists in 
projecting a matrix of markers composed of a number (49, 
for instance) of laser beam points on the examined surface 
and numerical recording of their shift to the left or right of 
the screen as the result of surface profile changes [2]. Since 
the areas situated between the matrix markers cannot be 
interpolated in this way, any surface changes taking place in 
these areas cannot be dimensioned either. As a consequence, 
the 3D mapping of the examined surface is only limited 
to laser marker points, which does not secure the required 
resolution of the measurement. 

3. The “Shadow” method has a form of single-fringe scanning 
and brings useful three-dimensional information only with 
respect to one plane of the examined surface profile within 
a very limited area, as a result of optical limitations of 
the applied “Shadow” type lens (,,FOV” — 50°, ,,DOF” - 
- 7+30 mm). An additional difficulty in the realisation of 
the 3D measurement with the aid of the “Shadow” method 
is the need for very precise linear positioning of the fringe 
with respect to the examined surface, which requires some 
experience. If the surface is not ideally flat, or the inspection 
probe lens is not directed perpendicularly, then a relatively 
large error can be recorded in the measurements which 
require referring the baseline to the points situated off the 
line, the line-point distance measurement for instance. 
For these reasons, and bearing in mind the fact that, like 
in the “Stereo” method, a huge-capacity processor is to 
be installed in the videoendoscope to solve equations of 
the mathematical model describing the 3D profile of the 
examined surface in this technique, the application of the 
“Shadow” method is rather limited to only dimensioning 
the detected surface defects [2]. 


3D MEASUREMENT TECHNOLOGY 


The endoscopic examination of inner spaces of the engine 
using the “PhaseProbe” of the videoendoscope Everest 
XLG3” can be realised in the standard mode or the measuring 
mode without change of inspection head lens. The standard 
examination is oriented on general qualitative evaluation of 
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Fig. 5. View of the computer screen after starting the programme GE Inspection Manager v.06 in the “Measurement area” mode 
(image parts which are not accessible for measurement as excessively distant, dark or screened by light reflection are marked red) 
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the technical state of a constructional element in the engine, 
along with the search for and identification of failures of 
surface defects. If any change in the constructional structure 
of the examined elements is identified, an attempt is to be 
made to evaluate it in quantitative way using the measuring 
mode for this purpose. As already mentioned, it does not 
require any additional actions concerning the operation of the 
inspection probe, except moving its lens as close as possible 
to the examined surface, while preserving at the same time 
a sufficiently large depth of field (like in the measurements 
performed in the “ShadowProbe” and “StereoProbe” options) 
and an optimal visual angle (45° is the recommended angle 
in case of the examination of shining or spotty surfaces). 
A remarkable advantage of the software of the videoendoscope 
equipped with “PhaseProbe” is continuous analysis and 
adjustment of the distance between the measuring lens and the 
examined surface. During the measurement of the phase shift, 
a range-finder strip is displayed on the LCD monitor of the 
video-probe handset (1 in Fig. 6a). Exceeding the measuring 
range is signalled by a single red strip, which means that we 
should not make attempt to measure. When the distance of 
the lens from the examined surface decreases, it results in 
the increasing number of the displayed strips, which, what 
is more, change the colour to green, which means that the 
measurement can be done. Scanning, started by pressing 
the push-button Enter on the manual panel, initiates the 
process of fast (lasting about 1-2 seconds) and multiple 
recording of the digital image. This is the crucial action of 
the technology of 3D measurement (as mentioned above), 
which lasts until the message “Processing” is displayed on 
the monitor screen. After the measurement is completed, the 
distance of the lens from the examined surface (at the point 
of the active cursor) is displayed on the monitor as the MTD 
(Maximum Target Distance) index. Unlike the “StereoProbe” 
and ”’ShadowProbe” methods, the phase shift method does not 
require additional correction of the measurement depending 
on this distance (“Accuracy Index”). Further processing of 
the recorded image is done using a specialist software, GE 
Inspection Manager v.06*), see Fig. 5, after completing the 
endoscope examination. 


7” Due to volume limitations of the present article, it does not include detailed constructional and parametric characteristics of the 
videoendoscope Everest XLG3. This information can be found on internet pages of General Inspection Technologies LP and dealers of 
Everest VIT endoscope equipment. Moreover, detailed descriptions of the methodology of the endoscope based diagnostic investigations 
of marine engines, which are also omitted here, can be found in numerous earlier publications by the author. 


8) By courtesy of “Everest VIT GmbH Poland”. 
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In the “PhaseProbe” method the following measuring 
options are available: 
- length, 
- multi-segment length, length of a broken line (perimeter), 
- distance of a point from the straight basic line, 
- depth (convexity), 

diameter of the area marked by a circular rule, 


Area=2.44 


- profile and profile map of the examined surface, 
as well as the 3D view of the cloud of the measuring points 
which map the examined surface, with options of its rotation 
with respect to X, Y and Z axis. 

Figures 6 and 7 show the processed results of the 
endoscope examination of the turbine engine rotor blades using 


Everest VIT 


Fig. 6. Results of the 3D measurement of the indentation on the turbine engine rotor blade trailing edge performed with the aid of the phase shift method. 
a) 1D measurement of indentation length — 3.03 mm, at MTD index — 16.79 mm, b) measurement of indentation profile — 0.98 mm, at MTD index — 16.79 mm, 
c) indentation profile sector — 0.98 mm, at MTD index — 16.79 mm, distance between cursors Z=7.07 mm, d) 3D view of the cloud of measuring points, 

e) 3D view of the cloud of points with the reference plane, f) dimensioned map of surface deformation profile 
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Everest VIT 


Area=0.000 


Fig. 7. Results of the 3D measurement of local material deposit (“eczema”) in the turbine engine rotor blade trough performed with the aid of the 
phase shift method. a) 1D measurement of deposit diameter — 0.86 mm, at MTD index — 10.01 mm, b) measurement of deposit profile — 0.14 mm, 
at MTD index — 10.40 mm, c) deposit profile sector, distance between cursors Z=1.6 mm, d) dimensioned map of surface deformation profile 


selected for 3D processing and dimensional analysis of the 
recorded numerical images had been earlier examined with the 
aid of the measuring probes “StereoProbe”, “ShadowProbe” 
and “LaserDots” [2]. Figures 6a and 7a show the results 
of standard dimensioning of, respectively, the length of 
indentation and point material deposit on the examines turbine 
blade. The examination was performed at rather large values 
of the MTD index: 16,79 mm and 10,01 mm. These results are 
very similar to those obtained in the past using other measuring 
methods. But completely new metrological capabilities have 
been brought by the 2D measurement of surface deformation 
profile, Figs. 6b, 6c, 7b and 7c, the three-dimensional cloud 
of measuring points, Figs 6d and 6e, and the dimensional map 
of surface deformation profile, Figs. 6f and 7d. A possibility 
to rotate the digital model of the examined surface which 
precisely models its real shape and geometrical dimensions 
makes the process of diagnostic inference easier and more 
efficient, and remarkably increases the reliability of the 
formulated diagnosis. 


FINAL REMARKS AND CONCLUSIONS 


A basic condition for formulating a reliable endoscope 
diagnosis of the technical state of working surfaces in 
a marine engine is an opportunity to perform not only 
qualitative but also quantitative assessment of the detected 


surface defects. New perspectives in this area are brought by 
digital endoscopy. Numerical image analysers cooperating 
with the measuring heads of “StereoProbe, “ShadowProbe”, 
“LaserDots” and here described “PhaseProbe” type provide 
opportunities for numerical processing of stereoscopic effects 
which in turn makes it possible to dimension the observed 
images in such a way that they give an impression of quasi 
three-dimensionality, with its depth, massiveness, and mutual 
distribution. 

Recent years show that further development of the 
endoscopic diagnostics of marine engines is absolutely 
determined by increasing technical capabilities of the more 
and more perfect and faultless measuring equipment. Precise 
and ergonomically designed digital endoscopes open new 
prospects for developing the diagnostic knowledge on the 
kinetics of working space wear in engines in operation, 
leaving less and less space for traditionally used optical 
endoscopes [1]. However, there is one sine gua non condition: 
the diagnostician should study in detail technical capabilities 
of the owned measuring videoendoscope and be able to use 
them effectively in the environmental conditions in which the 
diagnostic examination is performed. He also should avoid 
operating mistakes which would lead not only to the decrease 
of the life time and reliability of the very expensive measuring 
endoscope equipment, but also to worsened reliability of the 
diagnosis of the examined engine. 
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turbogenerators with the assistance 
of vibration diagnostics 
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ABSTRACT 


The paper discusses the problem of maintenance of the marine steam turbogenerator sets used in power 
systems on VLCC type oil tankers. The analysis of operation, monitoring and applied diagnostic methods 
of a selected turbogenerator working in the AT series waste energy system produced by Mitsubishi is 
presented. The analysis focuses on the applied maintenance method and its credibility and effectiveness in 
terms of technical assessment. The results of the performed maintenance actions and applied operational 
procedures making use of technical diagnostics and supported by a wider range of measurement capabilities 
of external service are presented and evaluated. The range of the study includes the readiness, reliability 
and safety of the use of turbogenerators being the result of the maintenance performed in accordance with 
the applied strategy of operation. 


Keywords: turbogenerator; vibration diagnostic; service; maintenance; condition monitoring; signal analysis 


INTRODUCTION 


Steam turbines in auxiliary power units of ship power 
systems usually perform the function of power generators 
used during a sea-passage as the main source of electrical 
energy and cargo pump drive. They are mainly used on vessels 
specialising in handling, storage and transport of crude oil, such 
as FPSO (Floating Production, Storage and Offloading) units 
for instance, and in waste energy recovery systems for VLCC 
(Very Large Crude Oil Carrier) tankers. They are produced by 
such manufactures as, inter alia: Mitsubishi Heavy Industries 
Ltd, Peter Brotherhood Ltd., Shinko Ind. Ltd [2]. 

During the operation of the marine steam turbogenerators, 
the processes of their intensive use resulting in the production 
of the required electric power take place, accompanied by 
directly related processes of maintenance and monitoring 
of technical state and suitability to perform tasks, as well 
as functional, technical and operational readiness, while 
preserving technical and economic efficiency under given 
restrictions and disturbances. These processes include: 

-  identification/determination and preservation of the required 
technical state, 
- technical/maintenance, prevention, and control services, 

complemented by periodical repairs [5, 9]. 


The above defined maintenance of steam turbines in marine 
power generation sets includes control methods, measures and 
systems used for controlling their functional and task-oriented 


readiness for rational and efficient use at an expected level in 
given conditions and time, according to their specifications. 
Here, of high importance is the diagnostics performed with 
the use of integrated and stand-alone systems and methods of 
technical diagnostics. It is performed in direct contact with the 
turbine and plays the superior role in maintenance, with the 
identification and analysis of its results allowing to control the 
turbine operation process. 

The maintenance strategy is adopted taking into account 
the existing conditions and resource supplies. For tankers, it is 
most common to use the planned maintenance strategy (PMS 
— Planned Maintenance System) which bases on the time- 
schedule (TBO — Time Base Overhauling) making use of the 
preventive maintenance (PM — Preventive Maintenance) and 
elements of technical diagnostics. At present, monitoring the 
correctness of turbogenerator operation consists in disassembly- 
free control of working processes, done via the evaluation of the 
generated diagnostic signals, including selected representative 
vibration signals, and their comparison with the reference 
values. The scope and credibility of the control depend largely 
on the used measurement equipment and qualifications of the 
personnel [1]. 

The aim of the article is to validate the proposed method 
of steam turbine maintenance with the modification of the 
set of parameters needed for monitoring the correctness of 
turbogenerator operation, and justify the need for monitoring 
the trend of changes of the parameters carrying the highest 
information capacity. 
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Reduction Gear 


Fig. 1. Longitudinal sections of the turbine (a), reduction gear (b) and generator (c) with the distribution of rotor 8 relative vibration sensor 10 and 
identification of bearings in the steam turbogenerator: 1) turbine rotor bearing, 2) pinion bearing turbine side, 3) pinion bearing generator side, 
4) turbine thrust bearing, 5) gear wheel bearing turbine side, 6) gear wheel bearing generator side, 7) gear thrust bearing, 

8) generator bearing turbine side, 9) generator bearing free end side 


THE RESEARCH OBJECT 


The procedures concerning the operation and diagnosis 
of ship steam turbines have been examined based on 
a turbogenerator with a 1100 kW horizontal, four-stage 
condensation steam turbine AT42M produced by Mitsubishi 
Heavy Industry Ltd, which drives a synchronous generator 
via areduction gear. These turbogenerators are installed in 
the waste heat recovery systems of a VLCC tankers and work 
with the exhaust gas boiler co-operating in some operating 
states with the auxiliary boiler. Leaving the turbine, the 
superheated steam from the exhaust boiler flows to a vacuum 
condenser cooled with the seawater. The cross sections of the 
turbogenerator with its basic components, their bearings and 
the distribution of the measurement points monitored during 
its operation are shown in Fig. 1 [3]. 

As a rule, the turbogenerator is equipped with a monitoring 
and control system consisting of thermometers measuring 
local temperatures of the bearings, see Fig. 1: points 1-9, and 
PTR-M type thermocouples (thermo-resistance bulbs) for 
their remote monitoring. Thermal and flow parameters of the 
working media are controlled by thermometers, thermocouples 
and pressure gauges. 

Moreover, the list of elements of the control and 
measurement system includes: pressure switches with the 
autostart of the second generator at low live steam pressure, 
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thermostats controlling the oil temperature with the auto start 
of an independent oil pump in case of low oil pressure, and 
turbogenerator limit controllers (safety systems) with the 
emergency stop in case of: oil pressure drop below the lower 
limit, abnormal/excessive vibrations, high vapour pressure at 
the inlet to the condenser, and/or an excessive turbine speed 
(overspeed trip), doubled by the manual emergency “stop”. 

The turbogenerator is also equipped with a system to 
measure the turbine rotor shaft vibration in relation to its 
bearings (relative rotor/ bearing vibration), which is the VM- 
5G model, produced by Shinkawa, with an eddy-current sensor 
shown in Fig. 1 (Section 10). The model can be, at customer’s 
request, extended by an additional sensor. This system allows 
a remote visualization of displacements (vibrations) and is 
connected to the general marine power plant monitoring 
system, with the alarm set at 50 microns, as a standard, and 
the stop function activated when the limit of 80 microns is 
exceeded [3]. 


OPERATIONAL SERVICE AND 
MAINTENANCE OF TURBOGENERATORS 


The service and maintenance of marine steam turbogenerators 
during ship operation has the form of continuous (semi- 
automatic) monitoring process, with the maintenance resulting 
from working hours. This process consists in simultaneous 


Tab. 1. Sample working parameters of steam turbogenerator working at 750 kW load 


Measured Alarm/ limits 


State parameter values of working media 


Emergency stop 
Main 
1 stage 
Gland 
Exhaust [MPa] 
Main 
Exhaust 
Bearings 


Steam pressure [MPa] 


Steam temp. [°C] 


Lube oil pressure [MPa] 
Lube oil temp. [°C] 


System 
LO IN 
LO OUT 
Water IN 
Water OUT 
Turbine (1) 
Pinion gear (2) 


Medium temp. in cooler [°C] 


Bearing temp. [°C] 


Pinion gear (3) 
Wheel gear (5) 
Wheel gear (6) 


Gear bearing temp. [°C] 


Turbine (4) 


Thrust bearing temp. [°C] Gear (7) 
ear 


Turbine (8) 


Generator bearings temp. [°C] at side of Feesn) 


RMS rotor’s relative vibration [um] Turbine (10) 
control, recording and archiving of operating parameters and 
working process parameters. The list of monitored measures 
with sample values registered at turbogenerator load of 750 
kW and nominal speed of 1800 rev/min, after a total of 57,200 
hours from the last overhaul (including safety thresholds for 
alarms and emergency stops) is given in Tab. 1. 

As a rule, most parameters are recorded every four hours and 
archived automatically by the power plant monitoring system. 
This rule, however, does not include vibration measurements. 
In the unmanned engine room, some parameters are entered 
in the engine room log once a day. A comprehensive report 
with turbogenerator parameters is analyzed in the technical 
department of the ship-owner once a month. This report also 
includes turbogenerator working/maintenance hours, and 
results of the lube oil analysis. Selected elements of such 
a report are presented in Tab. 1. 

The only element of vibration diagnostics which is part of 
the standard equipment of the turbogenerator is the measuring 
system of the bearing related displacement — the relative 
turbine rotor/bearing vibration. It aims at a continuous control 
of RMS vibration level of the running turbogenerator, with 
a possibility to switch the alarm on when the permissible 
values are exceeded, and stop the turbogenerator in case of 
exceeding the absolute limits. However, it is not possible to 
register vibration level changes as a function of time, with 
further trend analysis. 


OPERATIONAL OBSERVATIONS AND 
MEASUREMENTS 


In the turbogenerator operation process, an increase of 
relative vibrations from 20 to 30 microns, i.e. by about 50% 
compared to those observed in the previous period have been 
recorded. This fact inspired checking the vibration levels in 


full range of turbogenerator - TG load. The registered vibration 
levels confronted with the limits are shown graphically in 
Fig. 2. 


2 
a 


N 


—¢— TG Vibration before O/H 
—#— TG Vibration after O/H 


n 
© 


a 


ta 


Turbine rotor relative vibration - RMS [um] 


850 


800 700 600 500 

TG Load [kW] 

Fig. 2. RMS value changes of turbine rotor s relative vibration signal in the 
turbogenerator working at 1800 rev/min, for two technical states: before 


and after overhaul o/h 


400 300 200 100 0 


The RMS values of the turbine rotor vibrations as a function 
of load changes at constant nominal rotational speed increases 
throughout the entire load range. 

The vibrations at smaller rotor speeds than the designed/ 
nominal value were also inspected, see Fig. 3 for a quasi-static 
realisation of the turbogenerator start-up process. The RMS 
values of the turbine rotor vibrations as a function of rotational 
speed when the turbogenerator worked in the off-load mode also 
increased throughout the full speed range, which is particularly 
noticeable in the critical speed range of 800-1100 rev/min. 

Comparing the results in the two presented cases with 
the data from earlier trials/delivery measurements and 
with the results with similar measurements done later after 


POLISH MARITIME RESEARCH, No 1/2013 33 


30 


Area of the first critical speed 


Nm 
ma 


Turbine rotor relative vibration - RMS [um] 


0 
100 200 300 400 500 600 700 800 


900 1000 


ion before O/H 


tion after O/H 


1100 1200 1300 1400 1500 1600 1700 1800 


Turbine rotor revolution [revs/min] 
Fig. 3. RMS value changes of turbine rotor 5 relative vibration signal in the turbogenerator idling at various revolutions — for two technical states: 
before and after overhaul o/h, with the indicated area of the first critical speed 


turbogenerator overhaul (see Fig. 2 and 3), a significant increase 
of the vibration level is observed, which suggests worsening 
of the technical/dynamic state of the turbogenerator. Despite 
the elimination of external disturbances caused by different 
states of ship load, hydro-meteorological conditions and/or the 
work of other engines and equipment, repeated measurements 
confirmed the increased unserviceability of the turbogenerator. 
Unfortunately, it was not possible to determine its location 
and source. Taking into consideration the location of the 
measuring transducer and the nature of the wear observed in 
the turbogenerator bearings it was assumed that a very likely 
reason of the noticeable increase of the vibration level was the 
increased clearance in the turbine bearing (Fig. la, point. 1). 

The operating experience gained so far [8] suggests that 
the currently used operational methods and procedures are 
inadequate and insufficient. Their effectiveness and credibility 
depend on a proper assessment done by a human operator and 
his/her ability to use the archived information. 

The lack of the fault development trend analysis limits the 
diagnosis of the turbogenerator operation to a two-state form: 
good condition - the state of fitness/technical ability, or poor 
condition - the state of unfitness/technical inability. Such an 
approach makes it impossible to assess the time of nearing 
the state of unfitness even if there is only a partial state of 
fitness. 

Therefore the scope of the diagnosis performed for 
maintenance purposes needs to be extended by a trend-based 


diagnostic inference making use of an additional apparatus 
for measuring characteristic parameters of turbogenerator 
vibration [4, 6]. 


ASSESSING THE TECHNICAL STATE OF 
THE TURBOGENERATOR WITH THE AID 
OF VIBRATION DIAGNOSTICS 


Since the manufacturing and installation technology 
was believed to be correct with respect to the examined 
turbogenerator, further research was oriented on searching for 
operating reasons of the observed changes of turbogenerator’s 
dynamic state which manifested themselves by elevated 
bearing vibration amplitudes. The intensity and nature of such 
vibrations depend on, among other factors, unbalanced mass 
inertia, quality of the alignment, thermal deformations, and/or 
bearing operation conditions [6]. 

In this situation it was decided to perform additional 
measurements of turbogenerator dynamics. It was done by 
an external service team with the aid of the Schenk apparatus 
and the measuring instrument Yibroport with the module 
for RMS velocity measurement [7]. The absolute vibrations 
were measured using a piezoelectric velocity transducer with 
a magnetic base mounted on a flat surface of the machine 
bearing housing, see the diagram in Fig. 4. The measurements 
were done in the vertical (V), horizontal (H) and axial (A) 
directions. 


Absolute vibration velocity 
H measurement at points | — 5 [mm/sec] 
and measurement directions: 


V — vertical 
H — horizontal 
A — axial 


Relative vibration displacement 
measurement at point 6 [um] 


D Generator D 


4 5 


Fig. 4. Schematic distribution of vibration signal measurement points in steam turbogenerator: 1) and 6) turbine rotor bearing, 2) reduction gear pinion 
bearings, 3) reduction gear wheel bearings, 4) generator bearing turbine side, 5) generator bearing free end side 
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Tab. 2. RMS values of absolute bearing house vibration velocities, measuring points 1-5, and relative turbine rotor/bearing vibration, point 6, for various 
turbogenerator loads and nominal speed 1800 rev/min 


The measured RMS values of the absolute bearing house 
vibration velocities, the measuring points 1-5, and of the 
relative turbine rotor/bearing vibration, point 6, are collected 
in Tab. 2 for different turbogenerator loads and the rated speed 
of 1800 rev/min [7]. 

Tab. 3 collects the reference criteria for assessing the 
dynamic state a turbogenerator taking into account the 
acceptable and limited RMS velocity values [4, 5, 6]. These 
criteria make the basis for creating a set of quantities being 
comparison classifiers. The obtained results were compared with 
the measured reference values from Tab. 3. The all specified 
measuring vibration levels were within the acceptable range, 
below 3 mm/s, and their values increased with the increasing 
turbogenerator load. However, the observed vibration levels 
were slightly higher than those recorded by other indications 
for the turbine bearings, especially in the axial direction of 
2.74 mm/s. This trend was also reflected in the results of the 
measurements of the relative turbine rotor/bearing housing 
vibrations — 27 microns, recorded using a standard eddy-current 
displacement transducer installed by the producer. The results 
of the measurements have proved that the turbogenerator was 
in state of partial fitness with the reduced technical performance 
of the turbine bearing of concern. It was the reason why the 
verification of its technical state with the suggested need to 
replace the bearings was included to the scope of the incoming 


overhauling after utilization of the resource hours (TBO — Time 
Base Overhauling). 

In order to enhance the results of the diagnostic inference, 
the measured vibration signals were analyzed using spectral 
analysis and the fast Fourier transform FFT [9]. The results 
of the analysis of the 1s and 24 order harmonics for various 
turbogenerator loads are shown in Tab. 4 [7]. 

The values and amplitudes of the 1s and 2"4 order harmonics 
of particular frequencies measured at points 1-5 reveal the 
similar nature as in the time domain, with the elevated vibrations 
indicated by higher 1st order harmonics for the turbine bearing, 
at the 2»4 order harmonics close to zero. This may indicate an 
increased bearing clearance - point 1, Tab. 4 and Fig. 5 a-c. On 
the other hand, at point 2 Tab. 4 the harmonic components of 
the 1s and 2"¢ order do not reveal increased values and do not 
indicate changes in the dynamic state of the turbogenerator at 
this measuring point. 

The analysis of the harmonic amplitudes of vibration signals 
measured at points 3 and 4 Tab. 4 shows the increased 2"4 order 
harmonics which may indicate a slight (permissible) unbalance 
of the turbine rotor. However, by comparing the frequency of 
the vibration amplitude waveform at the measuring point 6, 
Tab. 4 and Fig. 6 we can notice that there are increased values 
for both the 1stand 2™4 order harmonics which may be associated 
with a slight (permissible) alignment deviation. 


Tab. 3. Ranges of absolute vibration velocity RMS values which classify the dynamic state of turbogenerator [mm] 


Small vibration. Low wear of bearing. Low noise level. State of operational suitability/technical efficiency. 
Good condition, no disability. 


Noticeable vibration, often caused by increased wear of bearings. Increased noise level. Needs to find a reason, 
plan inspection/repair in the near future. The need for additional measurements and observations. Status of partial/ 
incomplete technical efficiency. Sufficient state, small/partial damage. Plan the repair. 


High level of noise and vibration. Too high temperature of bearings. Loss of oil film. Wear. Need to stop and 
7.1 — 18 | perform inspection/repair. State of unfitness/technical malfunction. The state of need for improvement, dangerous 
damage. Do the repair. 


Very high level of noise and vibration. Damage to the machine, emergency stop. The need for repair. State 
of unfitness/technical malfunction. Unacceptable state, catastrophic damage. Immediately stop the engine. 
Turbogenerator should be stopped automatically by safety systems. 


Tab. 4. 1" and 2” order harmonics of the absolute vibration velocity signals recorded at measuring points 1-5 
Jor different turbogenerator loads and revolutions 1800 rev/min 
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Fig. 5. Harmonic analysis of vibration signals measured in turbine rotor bearing at 880 KW: 
a) vertical direction (V), b) horizontal direction (H), c) axial direction (A) 


The spectral analysis of the waveforms and the 1st and 24 
order harmonic components confirm the conclusion formulated 
after time domain analysis about the technical state of the 
turbogenerator, in particular the state of the turbine rotor 
bearing, and the need for its overhaul. It also made it possible 
to infer about the dynamic state of other bearing nodes — other 
measuring points. The formulated conclusions and operational 
decisions were verified during turbogenerator overhaul in 
a shipyard. 


TURBOGENERATOR REPAIR INSPECTION 


During the verification, traces of normal/operational wear 
were observed on both active and passive reduction gear 
wheels and in their bearings, and also in the generator bearing 
on the turbine side and in the turbine and gear thrust bearings. 
The results of relevant clearance measurements done in these 


bearings are given in Tab. 5, entries 2-8, the presented values are 
within the design ranges [3]. However, significant traces of wear 
were observed on the lower shell in the turbine rotor bearing, 
see Fig. 7a, which was confirmed by visual/organoleptic 
inspection and measurement of bearing clearance, see Tab. 5, 
position 1. These results slightly exceeded the acceptable 
values. In the generator bearing, traces of wear were observed 
on the edges of the outer working surface of the upper shell on 
the free end of the shaft, see Fig. 7b. This may indicate a slight 
misalignment between the generator and reduction gear shafts 
on the side of the passive wheel, as a possible consequence 
of relatively low precision of assembly. The measured values 
of the bearing clearance were equal to the design dimensions 
given in Tab. 5, position 9. Fig. 7 shows the observed traces of 
operational wear on the working surfaces of the turbogenerator 
bearing shells which were replaced with new ones during the 
dry-dock overhaul. 


Tab. 5. Design clearances and clearance measurements in AT42M type turbogenerator journal bearings 


Design clearances [mm] 


Measured clearances [mm] 


Bearing name 


Turbine rotor bearing (1) 


Turbine side 
0.100 — 0.166 


Generator side Turbine side Generator side 


Pinion bearings (2 and 3) Reduction gear 


0.328 — 0.365 


Wheel bearings (5 and 6) Reduction gear 


0.125 — 0.175 


Generator bearings (8 and 9) 


0.175 


Thrust bearings (4 and 7) 
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Fig. 6. Harmonic analysis of vibration signals measured in generator free end bearing at 880 KW: 
a) vertical direction (V), b) horizontal direction (H), c) axial direction (A) 


j J 


Wear traces 


Wear traces 


Fig. 7. Operational wear traces in turbogenerator type AT42M bearings: a) turbine rotor bearing Isteam exhaust side; 
b) bottom half of generator free end bearing 5 


Taking into account the values of the recorded clearances 
(Tab. 5 position 1 and 9), the state of the bearings 1 and 5 (Fig. 4), 
was classified as partially fit and conditionally accepted for 
further use/operation. However, due to the specific nature of the 
whole turbogenerator as the “critical equipment” for readiness 
task in a safe and trouble-free operation of the ship, and also the 
bearing themselves as critical spare parts on board, a decision 
was made to replace these bearings with new ones. 


Within the framework of comprehensive evaluation of the 
technical state of turbogenerator components, the measurements 
of turbine shaft runout in specific control planes were done. 
The measurements showed slight deviation of the measured 
values from zero, of an order of + / - 0.02 mm, thus proving 
the absence of shaft deflection due to thermal deformation. The 
measurements of the turbine rotor unbalance which were done 
at the dynamic balancer indicated the unbalance values not 
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exceeding 10 grams: 8.0 grams in the turbine bearing support 
area, and 6.4 grams for the reduction gear of the active wheel 
bearing. These values were reduced to less than | gram, 1.e. 
990 mg and 239, respectively [7]. 
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CONCLUSIONS 


The confrontation of the results of technical state verification 
of the examined turbogenerator with the analysis of the 
vibration measurements shows that they are compatible. 
This justifies the opinion about correctness of the diagnostic 
inference based on the measured vibration signals in both 
the time and frequency domains. Moreover, the vibration 
diagnostics turned out to be a useful and reliable source of 
information on the dynamic state of the turbogenerator and 
its components. 

This conclusion was confirmed by further measurements 
performed after completing the overhaul and assembling the 
turbogenerator and its components, including the coaxial 
connection of gear and generator shafts - checking of 
coupling alignment, which showed a significant decrease 
in the relative turbine shaft/bearing vibrations (the RMS 
change curves after overhaul in Fig. 2 and 3). Absolute 
vibration measurements after repairs were not performed 
—this decision was made based on the results of the overhaul 
and subsequent turbogenerator operational performance 
reports. 

The presented method of turbogenerator service and 
maintenance indicates that it is advisable to include the 
trend analysis of changes of its technical components 
as an important tool in decision-making processes. This 
would require supplementing the method by predictive 
maintenance models of turbogenerator operation. 
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Multi-stage magnetic-fluid seals for operating 
in water — life test procedure, test stand and 
research results 
Part Il 
Results of life tests of multi-stage magnetic - fluid seal 


operating in water 


Leszek Matuszewski, Ph.D. 
Gdansk University of Technology, Poland 


ABSTRACT 


The edge criteria for magnetic-fluid seals durability research were described in previous 
article [24] and now the tests results are presented. The tests were made with three magnetic 
fluids using various magnets quantity in magnetic assemblies. If values obtained in the 
repeated test deviated by more than 10 % from results of the first test an additional control 
test was made for final characteristics. The characteristics are: critical pressure, critical 
motion velocity and working life of the seal. The test stands were prepared for particular 
seals dimensions and results of tests are presented in the form of tables, cross section draws 


and diagrams. In conclusion one can see that multi-stage magnetic-fluid seals could be efficiently used 
in water for rotating shaft seals in a limited range of motion velocity and cycles quantity. Nevertheless 
a particular applications can be provided by magnetic-fluid seals only. 


Keywords: shaft sealing; magnetic fluid; utility water; seal durability; critical pressure; critical motion velocity 


INTRODUCTION 


This paper presents the most significant results of the 
research projects performed in the years 2009+2012. All the 
projects were carried out on the test stands located at the 
Faculty of Mechanical Engineering and Robotics, Mining and 
Metallurgy Academy (AGH), Cracow, Poland. 

As described in the first part of the paper [24] the research 
methodology of the performed projects ensured determination 
of static penetration pressure and limiting working speed of the 
seal for each of its tested versions. Each of the tests introduced 
important information which made it possible to draw seal 
life curve and determine the largest motion speed at which 
failure-free operation (lack of leakage) is possible up to the 
basic number of cycles (equal to 1 million rotations), assumed 
for the seal. 

This author conducted the research projects on the 
application of magnetic fluid (MF) seals operating in utility 
water in collaboration with the Laboratory of Seals and 
Application of Magnetic Fluids, AGH, [12, 13, 14]. The tests 
performed on the research stand adjusted to operation of seal 
in liquids [15] showed that application of some commercial 
magnetic fluids to MF seals being in direct contact with water in 
operational conditions typical for ocean engineering, mainly to 
driving systems, is possible. In the tests which had an utilitarian 
character, multi-stage MF seals representing models of real MF 


seals were used and their operational conditions were limited 
to the selected technical task. In this paper, due to a limitation 
of its volume, only some example diagrams are presented. 


1. RUN OF THE TESTS 


In the tests in question the MASTL1 — V2 test stand was 
used. The stand contains a special test head intended for 
realizing the assumed research aims, driving system composed 
ofan electric motor and its controller, as well as a measurement 
system composed of gauges and transducers for measuring 
pressure, motor’s torque and operational temperature. The 
stand is fitted with multi-channel measuring instruments for 
measurement data acquisition, processing and recording. 
Their detailed description was given in the first part of the 
paper [24]. 


1.1. Tests of seal critical pressure 


The tests of seal static critical pressure were performed in 
compliance with the procedure described in the paper titled 
“Multi-stage magnetic-fluid seals for operating in water — life 
test procedure, test stand and research results, Part i; [24]. 
The pressure in the test chamber was gradually increasing up 
to occurrence of penetration of MF seal. The sealed medium 
was pure utility water. Value of seal static critical pressure is 
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the basis for determining pressure values to be set to obtain 
the appropriate values of the relative pressure p,/p,,,; = 0.8, 
1.6, 2.4, complying with the test programme. The tests were 
conducted for three magnetic fluids: FLS 040.040, FLA 002.25, 
FLA 003.45, with using 14 magnets in the magnetic system, and 
additionally with 10 and 18 magnets for the fluid FLA003.45. 
In all the tests the nominal gap height for magnetic fluid was 
equal to 0,1 mm at 100u1 magnetic fluid dosage applied to each 
of the sealing stages. Each test was repeated. If values obtained 
in the repeated test deviated by more than 10 % from results of 
the first test an additional control test was made. 
In Fig. 1.1 the example diagram of pressure run during 

determination of statical critical pressure is presented. 

Measurement No.21.1 — FLA 003.45, Gap 0,1 

Statical penetration pressure 
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Fig. 1.1. The example diagram of determination of statical penetration 

pressure for the fluid FLA 003.45 at 0.1 mm gap height and three-stage 

sleeve: the pressure curve peak value represents the value of the statical 
penetration pressure (for the whole three-stage seal) 


1.2. Tests of critical motion speed of MF seal 


The tests of MF seal critical motion speed were performed 
to determine the lowest value of motion speed at which loss of 
sealing capability ofa seal occurs “immediately”. Information 
on critical motion speed value makes performing seal life test 
easier. The tests were conducted with gradually increasing 
speed in compliance with the procedure described in the first 
part of the paper, for the magnetic fluids: FLS 040.040, FLA 
002.25, FLA 003.45, and with the use of the seals of 0.1 mm 
nominal gap height and 100 u1 dosage of magnetic fluid applied 
to each sealing stage. According to the test programme relative 
pressure values were used. Each measurement was repeated to 
eliminate random results. 

For FLA003.45 fluid additional measurements were made 
at the relative pressure values p,, = p,/Py1 = 0.8, 1.6, 2.4 and 10 
and 18 permanent magnets used in the magnetic system. 

In Fig. 1.2 is presented the example diagram of run of 
determination of limiting penetration speed recorded during 
the tests. 


1.3. Life tests of MF seal 


The MF seal life tests were carried out to determine the 
seal motion speed at which MF seal life for the assumed seal 
operating conditions and gap geometry would be equal to 1mln 
cycles (number of rotations). To determine one-million-cycle 
life it was necessary to perform series of measurements at 
constant, but lower and lower for each successive measurement, 
seal motion speed. The series of measurements were started 
from a motion speed equal or close to a seal critical motion 
speed specific to a seal under the test. The tests were conducted 


40 POLISH MARITIME RESEARCH, No 1/2013 


Measurement No.18.1 — FLs 040.40, Gap 0,1X3 
Limiting speed test 
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Fig. 1.2. The example diagram of determination of the seal critical speed 
for the fluid FLS 040.40 at 0.1 mm gap height and three-stage sleeve. 
The instance of exceeding the critical speed is signalled 
by sudden drop of pressure 
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Fig. 1.3. The example diagram of pressure and temperature runs in function 
of duration time counted from the instant of disclosure of a leakage. The MF 
seal with FLA 003.45 magnetic fluid at 10 rps; the relative pressure of 2.4, 
the gap height of 0.1 mm with 3 edges; a) temperature change of the seal 
during its leakage, b) change of drag- to-motion moment of the seal during 
its leakage. 


for three magnetic fluids: FLA 002.25, FLS 040.040, FLA 
003.45, and with the use of the seals of 0.1 mm nominal gap 
height and 100 ul dosage of magnetic fluid applied to each 


sealing stage. The relative pressure values p, = P,/Piyt = 
= 0.8, 1.6, 2.4 were used. For FLA003.45 fluid additional 
measurements were performed with 10 and 18 permanent 
magnets used in the magnetic system in accordance with 
a respective programme. 

The life criterion will be satisfied if during the test at an 
assumed shaft rotational speed no leakage is observed after 
one million of cycles, then the test is stopped and the speed is 
considered to be corresponding to one-million-cycle life. In 
Fig. 1.3 the example diagrams of run of MF seal life test, is 
presented. In Fig. 1.3a is shown an increase of temperature at 
the instant of leakage, caused probably by penetration of water 
ofa higher temperature into cooled zone of the seal. In Fig. 1.3b 
is Shown a drop of moment of drag to motion of the seal at the 
instant of leakage, caused by loss of the magnetic fluid from 
the seal. During this test the leakage was very intensive. 


2. RESULTS OF THE TESTS 


The test results are collected in the form of tables, curves 
and diagrams. In the tables are given values of test parameters 
and measurement results obtained during the tests, and on 
the curves and diagrams are shown values of the achieved 
results in function of the parameters of the tested seals or of 
test running time. 

The test results for three kinds of the tested magnetic fluids 
are elaborated seperately for each of the standard research 
tests: 

e total statical critical pressure of a seal, 
e critical motion speed of a seal, 
e long-term operation speed of a seal. 


On the curves and diagrams are presented the obtained 
values in function of particular testing factors, as well as in the 
form of comparisons of various parameters of the tests. 

All the tests were conducted with the use of the same sleeve 
with three sealing stages, at 0.1 mm gap and 100 ul dosage of 
magnetic fluid applied to each of the gap. A variable design 
parameter was number of permanent magnets in the magnetic 
system: 14 magnets were used for three kinds of the tested 
magnetic fluids, and additionally 10 and 18 magnets in the 
tests with FLA 003.45 magnetic fluid. 


2.1. Test results of total statical critical pressure 
of the seal 


In Fig. 2.1 are shown the penetration pressure values for the 
seal with three sealing stages, obtained from the tests with three 
tested magnetic fluids. The penetration pressure measurements 
were made after one hour of fluid stabilization in the seal. 


Total statical critical pressure of the seal, 
edges quantity z = 3, magnets quantity i= 14, gap 6 = 0.1 mm 
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Fig. 2.1. Comparison of values of the total statical critical pressure for 
three tested magnetic fluids 


In the successive series of the tests were determined 
penetration pressure values for the seal of three sealing stages, 


obtained for FLA 003.45 magnetic fluid from the testing with 
the use of different number of permanent magnets in the sealing 
system. In all the tests 100 ul dosage of magnetic fluid was 
applied to each sealing stage. 


2.2. Test results of critical speed of the seal 


In Tab. 1 are collected the critical speed test results at 
different values of relative pressure for three magnetic fluids 
and 14 magnets used in the magnetic system of the seal. 


Tab. 1. Test results of critical speed for MF seal. Number of permanent 
magnets in the magnetic system: 14 


FLA 00.2.25 | FLA 003.45 | FLS 040.040 


Relative 
pressure P,/Pyrr 


Critical speed [rps] 


In Fig. 2.2 are shown the results of limiting speed tests for 
three magnetic fluids at the relative pressure p,, = p,/P,,1 = 0.8, 
with 14 permanent magnets. 


Critical speed values for three magnetic fluids. Relative pressure 
P,/Pus = 0.8, edges quantity z = 3, magnets quantity i = 14, 
gap ô= 0.1 mm, liquid volume v = 100 ul 
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Fig. 2.2. Comparison of limiting speed values for the tested magnetic fluids 
at the relative pressure py, = P,/Px = 0.8, and with 14 permanent magnets 
used in the sealing system 


In the case of two magnetic fluids the critical speed exceeded 
the value of 100 rps at the relative pressure p, = P Pir = 
= 0.8. Because of mechanical limitations in operation of the 
test stand it was not possible to increase the rotational speed of 
the seal in question. In the successive test series were obtained 
measurement results of critical speed for three magnetic fluids 
at the relative pressure p, = p,/P,,1 = 1.6 and 14 magnets used 
in the sealing system, and the critical speed values of the tested 
fluids were compared to each other at the relative pressure p, = 
= P/P = 2.4 and 14 magnets. 

When comparing the critical speeed values obtained for 
the magnetic fluids at the relative pressure pẹ, = p,/P,.1 = 2.4 
attention should be paid to significant differences between 
real pressure values set in the test chamber. During the tests 
with FLA 003.45 and FLS 040.040 fluids the pressure in the 
test chamber was over twice higher than that assumed in the 
tests with FLA 002.25 fluid. It indicates significant influence 
of absolute pressure value on operational efficiency of MF 
seals working in water. Influence of magnetic field on limiting 
speed for FLA 03.45 fluid was tested additionally at the relative 
pressure Py = p,/p,,,, = 1.6 and different number of magnets in 
the magnetic system of the seal. It was demonstrated that along 
with increasing number of magnets installed in the test head 
(i.e. increasing magnetic flield intensity) critical speed value 
for FLA 003.45 fluid also increases. 
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Tab. 2. Results of seal life tests in function of seal motion speed for FLA 002.25 magnetic fluid 


Kind of 


Stabilization time: 1h magnetic fluid 


FLA 002.25 


Nominal gap height: 0.1 mm 


Statical critical pressure 
of the seal: total value 
/value per one stage. 


Pic/Pior [MPa] 


0.0890245/0.0296748 14 


Number of magnets in the sealing 
system 


Critical speed 
[rps] 


Relative pressure (set pressure) 
Pw = P/Pxn (P,) [MPa] 


Rotational speed [rps] 


Number of rotations to leakage 
(time of operation [h]) 


100 1700 (0.00472) 
1840 (0.00639) 


0.8 


(0.023) over 100 


99936 (0.4627) 


2251800 (13.9) 
No leakage. exceeded 1mln rotations 


80 (0.05383) 


2450 (0.00972) 


171780 (0.2071139) 


257040 (2.38) 


388800 (5.4) 


648000 (18) 


100000 (55.6) 
No leakage. exceeded I mln rotations 


356 (0.02583) 


3000 (0.01389) 


5400 (0.05) 


34200 (1.9) 


Notations which appear on the below presented drawings, are explained on the example of Fig. 2.3: 


FLA 002.25 - symbol of the tested magnetic fluid 

Z=3 - number of stages of the test sleeve 

6=0,1 - height of the test gap, [mm] 

P,=0.023 - value of the set operational pressure, [Mpa] 


P/P = 0.8 - the relative pressure equal to 0.8, the ratio of the set operational pressure p, and the statical critical pressure per 
one stage, p,,,, (in the case of three-stage sleeve: 1/3 of the total penetration pressure) 
353532.6 - full distance passed by the seal up to leakage occurrence, [m]. 


2.3. Life test results of the seal with the basic 
number of permanent magnets in the magnetic 
system (equal to 14 pieces) 


In Tab. 2 are collected the life test results of the seal with 
FLA 002.25 magnetic fluid and 14 magnets installed in the 
head. Number of rotations to leakage constitutes a seal life 
measure. 

In Fig. 2.3 is presented the diagram of number of rotations 
to leakage in function of motion speed of the seal with FLA 
002.25 magnetic fluid at the relative pressure pọ, = p/P, = 
0.8. On the diagram is shown the trend line in the form of the 
function: 


y = -6.47|n(x) + 137.7 


In Fig. 2.4 is presented the diagram of number of rotations 
to leakage in function of motion speed of the seal with FLA 
002.25 magnetic fluid at the relative pressure p, = p,/Py4 = 
1.6. On the diagram is shown the trend line in the form of the 
function: 


= -7.88In(x) + 123.7 


In Fig. 2.5 is shown the diagram of number of rotations to 
leakage in function of motion speed of the seal with FLA 002.25 
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magnetic fluid at the relative pressure pẹ = p,/P,,,1 = 2.4. On the 
diagram the trend line is given in the form of the function: 


y = -22.94In(x) + 241.39 


The successive series of tests were conducted for 
different kinds of magnetic fluids, i.e. the FLA 003.45 and 
FLS 040.040, as well as with different numbers of magnets 
located circumferentially: 10, 14 and 18. Relevant tables and 
diagrams are not included in the paper due to limitation of 
its volume. 


Number of rotations to leakage in function of motion speed for 
Pr/ Prt = 0.8, z=3 and 5=0.1 mm, p = 0.023 MPa, magnetic fluid FLA002.25 
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Fig.2.3. Number of rotations to leakage in function of motion speed of the 
seal with FLA 002.25 magnetic fluid at the set operational parameters 
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Fig. 2.4. Number of rotations to leakage in function of motion speed of the 
seal with FLA 002.25 magnetic fluid at the set operational parameters 


Number of rotations to leakage in function of motion speed for 
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Fig. 2.5. Number of rotations to leakage in function of motion speed of the 
seal with FLA 002.25 magnetic fluid at the set operational parameters 


3. COMPARISON OF TEST RESULTS 


Series of the tests aimed at determination of MF seal life 
in function of magnetic field intensity, external pressure and 
kind of magnetic fluid, were performed. 


3.1. Seal life in function of magnetic field 
intensity in the system 


In Fig. 3.1 are presented the diagrams of number of rotations 
to leakage in function of rotational speed of the seal for different 
values of magnetic field intensity in the gap filled with magnetic 
fluid, dependent on number of magnets applied in the magnetic 
system of the seal. Results of the tests for FLA 003.45 magnetic 
fluid at the relative pressure p, = p,/p,,, = 1.6 and with 10, 14 
and 18 permanent magnets used in the sealing system, were 
compared to each other. 


Number of rotations to leakage in function of motion speed at 
different number of permanent magnets in the sealing system for 
Pr/Pxrt = 1.6, z=3 and 5=0.1 mm, magnetic fluid FLA003.045 
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Fig. 3.1. Influence of change in magnetic field intensity on number 
of rotations to leakage in the seal 


3.2. Seal life in function of set value of relative 
pressure 


In Fig. 3.2 are presented the diagrams of number of rotations 
to leakage in function of rotational speed for different values 
of relative pressure set in the tests with FLA 002.25 magnetic 
fluid. The tests were carried out at the relative pressure values 
Py = P/P = 0.8, 1.6, 2.4 and the same number of magnets 
equal to 14. 


In Fig. 3.3 are presented the diagrams of number of rotations 
to leakage in function of rotational speed for different values 
of relative pressure set in the tests with FLA 003.45 magnetic 
fluid. The tests were carried out at the relative pressure values 
Pw = P/P = 9.8, 1.6, 2.4 and the same number of magnets 
equal to 14. 


Number of rotations to leakage in function of rotational 
speed at different relative pressure values p;/Pkri = 0.8, 
1.6, 2.4 z=3 and 5=0.1 mm, magnetic fluid FLA002.25 
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Fig. 3.2. Number of rotations to leakage in function of rotational speed at 
different relative pressure values for FLA 002.25 magnetic fluid 


Number of rotations to leakage in function of motion speed at different relative 


pressure values in the sealing system for p,/py, = 0.8, 1.6, 2.4 z=3 and 
5=0.1 mm, magnetic fluid FLA003.045 
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Fig. 3.3. Number of rotations to leakage in function of rotational speed at 
different relative pressure values for FLA 003.45magnetic fluid 


In addition, were investigated and determined numbers of 
rotations to leakage in the seal in function of rotational speed 
for different values of relative pressure set in the tests with FLA 
040.040 magnetic fluid, i.e. at the relative pressure values p„ = 
= P/P = 0.8, 1.6, 2.4, and the same number of magnets equal 
to 14. The tests were repeated for different values of relative 
pressure set in the tests with another magnetic fluid, i.e. FLA 
003.45, and 10 magnets applied in the magnetic system. The 
tests were carried out at the relative pressure values pẹ = P/P = 
= 0.8, 1.6, 2.4. 


3.3. Seal life in function of kind of magnetic fluid 


In Fig. 3.4 are presented the diagrams of number of rotations 
to leakage in function of rotational speed for three magnetic 
fluids: FLA 002.25, FLA 003.45 i FLS 040.040, at the set 
relative pressure value p, = p,/p,,; = 0.8. The tests were carried 
out in the magnetic system with 14 magnets. 

Moreover, were investigated and determined numbers of 
rotations to leakage in the seal in function of rotational speed 
for three magnetic fluids: FLA 002.25, FLA 003.45 i FLS at the 
relative pressure value p, = p,/p,,; = 1.6 and with 14 magnets 
used in the magnetic system. 

The tests were repeated and numbers of rotations to leakage 
in function of rotational speed were determined also for three 
magnetic fluids: FLA 002.25, FLA 003.45 1 FLS 040.040 but at 
the set relative pressure value p,, = p,/p,,,; = 2.4. The tests were 
carried out in the magnetic system of 14 magnets. 
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Number of rotations to leakage in function of rotational speed 
for different magnetic fluids and the relative pressure value p,/Pxri = 0.8. 
Number of magnets 14, z=3 and 6=0.1 mm. 
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Fig. 3.4. Number of rotations to leakage in function of rotational speed for 
different magnetic fluids and relative pressure value equal to 0.8 


3.4. Rotational speed for seal life of one million 
number of rotations 


In Tab. 3 are presented the rotational speed values at which, 
during testing the seals, one million number of rotations was 
exceeded. The table contains results obtained for three magnetic 
fluids at different numbers of magnets, i.e. 10, 14, 18, in the 
magnetic system and different relative pressure values pẹ, = 
P/P = 0.8, 1.6, 2.4. 


Tab. 3. The rotational speed values at which one million number of rotations 
was exceeded, in function of number of magnets in the magnetic system, 
for three magnetic fluids and different relative pressure values 
Py “PPr = 0.8, 1.6, 2.4 


Rotational speed values 
[rps] 

Relative pressure p,, = p,/Pyr 

0.8 1.6 2.4 
FLA 002.25 below 5 
FLA 003.45 below 1 
FLS 040.040 2 
FLA 003.45 10 
FLA 003.45 


Number| Kind of 


of magnetic 
magnets 


fluid 


not tested 


not tested | not tested 

In Fig. 3.5 is presented a comparison of rotational speed 
values for three tested fluids, at which one million number of 
rotations was reached at the relative pressure pẹ, = P,/Piri = 
= 0.8. 
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Fig. 3.5. Comparison of rotational speed values at which one million 
number of rotations was reached for three tested fluids at the relative 


pressure Py = P,/Pyrt = 0.8 
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The tests were repeated for the tested fluids at the relative 
pressure values p, = p,/p,,,, = 1-6 and 2.4, respectively. 

Fig. 3.6 is shown a comparison of rotational speed values 
at which the condition of reaching 1 mln number of rotations 
was fulfilled for particular fluid FLA002.25 at three relative 
pressure values p,, = p,/P,,1 = 0.8, 1.6, 2.4. 


Rotational speed values corresponding to seal life of 1 million number of 
rotations, fluid FLA002.25, relative pressure Pr/Prri = 0.8, 1.6, 2.4 
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Fig. 3.6. Comparison of rotational speed values at which one million 
number of rotations was reached for FLA 002.25 fluid at the relative 
pressure values Py, =P/Pur = 0.8, 1.6, 2.4 


The tests were repeated for the magnetic fluids: FLA 
003.45 and FLS 040.040, at the relative pressure values p, 
= p/Pu = 0.8, 1.6, 2.4 and in configurations with 10, 14 and 
18 magnets. 

In Fig. 3.7 is presented an example comparison of rotational 
speed values at which the condition of reaching 1 mln number 
of rotations was fulfilled for FLA 002.25 magnetic fluid at 
three relative pressure values p, = p,/P,,; = 0.8, 1.6, 2.4 and 
with decreased number of magnets in the sealing system, 1.e. 
10 pieces. 


Rotational speed values corresponding to seal life of 1 million number of 
rotations, fluid FLA003.45, relative pressure p,/pyy1 = 0.8, 1.6, 2.4, 
magnets quantity z=10 
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Fig. 3.7. Comparison of rotational speed values at which one million 

number of rotations was reached for FLA 002.25 fluid at the relative 

pressure values Py, =P/Pur = 0.8, 1.6, 2.4 and with decreased number 
of magnets in the sealing system, i.e. 10 pieces 


4. DISCUSSION OF TEST RESULTS 


In the frame of the conducted experimental projects on the 
MF seal with three sealing stages operating in water, three kinds 
of magnetic fluid were tested in conditions of different values 
of factors which influence operational effectiveness of the seal. 
The three effectiveness measures were applied: 

e critical pressure of the seal (penetration pressure), 

e critical motion speed (,,immediate” penetration motion 
speed), 

* motion speed for a limited seal life (motion speed for one- 
million-cycles life). 


The two external factors and one design factor affecting 
seal operation were taken into account: 
e external factors: 

- water pressure in the test chamber, 

- seal motion speed, 
e design factor: 

- number of magnets in magnetic system of the seal. 


To estimate influence of the factors on the seal, 120 test 
runs lasting from about 1 hour up to over 30 hours, were 
performed. 

The test results were elaborated separately for each kind 
of the performed tests (critical pressure, critical speed, long- 
term operation speed), and then they were compared mutually 
to make it possible to assess operational features of the seals 
working with different magnetic fluids and in different external 
conditions. 


4.1. Critical pressure tests of the seal 


Results of the tests are presented in Fig. 2.1 and 2.2. In 
each case the total critical pressure p,,; of the seal with three 
sealing stages was tested by using 100u1 magnetic fluid dosage 
per one stage. 

In the tests of the seals with 14 magnets the highest critical 
pressure value equal to 0.217 MPa was achieved for FLA 
003.45 magnetic fluid of the saturation magnetization M, = 
45 kA/m [17]. A little lower value equal to 0,198 MPa was 
obtained for FLS 040.040 magnetic fluid for which the value 
of M, = (40 + 50) kA/m was specified by its producer [17]. It 
seems that for the tested fluid the value of M, ~ 40 kA/m is 
more realistic. For FLA 002.25 magnetic fluid of the saturation 
magnetization M, = 25 kA/m the critical pressure value equal 
to 0,089 MPa was obtained. 

The tests on FLA 003.45 magnetic fluid with the use of 10 
and 18 permanent magnets show that either at this smaller or 
greater number of magnets used in the magnetic system the 
critical pressure value was lower than that obtained for 14 
magnets. 


4.2. Critical speed tests of the seal 


The critical pressure is defined as the smallest value 
of rotational speed at which — at given constant values of 
parameters — a leakage occurs during continuous increasing the 
rotational speed from 0 to its maximum value of 100 rps. As the 
maximum value of rotational speed was limited by torque meter 
characteristics in some cases it was not possible to determine 
this value because it exceeded the 100 rps - limit. The critical 
speed was tested at each of three relative pressure values, for 
each of the magnetic fluids. For FLA 003.45 magnetic fluid 
the tests were also performed with the use of different number 
of permanent magnets (in configurations of 10, 14, and 18 
magnets). 

The test results are collected in Tab. 1 and presented in 
Fig. 2.3 through 2.6. In the tests, at the relative pressure values 
Pw = PY Pir = 0,8 and p, = p,/p,.1 = 1.6 the greatest limiting 
speed values were observed for FLA 002.25 and FLA 003.45 
magnetic fluids. In the tests at the relative pressure p, = P,/Pkri = 
= 0.8 the 100 rps speed value was exceeded for both the fluids, 
whereas at the relative pressure p,, = p,/P,,; = 1.6 the speed 
values amounted to 80 and 83 rps, respectively. However the 
magnetic fluids in question produced on the basis of the same 
primary liquid - SILOXAN, greatly differ to each other in 
saturation magnetization value: M,=25 kA/m and M,=45 kA/m, 
respectively. And, smaller limiting speed values were obtained 


in the tests with the use of FLS 040.040 magnetic fluid produced 
on the basis of silicon liquids. 

In the tests at the relative pressure p, = p,/p,,; = 2.4 the 
greatest limiting speed value equal to 76 rps was achieved in the 
tests with the use of FLA 002.25 fluid of the smallest saturation 
magnetization value. The limiting speed values for the 
remaining fluids are about six times smaller despite their almost 
twice greater saturation magnetization. It should be taken into 
consideration that absolute pressure values set in the tests with 
FLS 040.040 and FLA 003.45 fluids were significantly greater 
than that used in the tests with FLA 002.25 fluid. The value 
resulted from the assumed test concept according to which 
absolute pressure values are derived from relative pressure 
values: p,, = p,/p,,;- Since great differences appeared between 
statical penetration values for the magnetic fluids in question, 
pressure values set in the tests with FLA 002.25 magnetic fluid 
were consequently much lower than those set in the tests with 
the remaining fluids. 

In the limiting speed tests at the number of magnets: 10, 14, 
18, performed with the use of FLA 003.45 magnetic fluid at the 
relative pressure value p, = p,/P,,1 = 1.6 (Fig.2.6) the greatest 
limiting speed value of 91 rps was obtained for 18 magnets, and 
the smallest value of 75 rps for 10 magnets in number. 


4.3. Life tests of the seal 


Life of MF seals was the most thoroughly investigated 
issue in this research project. In most of the tests the number of 
cycles (number of rotations) performed by a seal up to leakage 
was taken as its life criterion. Full sets of the tests with 14 
permanent magnets in the magnetic system and at the relative 
pressure values p,, = p,/P,,1 = 0.8, 1.6, and 2.4 were performed 
for three kinds of magnetic fluids described in article [24]. 
The tests were aimed at checking relation between number of 
performed cycles and rotational speed set in test. 

For FLA 003.45 magnetic fluid were additionally performed 
tests with 10 and 18 magnets installed in the test head. One 
million rotations of the seal were assumed to be the basic 
number of cycles. For so defined basic number of cycles was 
determined the largest motion speed at which no leakage 
was observed in the tested seal. The quantity was called one- 
million-cycle life. It should be stressed that in some tests for 
the above mention speed a greater number of cycles without 
leakage was achieved, whereas in other tests it was decided 
not to continue running up to the basic number of | mln cycles 
as it would require to do long-lasting tests not planned in the 
project in question. 

The results of the seal life tests are collected in Tab. 2 
through 6 and shown on the diagrams in Fig. 2.7 through 3.6 

The seal life diagrams for three kinds of magnetic fluids 
tested at three relative pressure values pẹ = p,/Py = 0.8, 1.6, 
and 2.4 indicate large differences in values of the speed for 
which number of cycles to leakage, corresponding to them, was 
obtained. The differences appear both if to compare the tests 
for particular fluids at different relative pressure values and if 
they concern comparison of the results obtained from the tests 
of different kinds of fluids. 

The largest speed for long-term operation was reached in 
the test with FLA 002.25 fluid at the relative pressure value 
Pw = P/P = 9.8. It was equal to 45 rps, and in the test the basic 
number of 1 mln cycles was exceeded. The smallest speed for 
long-term operation was reached in the test with FLA 003.45 
fluid at the relative pressure value p, = p,/Pyy = 2-4. 

It was equal to | rps, and in the test the basic number of 
cycles was not reached. However it should be stressed that in 
the test with FLA 002.25 fluid in which the largest long-term 
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operation speed was obtained the pressure set in the test was 
equal to 0.0232 MPa, whereas in the test with FLA 003.45 
fluid, in which the smallest motion speed was obtained, the 
set pressure was equal to 0.1824 MPa. The experimental data 
obtained from all the tests were approximated by logarithmic 
trend lines of different coefficients depending on a given test. 

The seal tests under different magnetic field intensity 
(different numbers of permanent magnets in the magnetic 
system of the seal) were performed for FLA 003.45 magnetic 
fluid at the relative pressure value pẹ = p,/p,, = 1.6. Comparison 
of results of the tests indicates that the largest limiting speed 
of 91 rps was obtained in the test with the use of the magnetic 
system of the highest magnetic field intensity (18 permanent 
magnets). And, the values of operational speed for the basic 
number of 1 mln cycles without leakage were close to each 
other. From the comparison of the results the value of about 
15 rps was achieved. 


5. RECAPITULATION AND CONCLUSIONS 


According to the applied testing methodology it was 
decided, for each tested version of the seal, to determine the 
following: 

e statical penetration pressure of the seal, 
e limiting operational speed of the seal. 


Each of the measurements introduced important information 
which made it possible to reach the research aim as follows: 
e to determine seal life curve and the highest motion speed for 

which it is possible to achieve, without leakage, the basic 

number of cycles (rotations) equal to 1 mln as assumed for 
the tests in question. 


In the tests water pressure value set in the measurement 
chamber was a crucial factor. The factor was accounted for 
in the form of the relative pressure defined as the ratio of the 
set pressure and the critical pressure i.e. the seal penetration 
pressure per one sealing stage, determined in the statical test, 
(Pw = Py/Puct): 

The using such concept of setting the testing pressure of 
water made it possible to simply account for influence of the 
factor on seal life for particular magnetic fluids, but in the same 
time it made comparing life values of different fluids more 
difficult because of large differences in set testing pressures 
resulting from highly different values of the statical penetration 
pressure of the tested magnetic fluids. For instance the highest 
long-term operational speed equal to 45 rps was reached in 
the tests with FLA 002.25 magnetic fluid at the set pressure 
value of 0,0237 MPa (p,, = P/P = 9.8) resulting from its 
low penetration pressure value equal to 0,0297 MPa per one 
sealing stage. As in the analogous test at the value p,, = p/P, = 
= 0.8 for FLA 003.045 fluid (of the penetration pressure per 
one sealing stage equal to 0.072 MPa) it was required to set the 
water pressure of 0.058 MPa, the long-term operation speed 
of the seal with the fluid in question amounted to 20 rps only. 
The results show that the pressure greatly influences life of MF 
seal operating in water. 

The limiting operational speed of seal constitutes an 
important testing information as it makes it possible to estimate 
speed value from which determining seal life curve should be 
started - and consequently — determining the seal life for the 
basic number of cycles equal to 1 mln. However this parameter 
is not very important from practical point of view. 

The tests with various numbers of magnets in the sealing 
system revealed influence of magnetic field intensity on limiting 
speed of magnetic fluids, but they did not confirm a significant 
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influence of the factor on life of MF seals operating in water. 
Results of the life tests for majority of the tested magnetic 

fluids and seals show similar picture on the diagrams. At 
higher rotational speed values a short period of operation to 
leakage (small number of cycles) is observed. At lower motion 
speed values failure-free operation period is much longer and 
distinctly tends to the speed value which could be taken as that 
for long-term operation. Such character of seal life curve could 
be best described by a logarithmic trend line. 

The performed tests make it possible to offer the following 

conclusions: 

1. The tests confirmed correctness of the assumed testing 
methods. 

2. The test stand design is appropriate and allows to carry out 
the tests correctly. 

3. The tested magnetic fluids produced on the basis of liquids 
with strong hydrophobic properties (silicon liquid, siloxan) 
may be applied to rotational seals operating in water, 
however only in a limited range of motion speed. 

4. The tests confirmed that the assumed research hypotheses 
were right: 

e Motion speed influences life of MF seal operating in 
water: at lower speed values a significant increase of 
seal life is observed. 

e Water pressure influences MF seal life: at higher 
pressure values set in the tests a significant decrease of 
MF seal life is observed. 

5. For each of the tested fluids it was possible to determine 
a motion speed value at which the assumed life of 1 mln 
cycles (number of rotations) was reached without any 
leakage from the seal. 

6. The seal life diagrams for different motion speed values 
indicate that for each of the fluids can be observed a range of 
speed values below which a significant extension of failure 
— free operation period of the seal occurs. The process is 
probably associated with Kelvin— Helmholtz phenomenon 
of instability [2, 3, 11]. 

7. Value of statical penetration pressure depends on the 
saturation magnetization intensity of magnetic fluid, M,. 

8. Magnetic field intensity in the seal system only slightly 
influences the increasing of statical penetration pressure 
values of the seal. In the tests with the use of 18 permanent 
magnets, only a small decrease of the penetration pressure 
was observed in comparison with that obtained in the tests 
with 14 permanent magnets. 

9. Magnetic field intensity has a low influence on a rotational 
speed value at which the basic value of seal life equal to 1 
min rotations is obtained. 
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ABSTRACT 


In this paper the present state of grouping technology theory is discussed. Attention is 

paid to lack of publications on grouping methods based on linear operators instead of 

classical incidence matrices. A computation algorithm for generating three-dimensional 

groups has been proposed. Every group contains joined structures, working cells and final 

(output) structures. Operation of the algorithm has been compared with results gained from 

literature sources dealing with binary grouping. And, calculations based on real assembly 
sequencing of a set of ship hull sections, have been performed. 


Keywords: group technology; assembly sequencing; steel structures manufacturing 


1. INTRODUCTION 


Hulls of typical merchant ships belong to a group of steel 
stiffened shell structures (3S). Such structures can be found also 
in very large tanks, bridge elements, offshore units etc. Because 
of their dimensions and necessity to expand working frontage, 
production of 3S structures is usually carried out in accordance 
with the so called prefabrication subdivision. It means that 
final product is composed of many simultaneously produced 
prefabricated units (subassemblies). They usually are further 
subdivided into lower stage units. Consequently, a multistage 
assembly sequencing is formed. Number of stages depends on 
structure geometry and manufacturing conditions of a producer 
(a. o. — cranes at his disposal). 3S structures are as a rule 
produced in a few pieces or short series, that forces enterprises 
such as shipyards to keep agile production organisation and 
makes designing and production preparation work more labour 
consuming. The problem of a degree of producing in series is 
tightly connected with the notion of technological similarity of 
products. Products even if different visually and functionally 
may require similar production processes and — from the 
enterprise point of view — belong to one common series. 

One of the cardinal tasks of design office is to act in 
accordance with the principle of ,,design for production” 
by applying as much standardization of design solutions as 
possible, including assortment of rolled profiles and shapes 
of cut details. Such actions are of a low degree of structural 
complexity. 

A successive action aimed at minimization of unitary 
production character of 3S structures is searching for such 
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division planes for their assembly as to obtain the effect of 
a greater degree of serial production at least in early stages 
of assembly process. Two 3S structures may differ from 
each other in technology, for instance: by accessibility for 
automative welders, shapes of curvatures, assembling sequence 
and possibility of being manoeuvered. Simultaneously, their 
structurally simpler subassemblies may be very similar 
stiffened panels as regards their technology. 

The theory which makes it possible, on the level of models, 
to: 

- measure technological similarity, 
- state when a given product can be qualified to be an element 
of greater series, as well as 
- so design a given production system as it could realize 
production in a way as repeatable as possible 
is the classical theory of group technologies (GT). Its main aim 
is to distiguish, in the set of resources used in the process or 
set of products, some subsets called groups. In the literature 
(2, 3, 4, 5, 13] only working cells of different kinds or, simply, 
machines are usually distinguished as grouped resources. In 
the case of grouping the products we aim at maximization 
of their repeatability and reduction of design and production 
preparation processes [11]. Due to having at hand certain groups 
of products, designers of prototypes have at their disposal hints 
helping to maintain high producibility. 

Among profits resulting from application of group 
technologies to organizing production cells the following are 
numbered [9]: 

- widened standardization of working cells, applied 
instrumentation, set-ups of working parameters, 


- reduction of interoperative transport as a result of relocation 
of work objects within production cells, 

- facilitated scheduling of production, 

- shortened cycles of production, 

- lessened buffer stores and lowered freeze level of current 
assets, 

- elevated motivation of employees and simplified procedures 
of worktime estimating. 


The grouping is performed with use of: 

- visual assessment of elements of a set — on the basis of 
a subjective validation by experts, 

- classification and coding — on the basis of an assumed 
coding system of features (in the case of products: their 
geometry and production technology), 

- analysis of technological routes realized by products. 


The method of analyzing technological routes is advantagous 
as it makes it possible to group simultaneously a set of resources 
and products. Here the groups of the set of resources are called 
input ones. 

Identification of input and output groups permits also 
to design the so called production cells where each of them 
constitutes the pair (group of resources; group of products). 
Such pairs are called two-dimensional. Algorithms for two- 
dimensional grouping have been developed for a few dozen 
years [6, 1, 16]. 

The grouping of factors by using the method of analyzing 
technological routes is realized according to the following 
scheme [16]: 

1. The formulating of relations between resources and products 
by using the so called incidence matrix — to each of the 
resources an appropriate row of the matrix is assigned, and 
to each of the products — a column of it (or vice versa — it 
is a question of convention). 

2. The determining of mutual technological similarities between 
resources as well as between products. Consequently, two 
matrices of technological similarities are built — separately 
for each of the analyzed sets. 

3. The setting ofa threshold for technological similarity, above 
which it has to be considered that two elements of the set 
in question should be assigned to the same group. 


The two-dimensional grouping makes it possible to respond 
quickly to the question: how to organize production cells in 
the system if planned assortment of products and engineering 
techniques available to producer, are known? It is easy to 
reach the effect if the incidence matrix unambigously describes 
relations between given products and working cells. However 
in practice the forming of such incidence matrix is very labour- 
consuming and demanding comprehensive knowledge on 
technologies used in a given enterprise. 

An extension of the classical two-dimensional grouping 
method is the proposal of distinguishing, in the set of used 
resources, two separable subsets. As aresult, we have to do 
with three sets and perform analysis of similarity between 
products on the basis of their incidence with the first and 
second set of resources. In both the sets of resources we 
analyze similarity independently, on the basis of incidence 
with the second subset of resources as well as the set of 
products. From the analysis result three-dimensional groups, 
i.e. the ordered threes: (1s group of resources, 2"4 group of 
resources, group of products). The first concepts of the kind 
belong to Min and Shin [17]. Their proposals deal with analysis 
of technological similarities in the two subsets of resources: 
machines and operators as well as in the set of products. 
Successive elaborations have been a natural generalization 
towards n-dimensional models [14, 15, 18]. 


2. CONCEPT OF THREE-DIMENSIONAL 
GROUPING IN MANUFACTURING 3S 
STRUCTURES 


2.1. Input data, standardizing and scaling 


The three-dimensional grouping is proposed for application 
in manufacturing 3S structures with taking into account an 
ordered set of various working cells in which joining process 
of structures into higher stage structural units, is performed. 

The following notions are introduced: 

- S= {Si Sm — working cells, 


= Eie d 


- CTS ot — produced structures (put out the 
cells). 


— joined structures (put in the cells), 


The structures differ to each other as regards their geometry. 
And, it is not essential where butting contacts of plate sheets 
are placed, nor other features of technological importance 
(e.g. form of cut-outs and gasket plates). Only the functional 
features are taken into account, namely: arrangement of shells, 
stiffeners and brackets, ways of ensuring continuity of girders 
etc. Structures of various complexity degrees are considered. 
It may be single details or large-size spatial structures, e.g. hull 
block sections or even entire ship hull. 

To every output structure can be assigned elements of 
the set S and C which take part in its production. To every 
working cell can be assigned structures which are delivered 
to it. It leads to formulation of three matrices of incidence 
between relevant sets. The concepts have been published so 
far are based on binary matrices. However production of 3S 
structures is associated with a rather low diversity of cells. And, 
differences between structures amount to labour consumption 
of processes as well as quantity of lower-stage subassemblies 
of which the structures are formed. Practical implementation 
of grouping technologies into production of 3S structures 
faces the problem of an insufficient informing force of binary 
incidence matrices. For instance, structures which greatly 
differ to each other as regards requirements for equipment of 
working cells, can be built of identical subassemblies and with 
the use of identical welding methods. In this case, application 
of the binary grouping methods may lead to erroneous placing 
technologically different products into the very same cells. 
As a result low specialized production cells are obtained, that 
contradicts their essence. 

As asolution for the above mentioned problems it is 
proposed to apply a grouping model based on continuous 
incidence matrices. For every pair of sets of the input factors 
X and output factors Y the following matrix is determined: 
AYe M(XI x IYI, R*U {0}, where aX-Y > 0, if production of 
j-th element of the set Y requires to angage i-th element of the 
set X, ifnot: aX-Y = 0. Values of the non-zero coefficients a*-Y 
are standards of consumption of input factors per production 
unit of appropriate output factors. In the case of working cells 
the consumption standards usually represent the following: 

- amount of working time of a cell during production of unit 
of a structure, 

- amount of a concrete effect of work of a cell, generated 
during production of unit of a structure (e.g. one meter of 
welded joint, one square meter of painted surface etc). 


The continuous incidence matrix is a linear operator acting, 
in this case, between vector spaces of production factors 
[12]. Such operators find wide application to a.o. operational 
research, and they may find much wider use in solving the 
grouping problems in enterprise. 
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2.2. Measure of grouping quality 


Three-dimensional grouping the sets S, Ci and Cet is an optimization problem and amounts to searching for the division T 
of the Cartesian product S x Cin x Cout into y subsets: 


P= (Sx? xc™,7\= 


= lr, sS esx x0" imar), ay = fc Z A Com = 2} 
i=l i=l i=l 


The groups T; are really sums of the subsets S,, C= and Cout, however they are written in the form of ordered threes, that 
facilitates to identify elements of each of the group. 

The division I’ can be assessed as regards its quality. In 1986 Chandrasekharan and Rajagopalan [7] proposed to apply 
a criterion for the two-dimensional grouping of binary incidence matrix. This is maximization of a weighted sum of densities of 
submatrices and their neighborhood. And, the density is meant to be number of unities related to whole number of elements of 
the matrix in question. In 2007 Li [15] extended the criterion, for the multi-dimensional case, by assigning weights to particular 
incidence matrices. 

In the case of continuous matrices, differences in values of the elements should be accounted for additionally apart from the 
fact of occurrence of non-zero elements of incidence matrices in the space of groups and outside it. 

Let the weights of particular incidence matrices be equal to: @, ¢., @s_¢ respectively. Then it is proposed to take, 
as a quality measure for I division, the following criterion: . 


(2.1) 


We. 
out? Cin = Cout 


pout 


iS cn S-C cin cu 
Osc, WAT) +056, WAST) +06 6, WAT) 


E ( r) _ — max (2.2) 
Osc, + Osc, +O, -c.,, 


where: 
W(AX-Y) — grouping quality index for continuous matrix of incidence between the set X and Y. 


Let the set of indices of elements of the set X belonging to k-th group is marked i Then the following is defined: 


W( AT) =9-n (4°) +(1-9)-m (4°72) (2.3) 


where: 

n,(AX-¥, T) — sum of values of elements of the matrix AX-Y belonging to all groups, related to their total number, 

n(AX-Y, T) — sum of values of elements of the matrix AX-Y not belonging to any of the groups, related to their total number, 

q — weight of membership of elements in groups in relation to their possible being left outside the groups. In the 
literature q = 0.5 is often assumed [15]. 


n(4"7,r)=S Se et 2.4) 


DIDI 


k l#k iel jel} 


m| A“ T) = SS (2.5) 


Deel 
k læk 


The quantity n,(AX-Y, T) is the mean value of all elements within groups, and n,(AX-Y, T) is the mean value of elements 
outside the groups. As the elements a * are normalized within the range [0,1], both the quantities are also contained in this 
interval, consequently, the inequalities: 0 < W(A*-Y, T) < 1, 0 < E(I) < 1 are also satisfied. 

Additionally, every group is assumed to be consisted of sequence of non-empty sets. For instance, to the groups cannot be 
assigned only a set of output structures and cells without a set of input structures. The assumption is aimed at preventing against 
natural grouping input structures into one group, cells into second group and output structures into third group. 


2.3. Algorithm of the three-dimensional grouping for continuous incidence matrices 


The algorithm of multi-dimensional grouping was presented by Li in 2003 [14]. However this proposal is limited to binary 
incidence matrices. As a result it does not make use of cluster analysis methods but is the author’s solution strongly based on 
the binary description of elements of sets. 

Below are presented principal steps of the proposed algorithm. It is of a recurrential character. The grouping is realized by 
determining centres of groups having assumed number of elements and assigning to them the closest elements according to 
a selected measure. 
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The requirement of that every group should contain 
elements from all the grouped sets is satisfied in Step 2.5, in 
which one element of each of the sets is assigned to each of 
the groups and only in the next step unlimited assignment of 
remaining elements takes place. 

As already mentioned, it is worth normalizing the continuous 
incidence matrices within a selected interval. 


Start-up of the algorithm 
Step 1: 
Having at our disposal normalized continuous matrices of 
incidence between the sets: S, C! and Cot we form the following 


integrated matrix: 


~out 


(2.6) 


(m+n+p)x(m+n+p)} 


Further proceeding consists in solving the problem of 
two-dimensional grouping, where the cells and input and 
output structures are considered jointly as the so called set of 
elements: 

in out 
Peso UCT s fesla) CO 

The elements numbered from 1 to m represent working 
cells, the elements of numbers from m + 1 to m + n - input 
structures and the elements of numbers from m + n + 1 to m + 
n + p - output structures. Therefore we have: 


fi = Si N AO = Sino 
Jana m R A = Gs 


_ nut f. — nut 
Sed Z0 opes, m+nt+p Cp 


(2.8) 


Step 2: 

The division K of a set of elements F into clusters (one- 
dimensional groups) such that one working cell, one input 
structure and one output structure is present in each of them: 


K(F)= {K,cF: 


i=l,..„min{ i en 


l, (2.9) 


W(K, AS +Ø, K, AC" +Ø, K, nC" # 2) 


The division K should be performed for the following 
number of clusters: 


í c” ‘ (2.10) 


cf 


Such division for a given y is marked K”. Successive steps 
of the algorithm should be made for each of the variant of y 
value. 

The division into clusters can be realized by using one 
of the divison-based algorithms. Here a modified c-means 
method is implemented. The method in its classical form has 
been described a.o. by Hartigan in 1975 [10] as well as Tan, 
Steinbach and Kumar in 2005 [19]. The proposed modification 
is aimed at ensuring occurrence, in each of the clusters, at least 
one element out of each of the sets: S, Cin and Co, 


y=2,...,min| h 


Step 2.1: 

The forming of the cluster division matrix A” = (0) 
q=m +n +p is total number of factors. 

In its successive steps the algorithm assigns value 1 (unity) 
to components of the division matrix at intersection of cluster 
columns and rows of relevant factors. The initial division is 
made in accordance with the following algorithm: 


_| mene 


4 
for| j =1,2,....7 —1]: 


for[i=1+(j-1) r,s jer]: Ay <1; 
for[i=m+1+(j-1) r, m+ jer]: A% e1; (2.11) 


axy Where 


//rounding down 


melfamnete ii ment j-r]: Ai <1; 
for[ i =1+( (y-1)- 


ra 


att |: Ai, © l; 
e A el; 


for[i=m+n+1+(y-1) ryd |: A <1, 


return ^”; 


Step 2.2: 


The determining of the matrix of cluster centres, (a)... 


(2.12) 


Step 2.3: 

The determining of the matrix of distances between 
elements and cluster centres. 

For every i-th element and j-th group, D; value is 
determined. Here, it is possible to assume an arbitrary measure, 
not necessarily metric one. The following variants are proposed 
to be considered: 

- angular measure (scalar product of vectors): 


AT -a 
yv Y D= £ 
a j=ly y q 5 q (2.13) 
DY Ay Ay 
k=l k=l 
- Manhattan measure: 
v Y D, =|4°,- a) |+...+|40 | 2.14) 


i=l.. JH. 


- ties measure: 


al’ j (2.15) 


q 


v v D,= max Í 


izl, g Jal 7 kalno 


AM, -a eio 


- Mahalanobis measure: 
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In the below presented exemplary calculations the above specified measures have been mutually compared as regards their 


effectiveness. 


Step 2.4: 


The saving of the cluster division matrix after introducing a new variable: A = Av. 


The zeroing of the cluster division matrix: Ay = (0),..,. 


Step 2.5: 


The assigning of one factor from each of the three subsets to every cluster. This way three-element clusters y in number 
were built. Below is presented the algorithm of assigning the factors numbered from | to m, i.e. these belonging to the subset 


of cells: 


for [i =1,2,...,m]: 
for | j =1,2,...,7]: 
for |! =1,2....,7]: 


A, — D,; 


y 


for[i=1,2,...,7]: 
for| j =1,2....,7]: 
if[ A, =d |: 


Ay & l; 


//forming local matrix of distances 


(2.18) 


//assigning object to cluster 


AY) (max {A} +1 ; //ocal blocking j-th cluster (column) 
mxl 8. 


AM) eo 
De 
return A’; 


For the two successive subsets the procedure runs 
analogously. Changes deal only with the range of variability 
of factor’s number in the loops for[i = ...]. 


Step 2.6: 

In the preceding step it was ensured that in each of the y 
clusters one working cell, one input structure and one output 
structure was placed. Now it is possible to assign factors, 
without any limitations, by taking into account only the criterion 
of minimization of distances from the cluster centres a®: 


for[i =1,2,...,q]: 
for[ j} =1,2,...,7]: 


it| D, = min {D™} |; AY =1; 


return A’; 


Step 2.7: 

Comparison between the obtained cluster division matrix 
and the solution reached from the preceding loop. Here 
Frobenius norm is applied to difference between the compared 
matrices: 

i 


A 


Jaz -Â 


DAE ~ hy y 


i=l j=l 


(2.20) 
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(max {A} + I) ; //local blocking i-th object (row) 
(max {D} + I) ; //global blocking i-th object (row) 


Further action is dependent on assessment of value of the 
norm . The assessment is made by comparing its result with 
the assumed limit value 6 experimentally determined for 
a given case. 


ir [| 


otherwise: go to step 2.2; 


Step 3: 

For the determined y the cluster division matrix A’ is 
known. The clusters K{,..., K} are also given in the form of sets 
of the factors to which correspond unity values in respective 
columns of the matrix A’. As a result the division T into three- 
dimensional groups is obtained: 


V T, =(K7S8,K7 0C",K7 NC™) (2.22) 
i=l... 


Step 4: 

Comparison of variants of number of clusters and applied 
distance measures is performed by determining the division 
which maximizes value of the objective function : 


l=: E(T)— max 


The end of the algorithm 


A . . 
A Al < Ji go to step 3; Gm 


(2.23) 


The investigating of many variants of number of clusters 
may seem time-consuming. However it should be taken into 


account that in practical building 3S structures we deal with a total number of different working cells, which amounts to a few 

dozen only. This is much less than that of possible variants of input and output structures. 

The algorithm makes it possible to search for similarities between produced 3S structures, working cells and structures which 
are joined together in the cells. As a result, is obtained information which deals with the following: 

- structure of working cells, i.e. which operations have to be realized by the cells and how many working cells should realize 
given operations, 

- kinds of the structures which will be most often delivered to the cells to be joined together there, which is connected with 
ensuring appropriate transport means, as well as with possible location of the cells and manufacturers of subassemblies in 
close neighbourhood, 

- indication of a cell which has to be responsible for production of a given structure in the case when the structure is qualified 
to belong to one of the groups. 


The described problem may seem to be solvable by independent realization of the two two-dimensional grouping: output 
structures — cells, input structures — cells. However it should be observed that possibility of building a given output structure 
by a production cell is strictly dependent on a degree of complexity of input structures which have to be joined together. Two 
cells may produce the same structure, however one of them may be specialized in assembling complex subassemblies whereas 
the other may realize the process of assembling single details to form simple structural units. Therefore the output structure will 
be assigned to three-dimensional group only if it itself shows technological similarity and also the input structures and working 
cells coincidental with it do the same. 


Technological database: Production task: 


Task assigning to one 
of the groups 


Database 
updating 


A 


Organization of production 
‘ cells 


x | Production 


Fig. 2.1. Location of the proposed algorithm into production planning process of 3S structures 


3. EXEMPLARY CALCULATIONS 
3.1. Example 1 — the binary problem proposed by Li 


In 2007 Li [15] considered the following 3-D binary problem (let’s call it: „Li” problem): 


P, P, P, P, 0, 0, 0, 0, 0; 0, 0, 0, 0, 9, 
m (1 1 m, 1 1l p,{l 1 1 
AYP =m, | 111ļ|,4 °? =m |1 1 1ļ},4°=p]| 1 11 
m, 1 m, 1 1 pil 1 1l 
m, 1 1 m, (1 P, l 1 


and formulated its pseudo-optimum solution in the form of division into the two groups: 


r, =({m,,m,},{P,, Ps} .{0,,0;,05}) 
D, =({m >m}, {Pi Ps} +{02+04}) 


The grouping can be represented by the sorted incidence matrices with the distinguished submatrices of groups as follows: 


P, P, P, P, 0, 0, 0, 0, O, 0, 0, 0, 0, 9, 


m,{|l 1} 1 m, (|111 p, 1] 11 
AM? =m, | [ll AN Oo =m, | LI Aon =P, Ij 11 


m, 11 m, 11 P| 111 
m, 1 m, 11 P, | 111 
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The quality index of the above mentioned grouping, under assumption of the weights q = 0.5 and œ; = œ, = @,, amounts to: 
ETY) = 1/3(0.875 + 0.95 + 0.05) = 0.625. 


The Li’s solution is compared now with the results obtained by means of the proposed algorithm and distance measures 
described by Eqs. through. As can be observed in this case the division can be done into 2, 3 or 4 groups. Tab. 3.1 presents the 
divisions which provide the best values of the quality indices for particular measures. 


Tab. 3.1. Results of solving the problem ,,Li” 


Measure Division 


Manhattan 


Euclidean, {m,}.{P,P2»P3}.{0})s 


ahal: ig l f 
Mahalanobis m,,m,,m,},{ p,},{0),0;,0,,05}) 


Angular 


}), 


Chebyshev 


As can be observed, the Chebyshev measure provides the division of the quality index higher than the solution obtained by 
Li. The division can be represented by the following sorting of incidence matrices: 


P, P, P, P, 0, 0, 0, O, 0, 0, 0, 0, O, 0, 
m, 1 1 m, l 1 P, l l 1 

Amn =m, || 111|, A =m, Ann = P, il 
7 l m, p, | {111 


| 
m, 1 [] m, P, 


It can be hence stated that the proposed algorithm is potentially useful to grouping the sets whose incidences are expressed 
by binary variables. Certainly, quality of obtained results depends on assumed weights. Let’s observe that the increasing of 
the weight q would cause the increasing of superiority of the solution obtained with Chebyshev measure over the result of Li. 


Simultaneously, the decreasing of the weight œ, would result in improving the solution of Li and decreasing the quality of our 
solution. 


Example 2 — extension of the binary problem of Li 


In order to present action of the grouping algorithm in the case of continuous incidence matrices the binary problem discussed 
in Example 1 will be modified. It is done by assigning random values from the interval (0,1] to non-zero elements of incidence 
matrices. As a result the following problem is obtained: 


P, Pr P; Ps O, 0, 0; 0, 0; O, 9% 0, 0, O; 
m, (01 0.5 m, 0.2 0.3 p, (0.2 0.1 l 
A”? =m, 021 09|,A”° =m, |1 07 05|,4°° =p, 0.9 0.3 0.8 
m, 0.6 m, 0.3 0.1 p, | 0.5 07 0.6 
m, 0.5 0.3 m, \ 0.6 P, 0.9 0.4 


Tab. 3.2 presents the best results obtained with the use of the grouping algorithm. 
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Tab. 3.2. Results of solving the problem ,,Li” extended to continuous form Also in this case the Chebyshev measure turns out to be 


Quality most effective in grouping the elements. This time the division 
into four specialized groups is the best. In matrix form it can 
be presented as follows: 


-0 
AM =m, |] me > 
Manhattan j > 5 m, [1] 0.7 0.5 
s ™% {o3 o [| 
0, 0, o, 0; 0; 
Mahalanobis 2 2 i DP 0.3 0.8 
P-O 
| AP? =» [] 02 011 
Pa| o9 04 [| 
Chebyshev 


Let an example of building five kinds of input structures shown in Fig. 3.1, be considered. The structures have been not 
subjected to any analysis of possible, at least partial, unification of their sets of lower-stage subassemblies. Intentionally, real 


rough data are considered. 
2 
5 


Fig. 3.1. Considered output structures 


3.2. Example 3 — practical problem — the welding of steel structures 
3.2.1. Problem description 


1 
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On the basis of technical documentation of subdivision of the output structures into prefabrication steps, a list of 37 input 
structures (subassemblies) was elaborated. They will be joined together by welding processes. The continuous incidence matrix 


Act. c™ after its normalization within the interval [0,1], takes the following form: 


out out out out 


C c C 


1 5 1 5 
c” [0.33 0 000 ci” 10 0 0.04 0 0 c" cs 
l 0 0.04 c |0 0 0 0.04 0 
0.33 0 0.04 0 0.15 
0.15 0 0.04 0 0.11 
0.04 0 0.19 0 0.15 
ge- | 0.04 0 0.19 0 0.04 
0 0.04 0 0.19 0 0.04 
0 0.04 0 0.19 0 0.19 
0 0.04 0 0.19 0 0.04 
0 0.04 0 0.52 0 0.19 
0 0.04 0 0.52 0 0.04 
0 0.07 0 0.19 ci” [0 0.33 
c | 0 0.04 c”, |0 0.04 


The third set of elements is formed from welding operations which differ to each other as regards welding position (Fig. 3.2), 
type of weld (B — butt, T — fillet) and accessiblity for welding machine (M — semi-automatic, R — automatic). 
Not all combinations of the features take place hence the set S has eight elements only — see Tab. 3.3. 


Element of 


meats Position 


Type of joint 


Tab. 3.3. Symbols of welding operation types 


Degree of 
automation 


Fig. 3.2. Welding positions: O — overhead, V — vertical, H — overside, 
horizontal, D — floor 


On the basis of an analysis of output structures, lengths of the welds made in the distinguished welding positions, were 
determined. Also, for each ofthe input structures lengths of the welds for its joining with other input structures, were determined. 


As a result the two successive normalized incidence matrices were obtained: 


E H Co a 
s (0 0 0 0000 0 0 0 0 0 0 033 0 
0 0.05 0.03 0.02 
. 0 
AS = | 0.01 0.02 0.01 0.18 0.02 0.05 0.01 0.04 0.03 0.03 
0 
0.01 0.01 0.01 0.12 1 
0.01 0.01 0.01 0.12 1 0.24 0.03 0.02 0.06 0.06 0.02 0.03 
S \ 0 
Cis ĉj Č% G Cp C3 
s, 0.02 0.03 0.03 0 0 0 0 0 0 0 0 0 0 0 0 0 
0 0 
0.03 0.04 0.03 0.04 0.02 0.05 0.04 0.02 0.01 0 
0 0.003 0.01 0.02 0.01 
0 0.06 0.03 0 
0 0.02 0.03 0.03 0.29 0.02 
0 0.1 0.003 0.02 0.01 0.19 0.32 0.02 
S 0 0 
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in 
17 


0.06 0 


C. 


0.14 0.18 


0.36 0.19 0.03 0.06 


0.1 0.02 0.08 0.01 


0.06 0.02 0.03 
0.06 0.02 0.09 0.02 


„out „out 


s,(0 0 017 0 0 
0 0.02 
0 1 

AS" = |029 0.04 0.01 0.08 
0 0.06 
0.41 0.2 0.13 
0.41 0.1 0.04 0.2 

s \ 0 0.004 


3.2.2. Results of the grouping 


With the use of the algorithm described in Chapter 2.3 the 
divisions into 2, 3, 4 and 5 three-dimensional groups were 
obtained according to the condition. The following values of 
weighing coefficients: q = 0.5, @s -c = Os _¢ 
were assumed. 

Five measures of distance between elements were tested, 
i.e. the angular, Manhattan, Euclidean, Chebyshev and 
Mahalanobis one. For four numbers of groups and five distance 
measures, twenty divisions and quality indices corresponding 
to them, were obtained — Fig. 3.3. 


=O. = 
out Cin- Cout 


Manhattan 
Euclidean 
Chebyshev 


Mahalanobis 


Number of 
groups 


N 


Measure 


Fig. 3.3. Quality indices of divisions for different numbers 
of groups and distance measures 


The best solution was obtained for the division into five 


groups by using Euclidean measure: (E(T) = 0.545). The groups 
are as follows: 


in „in „in out 
L; =({s,}.{e8 6 ee }) 
in in in in in in 


in 
T,= (ORG KAAT R A A 


in in in in 


Ea ra T) 
D (sheda) 
Pas (besesr fee fer) 
T; = Maea) 


If the grouping between pairs of the sets are taken into 
account then it can be observed that the greatest mean value 
in the groups (Eq. ) is always reached for the division into five 
groups. However effectiveness of the measures shows different 
values for different pairs, e.g. for: 

- §-Cin — angular measure is the best, (n,(AS- ©", T) = 


= 0.0515) 

- §-Ci— Euclidean measure is the best (n,(AS-°", T) = 
= 0.2857), l 

- Cin-Cout— Chebyshev measure is the best (q, (A-0, T) = 
= 0.1131). 


As can be observed, during analyzing we are forced to take 
into account various measures, since their effectiveness depends 
on values of co-ordinates of grouped elements. 

In the above discussed example can be observed a tendency 
for breaking up the division into as many groups as possible. 
It results from a low density of incidence matrices which 
constitute input data. And, such density directly results from 
the breaking-up of the set of input structures. Reduction in size 
of the set Ci" is possible by standardizing the input structures. 
Influence of such approach on results of grouping will be 
discussed in a separate publication. 


CONCLUSIONS 


In this paper has been presented the algorithm which on 
the basis of three assumed sets, groups their elements in such 
a way as locate at least one element out of each set into each 
of the group. In the algorithm the method of cluster analysis of 
k-means as well as five different measures of distance between 
multi-dimensional space points, was applied. It was elaborated 
with the aim of its applications to organization of production 
of 3S structures. To this end, were distinguished three sets of 
elements occurring in enterprises which produce structures of 
the kind. The sets as well as continuous matrices of incidence 
between them form input data for the proposed grouping 
algorithm. These are: set of working cells, set of input structures 
(subassemblies) and set of output structures. 

The grouping division obtained this way can find application 
during designing production cells. Certainly, the above 
determined groups are only an indication as to mutual location 
of working cells and their equipment. Three-dimensional group 
can be interpreted as an indication for very close location of: 
- working cells for production of concrete input structures 

(or points of delivery of the structures from outside of the 

enterprise), 

working cells for joining the structures into larger objects 

(output structures) as well as 
- working cells which use the output structures for further 

joining them, or possible cleaning, painting, fitting, 

launching or also loading and dispatching. 


Full modification of a non-optimum spatial organization 
of working cells may be impossible. It is hard to expect that 
because of results of grouping whole workshop buildings would 
be moved. However enterprise logistic system is much more 
adaptable in this respect. 

The group may be hence interpreted also as an indication for 
elevation of enterprise investment priorities in certain areas. 

It can be: 

- extension of internal roads between selected cells; 

- purchase of transport means dedicated to certain structures, 
acc. their mass, gabarites, shapes; 

- ordering, to subcontractors, some workings which especially 
disturb determined production flows within enterprise. 
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pe 


A general conclusion which can be drawn on the basis of the 
rformed calculations is that the problem of designing group 


technologies becomes greatly complicated in the moment of 
relating it to real production problems. Theoretical methods in 
this area have been developed for tens of years however reality 
seriously challenges both theoreticians and engineers aimed at 
implementing the methods into practice. 


BI 
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ABSTRACT 


The process improvement operations are very significant in shipbuilding industry as the other industries. 
In recent years, the shipyards attempt to improve their processes by examining their current production 
system and to reduce the cycle time of the interim product so that they can keep their competitive power. If 
the cycle time of the interim products is decrease, that may cause to increase the annual production capacity 
and market share of the shipyard. In order to do this, the shipyards have to analyse their own production 
system and carry out some improvements on it. In this study, the effects of carrying out the outfitting and 
piping assembly operations in earlier work stations were investigated by using a methodology presented 
here. The steps in the methodology were applied to a double bottom block ofa container ship. It was shown 
that if the outfitting and the piping assembly operations are carried out in earlier stations instead of the 
block assembly station, this may increase the throughput by 33%. The results of the study were discussed 
in the final section. 


Key words: Pre-outfitting assembly; pre-piping assembly; simulation; ship production 


INTRODUCTION 


Shipbuilding is a global industry competing in the world 
extent [1]. In recent years, the shipyards attempt to improve 
their production processes in order to yield advantages against 
their competitors by manufacturing the ship as soon as possible. 
Geoje Samsung Shipyard is able to manufacture 40 ships per 
year in South Korea and it is also known as one of the most 
efficient shipyards in the world [2]. As compared with the other 
shipyards, Samsung shipyard has a great competitive advantage 
against its competitors. So, the point is that the shipyards have to 
investigate their production system and do some improvements 
in order to keep the competitive power. 

One of the latest trend in manufacturing process is to reduce 
the cycle time of the product. In order to reduce the cycle time, 
it is needed to make process improvements on the current 
production system. Eker [3] reported some improvement 
suggestions with regard to the processes in the painting shops 
of the Sedef Shipyard and Schalekamp Shipyard. In the study 
[4], the format of the design drawings coming from design 
department was changed and the cutting process was carried 
out in shorter time and in more effective way. The process 
improvement operations was also applied on the scaffold area 
and material stock area [5] and the layout of the work stations 
[6] and [7]. In Todd Pacific Shipyard, as a result of improvement 
process, the moving distances of the forklifts were shortened 
in the rate of 50% and the oil wastes occurring during the 
operations of the forklifts were also reduced [8]. As can be seen 


from the above works, the process improvement applications 
are able to be applied in many fields in shipyard. 

The effects of the improvements on the production system 
can be seen by simulation. Simulation has a great importance 
for the production companies. In the competition environment, 
the changes on the production system and the effects of 
these changes are very significant in terms of the company 
performance. Simulation has been applied in many industry 
for years as it provides a great advantage for the planner. In 
shipbuilding industry, it has a great deal of application fields 
such as layout, production processes and so on. In the study 
that has been still going on by Michigan University and Seoul 
University, it is aimed to simulate all the activities in shipyard 
and to see the effects of the changes on the production system 
[9]. In the other study, the stations forming the sub assembly 
line were modelled by using simulation and after the system 
was simulated, a robot was settled in the production line and 
the rate of productivity was determined [10]. Shin [11] aimed 
to settle the work stations of the shipyard in the optimum 
way. That’s a layout application of simulation. Alkaner [12] 
considered a profile cutting station and the processes of the 
profile cutting station were determined. Then, these processes 
were modelled in a simulation program and the effects of 
changing the resources were investigated. In the study [13], 
the panel production station was considered as a bottleneck 
station. Processes were determined and modelled in simulation 
program. Then, by doing some changes on the processes, the 
completion time of the panel cutting station was tried to be 
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optimized. In this study, Arena simulation program [14] was 
used for modelling the work flows between stations as the 
modules of it are very appropriate for modelling the production 
activities. The activities can also be modelled easily by using 
the modules. 

In many shipyards in Turkey, the whole assembly operations 
with regard to outfitting and piping are carried out in the block 
assembly area. In the block assembly area, after the steel 
assembly operations of the block are over, the outfitting and 
piping assembly operations are carried out. Therefore, the work 
load of the block assembly area is very high and this may cause 
to a new bottleneck in the production system. In this study, the 
outfitting and piping assembly operations of the block will be 
carried out in earlier work stations and some changes will be 
applied on the current production line. As a result, the effects 
of all these changes on the current production system will be 
determined. The phases of the methodology are followed one by 
one and the effects of the pre-outfitting and pre-piping activities 
on the entire production system are seen. 

In literature, the papers with regard to pre-outfitting are 
very limited. In the papers, the modular outfitting concept 
is mentioned. Fafandjel [15] mentioned that the ship is 
manufactured in shorter time by using modular outfitting 
concept. Yu and Ishida [16] searched how to determine the 
effectiveness of the modular outfitting concept. Baade et al. 
[17] mentioned the advantages of the modular outfitting in 
his paper. 


METHODOLOGY 


Fig. 1 presented the methodology implemented in this study. 
The methodology consists of nine steps. In the first step of the 
model, the product, which is produced in shipyard production 
line, is defined. Then, the workstations, that are involved 
in product fabrication, are described. In the step 3, detailed 
process analysis of current production system is performed. 
After comprehensive process analysis, simulation model of the 
current production system is modelled in ARENA simulation 
software (in step 4). In step 5, simulation model is run along 
a specific period and production quantity (throughput) of the 
production system is achieved. Then, some suggestions on 
current production system are made in step 6. After that, these 
suggestions are applied on current simulation model in step 7 
and the new production case is obtained. In step 8, the effects 
of the changes are discussed and evaluated. Finally, in the last 
step (step 9), comparison of the current and new production 
cases is performed. 


APPLICATION 
Defining of product (Step 1) 


Ships are fabricated in blocks which constitute ship 
structure. There are various sort of blocks such as bow, aft, 
side, bulkhead and double bottom block. In this study, a double 
bottom block of a container ship is illustrated as an example. 
Double bottom block is one of the main interim products in 
shipbuilding. The reason of considering the double bottom 
block in this study is that almost all of the work stations are 
involved in its production activities. 

As the structure of a double bottom block is examined, it can 
be seen that it has different production phases, which are called 
as production stages. Table | shows these production stages and 
the structures representing these production stages. 

A double bottom block is built by coming together the 
production stages mentioned in Tab. 1. In the first phase of the 
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Defining of product (Step 1) 


Description of work stations (Step 2) 


Detailed process analysis of current 
production case (Step 3) 


Simulation model of current 
production case (Step 4) 


Determination of throughput of 
current production case (Step 5) 


Suggestions on the current case (Step 6) 


Application of the suggestions on the 
current case (Step 7) 


Evaluation of the new case (Step 8) 


Comparison of the current and 
new cases (Step 9) 


Fig. 1. Application steps of the study 


double bottom production, single section parts (A) and single 
plate parts (B) are fabricated. These parts are cut from the plates 
and profiles which have standard-dimensions and they have 
specific dimensions after cutting process. Then, they are fitted 
and minor assembly is built (C). Two or more minor assemblies 
constitute sub assembly (D). Flat plates are sub-merged welded 
and flat plate assembly (E) is fabricated. When the profiles are 
fillet welded on flat plate assembly, the flat plane assembly (F) 
is built. Minor and sub assemblies are welded on the flat plane 
assembly (F) and major sub assembly (G) is manufactured. 
Curved panel assembly (H) is manufactured on pin jigs. In 
block assembly area, sub unit assembly (J) and curved panel 
assembly are mounted and welded, finally a double bottom 
block (K) is built. 

There are also some outfitting equipments in double bottom 
block just like manholes, bottom plugs, zinc, vertical ladder 
and doubling plates. These outfitting equipments are mounted 
into block structure’s body. Tab. 2 shows these outfitting 
materials. 


Descriptions of work stations (Step 2) 


Ship production is extremely hard job since it includes 
a great deal of process. A ship is manufactured by performing 
thousands of work activities. In order to manufacture a ship, 
various types of workstations are needed. Every workstation 


Tab. 1. Production stages of a double bottom block and structures 


aion | ard a a Structures representing production stage 
Stage stages 


Single section part 


Single plate part 


Minor assembly 


Sub assembly 


Flat plate assembly 


Flat panel assembly 


Major sub assembly 


Curved panel assembly 


Sub unit assembly 


Unit assembly 


has a task for ship production. Tab. 3 shows the work stations 
which have function in double bottom block production. 

In edge cutting station (11), the edge cutting operation of 
ship hull plates is carried out. Edge cutting operation is the 
contour cutting of flat plates. And the plates which are subject 
to edge cutting constitute the panel structure. The edge cleaning 
operation of the ship hull plates, which are cut in edge cutting 
station (I1), is carried out in edge cleaning and sequencing 


station (12). There are some materials and slags on the edge 
surfaces of the plates after edge cutting. Using a grinding 
machine, these materials and slags remove from the edge 
surfaces of the plates. In I2 station, the plates are also sequenced 
in accordance with the process turn. The hull plates are welded 
and the panel structure is produced in panel production 
station (13). In panel cutting station (14), the panel which is 
manufactured in panel production station (13) is subject to 
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Tab. 2. Outfitting materials used in double bottom block 


Outfitting name Outfitting structure 


Manhole 


Bottom plug 


Vertical ladder 


Doubling plates 


Tab. 3. The work stations in the production process 
of the double bottom block 


Station no Station name 


Edge cutting 


Edge cleaning and sequencing 


Panel production 


Panel cutting 


Profile spot welding 
Profile TIG welding 
Section spot welding 


Section TIG welding 


Grinding 


Profile cutting 


Profile bending 


Nest cutting 


Pre-fabrication 


Jig 
Plate bending (Press) 
Block assembly 


counter cutting in accordance with its dimensions. The profiles 
are assembled on the panel by spot welding in profile spot 
welding station (15). The profiles are welded by TIG welding in 
profile TIG welding station (16). The minor and sub assemblies 
are joined on the flat panel assembly by spot welding in section 
spot welding station (I7). The minor and sub assemblies are 
welded on the flat panel assembly by TIG welding in section 
TIG welding station (18). Grinding station (19) is the last station 
of the panel line. In this station, the grinding operations of the 
flat panel and major sub assemblies are performed. The cutting 
operations of the profiles are performed in profile cutting 
station (110). Standard-dimensioned profiles, which are sent 
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to profile cutting station, are cut with specific dimensioned 
profiles. The bending operations of the profiles are performed 
by box machine in profile bending station (111). The bending 
profiles are used in curved panel. Nest cutting station (I12) is 
the heart of the shipyard production system. In this station, 
single plate assemblies are manufactured. Minor and sub 
assemblies are produced in pre-fabrication station (113). Curved 
panel assemblies are produced in jig station (114). Jig structure 
consists of curved jigs. The curved panels are lied down the 
jig structure and the curved profiles are welded on this curved 
plates. In plate bending station (115), the bending operations 
of the plates, coming from nest cutting station, are performed. 
Therefore, the flat plates are transformed to the curved plates. 
The structures and parts produced in previous work stations 
are sent to block assembly station (116) and the block structure 
is formed by assembling the corresponding parts. Fig.2 shows 
the work flows between work stations. 


Fig. 2. Work flow in the shipyard production system 


Detailed process analysis of current production 
case (Step 3) 


In this step, the comprehensive process analysis of the work 
stations are carried out. Every work station has a considerable 
number of work activities. At the beginning of the detailed 
process analysis, the work activities are determined and then 
the durations of each work activity are achieved. Because there 
are over 2500 work activities, the process analysis of only panel 
production and panel cutting stations are presented here. 


Process analysis of panel production station (13) 


In panel production station, panel structure is fabricated by 
submerged arc welding. Fig. 3 shows the general arrangement 
of the panel production station. In the first step, the plates 
are fixed to panel line by using conveyor fixing mechanism. 
Then, they are sent to tolerance plate’s welding area where 
the tolerance plates are welded with tack welding. Submerged 
arc welding operation starts from the tolerance plates because 
the welding becomes more stabilized. After that, the plates are 
transferred to submerged arc welding machine with conveyor. 


Tab. 4. Detailed process analysis of panel production station (I3) 


Repetition Activity 
Activity description number of duration 
activity (min.) 


The overhead crane goes to the flat plate to be transported 


The overhead crane comes down to the flat plate surface 


The overhead crane holds the flat plate 


The overhead crane lifts the flat plate 


The overhead crane transports the flat plate 


The overhead crane takes down the flat plate 


The overhead crane leaves the flat plate surface 
The overhead crane goes up the flat plate surface 


YSI uin ATR] wl um] =e 


The operator walks to the start point of the conveyor line 


jak 
© 


The worker walks to the conveyor line plate fixed mechanism. 


jà 
j 


The operator fixes the plate on the conveyor line 


= 
N 


The worker runs the conveyor line 


— 
w 


The operator walks to the welding tolerance plates 


j= 
& 


The conveyor transports the plate to the tolerance plate welding area 


= 
Nn 


The alignment of the plate on the tolerance plate welding area 


jæi 
a 


The operator brings the tolerance plates to the flat plate 


— 
N 


The operator puts the tolerance plates near the corners of the flat plates 


— 
oo 


The worker walks to spot welding machine 


— 
O 


The worker runs the spot welding machine 
The operator sets up the spot welding machine 


The operator fixes the tolerance plates on the corners of the flat plates with 
spot welding 


The operator walks to the other corner of the flat plate 


The operator goes to the TIG welding machine 


The conveyor transports the flat plate to the TIG welding area 


The worker walks to the TIG welding machine in order to check the position 
of the flat plate 


The worker checks out the position of the flat plate 


The operator runs the TIG welding machine 


The flat plate is fixed to the TIG welding machine with clamps 


The worker checks out whether the plate is fixed on the TIG welding 
machine conveniently or not. 


The worker brings the spot welding machine near the TIG welding machine 0.72 


The operator comes down to the plates for spot welding 0.53 


The spot welding operation of the plates prior to TIG welding 115.358 


The worker goes to the torches of the TIG welding machine 0.32 


The operator sets up the TIG welding machine before TIG welding process. 41.664 
TIG welding of the plates 272.19 
The conveyor transports the panel outside the TIG welding station. 1.728 


Here, the plates are welded with submerged arc welding and panel fabricated in I3 gets to Buffer Area 2. Then, the panel 
they are then sent to buffer area. As a result, flat plate assembly is transferred to panel cutting machine with conveyor. Panel 
(production stage E) is fabricated. Table 4 illustrates the detailed cutting machine performs counter cutting operation. But, 


process analysis in panel production station. before this, blasting operation is fulfilled. Then, automatic 
marking operation is done. Therefore, the alignments of the 
Process analysis of panel cutting station (14) piece parts can be easily performed. After that, the counter 


cutting operation is carried out. Finally, the panel is transferred 
Counter cutting of the panel is performed in this station. to Buffer Area 1. Detailed process analysis of panel cutting 
Fig. 4 shows general arrangement of panel cutting station. The station is illustrated in Table 5. 
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ESS" Submerged Arc Welding =: 
psa Machine i: 


ef ae Tolerance plate’s welding 


; [its Conveyor fixing mechanism =: 
Fig. 3. General arrangement of panel production station (I3) Fig. 4. General arrangement of panel cutting station (I4) 
Tab. 5. The process analysis of panel cutting station (14) 
Repetition Activity 


Activity description number of duration 
activity (min.) 


N 


The operator transports the panel to the panel cutting station 


The operator goes to bring the cutting program CD 


The operator takes the cutting program CD 


The operator goes back the cutting station with CD 


The operator loads the Program CD to the cutting machine 


The cutting machine moves on the counter line of the panel 


The cutting machine signs the reference points on the panel 


The cutting machine loads the reference points on itself. 


0.5 
139.937 
101.457 

2 
17.583 


The operator runs the panel cutting machine 


Blasting and marking processes 


Cutting process 


The worker takes the marking pen 


NÍIÍINIÍININIÍINIÍNINIÍNJINININIJN 


Manuel marking process 


Taking the outfitting parts and bringing them to the cutting 
station. 


— 
N 


The workers put the outfitting parts on the panel. 


The worker runs the spot welding machine 


Spot welding process 


The worker runs the grinding machine 


Grinding process after the spot welding process 


The worker runs the TIG welding machine 


TIG welding process 
Grinding process after the TIG welding process 


The conveyor transports the panel outside the panel cutting 
station 
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In the same way, the detailed process analysis are performed for the other work stations. Because there are over 2500 work 
activities, it is impossible to present here. After the work activities and their durations are determined, the completion times of 
the work stations are calculated by considering the serial and parallel work activities. In Table 6, the completion times of the 
work stations are given. For instance, the process time of edge cutting station is 190 minutes. In other words, the edge cutting 
process of the plates, which constitute the panel of the double bottom block, takes 190 minutes. 


Tab. 6. Station completion times of the current case 


Station name Completion time (min.) 


Edge cutting 190 
Edge cleaning and sequencing 203 


Panel production 622 
Panel cutting 


Profile spot welding 
Profile TIG welding 
Section spot welding 


Section TIG welding 


Grinding 


Profile cutting 


Profile bending 


Nest cutting 


Pre-fabrication1 


Pre-fabrication2 


Jig 
Plate bending (Press) 


Block assembly 


Simulation model of current production case (Step 4) 


In this step, simulation model of the current production system is created. Fig. 5 depicts the simulation model of the current 
case and Tab. 7 shows the modules which are used in simulation model. 
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Fig. 5. Simulation model of the current case 
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Module name 


Createl 


Tab. 7. Definition of module numbers 


Module name 


Assign13 


Module name 


Process of profile bending 
(111 station) 


Process of plate edge cutting 
(I1 station) 


Assign14 


Assign 27 


Process of plate edge cleaning 


Seperate 1 


Assign15 


Batch3 


Seperate 5 


Assign 28 


Assign | 


Process of section spot 
welding(I7 station) 


Assign 29 


Assign2 


Process of section TIG 
welding (I8 station) 


Batch 5 


Assign3 


Assign16 


Process of pre-fabrication 1 


Assign4 


Assign17 


Assign 30 


Delay1 


Process of grinding 
(19 station) 


Separate 6 


Process of panel production 
(13 station) 


Decidel 


Batch4 


Assign18 


Assign 31 


Assign 32 


Assign5 


Assign6 


Create2 


Process of profile cutting 
(111 station) 


Batch 6 


Process of pre-fabrication 1 


Delay2 


Assign19 


Assign 33 


Process of panel cutting 
(14 station) 


Separate 2 


Process of plate bending 
(115 station) 


Decide2 


Assign20 


Separate 7 


Assign7 


Decide4 


Assign 34 


Batch1 


Separate 3 


Batch 7 


Assign8 


Assign 21 


Process of jig (114 station) 


Delay3 


Assign 22 


Assign 35 


Assign9 


Create 3 


Assign 36 


Assign 10 


Process of nest cutting 
(112 station) 


Assign 37 


Batch2 


Assign 23 


Assign 38 


Assign! 1 


Separate 4 


Batch 8 


Assign 12 


Assign 24 


Process of block assembling 
(116 station) 


Process of profile spot 
welding (15 station) 


Decide 5 


Dispose 


Process of profile TIG 
welding (16 station) 


Assign 25 


Decide3 


Assign 26 


In ARENA simulation model, various modules are employed to create the simulation model of the production system. Tab. 8 
represents the modules and their definitions. 
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Create 


Tab. 8. Description of modules 


Description 


Describes the arrival time and quantity of the 
materials entering into the production system 


Process 


Describes the work stations in the production 
system 


Assign 


Makes assignments the products leaving any 
work station and it is also used in assembly 
operations 


Decide 


This module orients the products to the work 


Determination of i Ht det of current 
production case (Step 5) 


After the simulation model is created, the model is run 
along 720 hours. The current case has a capacity of producing 
18 double bottom blocks at the end of the operation time of 720 
hours. In other words, the current case can produce a double 
bottom block in 2400 minutes. 


Suggestions on the current case (Step 6) 
In the current production line, the assemblies of the outfitting 


and piping are carried out in the block assembly station. After 
the completion of the steel works, the outfitting assembly and 


stations where they will go to 


then the piping assembly are carried out. Therefore, all these 


Batch 


Used in assembly operations. It combines the 
products for assembly operation 


operations are done in serial way, which means more work 
load and longer lead times. If we can carry out the outfitting 
and piping operations in earlier work stations and distribute the 


Seperate 


Divides the main product into sub products. It is 
usually used for cutting operations 


work load, we can reduce the work load and the completion 
time of the block assembly station. By doing this, the waiting 


Dispose 


No. of 
suggestions 


Shows the exit of the production system 


Workstation 


Panel cutting 


time in front of the block assembly area is lowered and the lead 
time of the double bottom block may be shorter. Tab. 9 depicts 
the suggestions on the current case. 


Tab. 9. The suggestions on the current case 
Description of suggestions 
The welding operations of the manholes on tanktop panel were carried out. These welding 


operations are one side welding. In this way, the one side welding operations of the 
manholes will be finished 


Panel cutting 


The one side welding operations of the bottom plugs (two pieces) on hull panel 


Panel cutting 


The welding operations of the zincs (six pieces) 
on hull panel were carried out and finished 


Pre-fabrication 


The assemblies of vertical ladders (three pieces) were carried out on the sections 


Pre-fabrication 


The assemblies of the zincs (thirty-two pieces) were carried out 


Pre-fabrication 


Assemblies of pipe systems (twenty-six) on the module 


Pre-fabrication 


Manufacturing the module in pre-fabrication station (113) 
in order to assemble the piping systems 


Jig 


Assembly of pipe systems (nine) on the curved panel assembly in jig station 


Block assembly 


Assembly of one vertical ladder when the block is upside down 


Block assembly 


Completing of the welding operations of two bottom plugs when the block is upside down 


Block assembly 


Completing of the welding operations of four manholes when the block is in flat position 


Block assembly 


Assembly of twenty doubling plates when the block is in flat position 


Block assembly 


Assembly of two bilge wells when the block is in flat position 


Block assembly 


Assembly of four zincs when the block is in flat position 


Block assembly 


Carrying out the outfitting and steel operations 
in parallel way in the block assembly station 


Block assembly 


In the block assembly station (116), in the steel operations, there are fourteen TIG welding 
and spot welding workers in the current case. In the new case, the quantity of the TIG 
welding and spot welding workers are sixteen for steel operations. The TIG welding and 
spot welding operations are carried out by the same workers. And also, two spot welding 
and TIG welding machines which are on idle are added to the station 
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Application of the as oe on the current 
case (Step 7) 


In this phase, the suggestions given in Step 6 are applied in 
the current simulation model in Fig. 5. Panel cutting, section 
spot welding, pre-fabrication, jig and block assembly stations 
were influenced by the suggestions on the current case and the 
completion times of these work stations are calculated again. 
Table 10 shows the completion times of each work stations in 
new case. 


Tab. 10. Station completion times of the new case 


Completion time 


Station name 7 
(min.) 


Edge cutting 


Edge cleaning and sequencing 


Panel production 


Panel cutting 


Profile spot welding 
Profile TIG welding 
Section spot welding 


Section TIG welding 


Grinding 


Profile cutting 


Profile bending 


Nest cutting 


Pre-fabrication1 
Pre-fabrication2 
Jig 
Plate bending (Press) 
Block assembly 


Evaluation of the new case (Step 8) 


In the evaluation of the new case, Arena simulation program 
will be used. The new case manufactures 24 double bottom 
blocks at the end of 720 hours. That means the new case can 
produce a double bottom block in 1800 minutes. 


Comparison of the current and new cases 
(Step 9) 


Fig. 6 shows the quantity of the double bottom block 
manufactured in both cases at the end of 720 hours. While the 
current case produces 18 double bottom blocks, the new case 
produces 24 double bottom blocks. 


30 


BLOCK QUANTITY (UNITS) 
rR 
Oo 


CURRENT CASE NEW CASE 
PRODUCTION CASES 


Fig. 6. Block quantity manufactured in the current and new cases 
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Fig. 7 shows the cycle times of the double bottom block for 
both cases. While a double bottom block is produced in 2400 
minutes in the current case, in the new case, a double bottom 
block is produced in 1800 minutes. 


3000 
2500 +4 
2000 + 
1500 +4 
1000 
500 


CYCLE TIME (MINUTES) 


CURRENT CASE NEW CASE 


PRODUCTION CASES 
Fig. 7. Cycle times of the double bottom block for the current and new cases 


It should be noted here that, no extra expenditures are paid 
during the transition from the current case to the new case. All 
the changes on the current case were performed by using the 
possibilities of the shipyards without purchasing any workers or 
equipments. Table 11 presents the need of additional equipment 
and manpower. 


Tab. 11. Need of resources after changing 


Suggestion | Work Station | Additional 
no. equipment 


Additional 
manpower 


Panel cutting No need No need 


Panel cutting No need No need 


Panel cutting No need No need 


Pre-fabrication No need No need 


Pre-fabrication No need No need 


Pre-fabrication No need No need 


Pre-fabrication No need No need 


Jig No need 


No need 


Block assembly | No need No need 


Block assembly | No need No need 


= | = 
ZiSleslælxȚxal oala] ajen 


Block assembly | No need No need 


Block assembly | No need No need 


Block assembly | No need No need 


Block assembly | No need No need 


Block assembly | No need No need 


No need No need 


Block assembly 
CONCLUSIONS 


In this study, some improvement suggestions with regard 
to outfitting and piping assembly operations were carried out 
on the shipyard production system. The improvements were 
applied on the block assembly station that was a bottleneck. 
The main changes on the current case were the pre-outfitting 
assembly operations and pre-piping assembly operations. 
According to simulation work, in the current case, the system 
can manufacture 18 blocks at the end of 720 hours. After the 
changes were applied on the current case, the new production 
case can produce 24 blocks at the end of 720 hours. So, the 
changes on the current case provide an improvement rate of 
33% without additional expenditure. In other words, when the 
assembly operations of the outfitting and piping are carried 


out in earlier work stations, the throughput increases. The 
pre-outfitting and pre-piping assembly operations provide less 
cycle time. The authors strongly recommend the shipyards to 
perform outfitting and piping assembly operations in earlier 
stations instead of block assembly area since the cycle time of 
the block is considerably decreased. 
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On the ship course-Keeping control system 
design by using robust feedback linearization 
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ABSTRACT 


In the paper the problem of ship autopilot design based on feedback linearization method 
combined with the robust control approach, is considered. At first the nonlinear ship model 
(of Norrbin type) is linearized with the use of the simple system nonlinearity cancellation. 
Next, bearing in mind that exact values of the model parameters are not known, the ensuing 
inaccuracies are taken as disturbances acting on the system. Thereby is obtained a linear 
system with an extra term representing the uncertainty which can be treated by using robust, 
HĪ optimal control techniques. The performed simulations of ship course-changing process 
confirmed a high performance of the proposed controller despite the assumed significant 


errors of its parameters. 


Key words: feedback linearization; robust control; differential games; ship autopilot design 


INTRODUCTION 


The feedback linearization (FL) method [9, 10, 11, 13] 
consists in such transformation of a given nonlinear system 
which results in a new, linear time-invariant one. Here, by the 
transformation we mean the application of a proper control law 
combined with possible change of the system coordinates. Once 
a linear system is obtained, a secondary control law (or sub- 
control) should be designed to ensure that the overall closed- 
loop system performs according to the specifications. 

In the simplest case of the system the FL method is reduced 
to the ordinary cancellation of nonlinearity by means of 
a properly selected control function. 

One of the main drawbacks of FL method relates 
to inaccuracies arising during cancellation of system 
nonlinearities. The thereby obtained transformed system is in 
fact not perfectly linear and, moreover, these imperfections may 
often prevent the use of efficient techniques of linear systems 
synthesis. An effective way to solve this problem consists in 
combination of feedback linearization method with the robust 
control techniques. In this paper H. optimal control theory 
within the state space framework is applied, i.e. the problem is 
considered from position of differential games theory [4, 5, 8, 
16]. In this view the model uncertainties are considered as an 
action of adversary player (or opposing nature) while our part 
is to invent a control strategy that is the best in terms of some 
given quality criterion (cost functional). In other words we are 
trying to minimize the cost assuming the ‘worst-case action’ of 
our opponent player (disturbances). Such an approach allows 
devising a controller which taking into consideration system 
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parametric uncertainties, guarantees at the same time a good 
process performance. 

In the paper, besides presenting a relevant portion of the 
above stated theory, its usefulness to ship course-keeping (or 
changing) problem is considered. 

It is known that the PID (proportional-integral-derivative) 
controllers traditionally used in the field still have many 
shortcomings. 

A disadvantage of a PID controller is that it can provide 
optimal performance only at the operating point for which 
it is designed. The ship parameters vary significantly with 
operating conditions such as e.g. forward speed of the vessel. 
Under the varying operating conditions, is tedious and difficult 
to determine properly the fixed parameters of the controller that 
results in good performance. Furthermore, the PID autopilots 
can cause difficulties when the ship makes large maneuvers 
involving non-linear dynamic behaviour. To avoid the problems 
of fixed structure PID autopilots, adaptive autopilots were 
introduced in the 1970s and have remained a major area of 
research until now [2, 17]. 

An alternative and promising research direction is the 
autopilot design from position of robust control theory for 
nonlinear systems. Therefore the main motivation of this paper 
is a proposal of the robust ship autopilot design which, based on 
recent H,, theory, is able to finely cope with above mentioned 
inconveniences. 

The paper is divided into five sections and ended with 
conclusions. The second section presents the system class 
definitions, basic concepts as well as the transformation of 
considered system into standard differential game form. In 


the third section a course-keeping problem and the Norrbin’s 
ship model structure are given. In the fourth section the two 
cases of robust ship autopilot design are presented (with and 
without load disturbances taken into account). The last, fifth 
section includes a short description of simulation tests and 
their results. 


BASIC CONCEPTS AND DEFINITIONS 


Let us consider a nonlinear system in Brunovsky form 


[11]: 


X, =X, 
X% =X; 
(la) 
X, = f(x, 0f) + g(x, 0s)u 
y=xX, (1b) 
where: 
xe R — the state vector, 


Ofe Rk, Oee R! — vectors of system parameters, 
ue R — the control input, 
yeR — the system output. 


In the case of exact system model, i.e. under assumption 
that the nonlinear functions fand g of the model (1a) are known 
(and g(x, 02) #0, Vx € R») the insertion of the simple output 
feedback linearizing controller [9, 10, 13, 18]: 


TE f(x,0 +v 
g(x, 0°) 


in the system (la) results in exact cancellation of both 
nonlinearities (f and g), which yields: 


yo=y 6) 


(2) 


To find the control v(t) stabilizing this linear system 
a standard pole location technique can be used. If v is chosen 
as: 


vey uy (4) 


and the coefficients p; are such that [(s) = st + ps) +... u; is 
a Hurwitz polynomial in the Laplace variable s, then the output 
y(t) and its derivatives converge to zero asymptotically, because 
the closed-loop dynamics is reduced to the equation: 


yO + wyed +... by =0 (5) 


which, by virtue of the choice of coefficients u; is asymptotically 
stable. Let’s notice that also the internal stability i.e. the state 
x — 0 as t > œ% is obtained. 

As the true system parameters 0, 02 are unknown and only 
some their estimates f, 6% are at our disposal the nonlinear 
functions of the system (la) are f(x, 0°) + g(x, 6®)u and the 
control law (2) takes now the form: 


pert (x64 +v 

g(x, 6°) 
Since the insertion of the control (6) into the system (1a) 
no longer guarantees exact cancellation and thereby a resulting 
system linearity (like in the former case (2)) we will try to solve 
this problem from position of differential games. To transform 


our system into a suitable form we perform the following 
calculations. 


(6) 


Adding to and subtracting from the last equation of the 
system (la) the sub-control v we get: 


x,=f-f+(g-g)ut+v (7) 
m 2 WwW 
where the formula v = f + gu obtained from (6) was also used 
here. 
Now denoting w = f— f + (g — g)u and treating w as system 


disturbances we may rewrite the equations (1a) into the general 
matrix form [3]. 


x =Ax+Bv+B,w (8a) 
y= Cx + D,,u (8b) 
01 0 - 0 0| 
0 0 1 0 0 
where = , B,=B,=}: (9) 
0 1 0 
0 0 0 0 1 


To _ T = 
DiC, =0. DD, =I 
Let’s observe that for optimization purposes a new output, 
(8b), is defined. 


The matrices C, and D,, should be so selected as to achieve 
the proper weighting in the following cost functional: 
2 Z y2 = 
Jv.) =y -v O= 
"i (10) 
= (FOYO -y W Od 
0 
where y > 0 is a properly chosen constant (see below) called 
the performance bound. 
The last two conditions of (9) are assumed to avoid cross 


terms in the functional (10) so that to get an equivalence to the 
corresponding Lar criterion, i.e.: 
T T 
Q=C'C,, R=D'D,=1 
(compare (33), (34)). 

We have defined a differential game of two players where 
each of them is trying to influence the system by proper (profitable 
for him) choice of his strategy — v and w, respectively. 

We (as the first player) are trying to minimize (10) with 
respect to v while the disturbances (the second player) is 
assumed to maximize the cost (10). 

Assuming the commutativity of min-max operators as 
well as that the optimal strategies v*, w* of both parties exist 
we have: 


(11) 


min max J,(&,V, w) = J, v; w) 


(12) 


Let’s observe that minimization of the maximized, by 
the second player, cost (10) refers to the situation where we 
are doing our best (optimally stabilizing the system output) 
assuming the worst-case realization of system disturbances. 

One can prove [3] that if for a given y > 0 w* maximizing 
the cost (10) exists, this fact can be interpreted in terms of the 
condition: 


Io. <1 03 
where: 


K — applied controller which denotes boundedness of the 
infinity norm of the closed-loop transfer function G,,, from the 
disturbances w(t) to the output y(t). This guarantees also the 
system L,-gain stability (from w to y) [1]. 
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Minimization of (10) with respect to v, under assumption 
that the minimal bound y of is obtained, allows in turn to find 
the optimal controller K*. 

The explanations can be summarized into the following 
condition: 


min [Cn] = [GK 9-1 09 
where: 
y* — the smallest y in the condition 
K — aset of available controllers. 
The formulas for suboptimal strategies are [3, 1]: 
v(t)= — By P(t) x(t) 
(15) 


w(t) =y" By, P(t) x(t) 


where matrix P(t) is a solution of the Hoo-like algebraic Riccati 
equation (RE), 


PA+A'P-P(B,B'—y°B,B!)P+ 


h . (16) 
+C/C,=0 


such that: 
A-(B,B! —y°B,B.)P (17) 


is stable, i.e. all of the eigenvalues of this matrix have negative 
real parts. 

In practice, to solve the above formulated problem, we have 
to use some iterative procedure which consists in repeating the 
required calculations for different values of y in order to choose 
the smallest one. For the properly selected y (i.e. for y = y*) the 
strategies (15) are optimal. 


COURSE-KEEPING PROBLEM AND THE 
NORRBIN’S SHIP MODEL STRUCTURE 


The problem of course-keeping is a task of designing 
an automatic control aid (autopilot) which is able, by using 
appropriate rudder actions, to control the ship motion as to 
maintain a pre-assigned constant heading. This problem is 
often considered as consisting of two sub-problems. The first 
one concerns the ship controlling along the straight line at small 
course deviations. The second problem concerns the proper ship 
control during large manoeuvres, i.e. the problem of quality of 
transition process as a reply for a course step-change. 

In order to synthesize a course-keeping controller we apply 
the following Norrbin’s [7, 12] ship model general structure: 


Tý + FQ) = kô (18) 
where: 
w — actual ship course angle (heading) 
ò — rudder deflection angle as a control variable 
T, k — unknown model parameters 
F() — unknown function with assumed structure 


In the ‘classical’ approach to ship control the structure 
of the nonlinear function F is (according to Norrbin model) 
often assumed in the form of a polynomial of the third order. 
Generally it may be assumed as follows: 


F @) = aa}? +a? + ay + ao (19) 

For the ship with hull symmetry we have a, = 0. The bias 

term a, is frequently taken as null, being conveniently treated 

as an additional rudder-offset which can be made null by 

the adequate selection of the integral action in the autopilot 
design. 
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The most commonly used structure has therefore the 
form: 


F Qj) = aay +ayy (20) 


Now, assuming that the structure of the function F has been 
predetermined the coefficients a; are usually identified during 
sea trials [12]. After some ship-circulation tests we apply 
regression analysis to the obtained data (y in function of ô), 
for each of the mentioned structures separately, and then we 
search for the best fitting solution. 

Since each of the tests should be performed for prescribed 
sailing conditions, e.g. different ship load, trim and velocity 
it follows that the problem of sufficiently general model 
building is a laborious and expensive task. For this reason in 
practice linear models are preferred [14] as being simpler for 
identification as well as we have, in this case, a number of 
linear synthesis methods at our disposal. For several control 
tasks however, especially for strongly nonlinear objects, 
the linear models are insufficient. The control algorithm 
obtained from such model leads not only to the deterioration 
of control performance but may also produce an unpredictable 
destabilization of system. Due to the above mentioned facts 
we propose an approach which, while dealing with nonlinear 
models, avoids, at the same time, to cope with the demanding 
identification task. 


ROBUST SHIP AUTOPILOT DESIGN 


In this section, based on the above given theoretical facts 
as well as on the presented Norrbin’s model structure, we 
formulate and then solve the problem of robust ship course- 
keeping control synthesis. It is assumed, for simplicity (but 
without loss of generality), that the desired (reference) course 
value equals to zero: y, = 0. 


The case without additional load disturbances 


To apply the above mentioned theory let us first define the 
course—keeping problem as a differential game. To this end we 
rewrite the Norrbin’s model (18) in the state space form: 


=r 
j (21) 
r= D(r)+cd 
where: 
= —F()/T, 
c= kT 
r — course turning rate (angular velocity). 
According to (20) the function ® is: 

Diy) = baj’ + bay (22) 
where: 
b; =a/T, i= 1,3 


To make use of simple feedback linearization method we 
apply (compare (6)) the controller: 


(23) 


ae ey - (bay? + bay) +v 


A 


Cc c 


Now, following the procedure given in section 2, we can 
write Eq. (21) in the form: 


y =r 
t= v+ {b, -bj +b -bý +(c- 98} 
S y 


w 


(24) 


or simply as follows: 


oH th 


where w denotes errors caused by the model parameter 
uncertainty (considered here as a control strategy of the 
opposing nature). 

The system output may be assumed as follows: 


rb lookebte] 


where the constant A is a criterion weight. 
The matrices which define the system in question are as 
follows: 


01 0 a 0 0 
A= B,=B,=|.| C =|" | D =| | en 
0 Oo}, i). (oa) p 


Let’s note that by the proper choice of the output matrices 
C,, D,, a desired cost criterion (compare (11)) is obtained. 

To find the proper numerical solution P(t) of Riccati-like 
equation (16) we can use Matlab function care [6] for different 
values of the constant y > 0. As we want to select y to be 
minimal (provided the solution RE exists) we have to apply 
this procedure repeatedly. This way we get the proper control 
strategy (see (15)) which will be used later in the process 
simulation. 


The case with added load disturbances 


During the routine work of the autopilot (i.e. under specified 
steady-state conditions, e.g. straight route mode) besides 
modelling errors we have also to take into account the load 
disturbances. The disturbances lead to the steady heading 
error which can be typically compensated by integral action 
introduced to the control system. 

To cope with the steady-state error we will consider in this 
section the problem of robust regulator design equipped with 
integral action. 

Let us first define an extra error in the form of integral of 
the course error: 


e= fva 


0 
Based on the model (25) we define the augmented system by 
adding the error equation (28) (written in differential form): 


wer wl {0 1 Olly} J0 0 0 
t=v+w+d |rl=/0 0 Of rj+{ljv+{l}w+]1 jd (29) 
e=y lè] [1 0 oje] [0 0 0 


(28) 


where d are constant or slowly varying load disturbances. 
The complementary output equation is: 


yı A O0O0ly] 10 Aw 
y=ly,;=|000//rj+{liv=| v 
y.) |OO01};e] 0 e 
Thus the matrices that define our system are as follows: 
010 0 200 0 
A=|0 0 0|,B,=B,,= 1C, = 00 0),D,,= 1| G1) 
10-0 0 001 0 


(30) 


Let’s notice that during controller synthesis process the 
disturbances d are not directly taken into account. However, 
thanks to above presented system structure, any appearance 
of disturbances during system operation leads to controller 
counteraction which prevents steady-state error. 


COURSE-CHANGING PROCESS 
SIMULATIONS 


The standard method of assessing the control system quality 
is based on analysis of the transition process as a response to 
the step input. Thus in the following simulations we will test the 
ship behaviour after step change of the course set-point under 
different controller structures as well as various criteria data. 


Ship simulation model 


According to [15] we assume, as the actual parameters of 
the ship model (21), the dynamic manoeuvring parameters of 
the m/s Compass Island. The units of time, length and angle 
are one minute, one nautical mile and one radian, respectively. 
The parameters were determined as follows: 


b, =— 0.62 min, b, =— 1.064 min“, 
b, = b, = 0, c = 3.553 rad/min 


The parameters are, of course, not known to the control 
system designer thereby during the process simulations we 
will take another set of their values — possibly, substantially 
different. 

It is assumed also that the ship’s velocity during the steady- 
state control process phase is constant u = 0.25 nm/min and 
that all the time we have constant propelling force. 

The ship has the following characteristics: the maximum 
rudder angle is 35 deg., gross register tonnage of 9214 RT, 
deadweight of 13498 t, length of 172 m., draught of 9.14 m, 
one propeller, and maximum speed of 20 kn. Let’s note that the 
assumed parameters make the ship directionally stable [12, 14] 
but we can get similar results for the ship which is directionally 
unstable (e.g. for the opposite sign of the coefficient b,, i.e. 
b, =— 1.064 min). 


Simulation results 


(32) 


In order to examine the performance of the robust controller 
in question a corresponding Lqr regulator [3] has also been 
tested for comparison. The matching criteria matrices for the 
two cases were chosen the same. 

As an approximate counterpart values to the actual ship 
parameters (32) (unknown for designer) we assume: 


b, =— 1 min, b, =— 2 min", 


a 


b, = b, = 0, c= 2 rad/min 


The parameters are essential for the controller design as they 
contain some uncertain data of the ship’s dynamics. 
The remaining numerical data are as follows: 


y* =1.05,4=2 


and 


(33) 


where the criteria matrices Q and R refer to the corresponding 
Lar problem. Let’s notice that the matrices are exactly the same 
as those which occur in the matching robust case (compare 


(26)). 
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Fig. 1. Ship headings versus time for: a) robust controller, 
b) Lar regulator (30° turning) 


The following figures show the situation where the ship 
is moving for 3 min with steady speed (0.25 nm/min) along 
a straight line and then we apply a 30-degree step heading 
change command (in the case given in Fig. 3 and 4a disturbing 
torque is additionally applied). 

The graphs presented in Fig. 1 and 2 show respectively 
the plots of headings and rudder deflections versus time for: 
a) robust controller, b) Lqr regulator. 

In Fig. 3 and 4 are presented corresponding plots but for 
the case where the controllers are provided with integral 
action. 

Let’s notice that despite the fact that the order of parameter 
errors is here about 100% the robust controller performance 
(Fig. 1a) is very good. On the other hand the corresponding Lqr 
controller operates close to the stability margin (see Fig. 1b). 

Moreover, according to the author experience, the plant- 
model parameter errors for the robust controller can be made 
substantially higher without noticeable process performance 
degradation. 

Fig. 3 and 4 show corresponding graphs, however, with 
the difference that the controllers: robust and Lqr one are now 
equipped with integral action. The intensity of the constant 
load disturbances d acting on the ship, in the form of torque, 
corresponds to the compensating deflection of the rudder by 
11 deg (see Fig. 4). 

The numerical data for this case are as follows: 


y°=1.2,4=7 


74 POLISH MARITIME RESEARCH, No 1/2013 


Deflection [deg] 


“0 2 4 6 8 10 12 #14 16 18 20 
Time [min] 


Deflection [deg] 


U52 4 6 s l0 12 l4 l6 I$ 20 
Time [min] 
Fig. 2. Rudder deflections for: a) robust controller, 
b) Lqr regulator (30° turning) 
and 
2 

x00 
Q=|0 0 0,R=1 (34) 

001 


where the criteria matrices Q and R refer, as before, to the 
corresponding Lar problem. 


As it can be seen, also in this case the robust controller 
performance is still better than its Lqr counterpart. 

Let’s observe that the plant-model mismatch disturbances 
w are diminishing with the progress of a transient phase of the 
process regulation (compare (24)). Thus the proper effect of 
integral action is manifesting itself during the system steady- 
state, i.e. just after the process transient phase. 

By manipulating the criteria weighting à we can adjust the 
required trade-off between precision of the output variable 
control (y) in relation to the control effort (v) and/or the integral 
of course error (e). 

Taking e.g. 4 = 7 we highly penalize the course error in 
relation to the control effort v as well as to the integral error e. 
To increase the number of degrees of freedom in this respect 
we could introduce some additional weighting related to the 
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Fig. 3. Ship headings versus time for controllers with integral action: 
a) robust controller, b) Lqr regulator (30° turning) 


state variable e. High penalization of this variable results in 
a rapid decrease of the course steady-state error, however, at 
the cost of increasing the course overshot. 


CONCLUSIONS 


The paper describes a proposition of a novel technique for 
robust control system design, which can be used to the ship 
motion control. 

The technique seems to be very attractive for designing the 
ship autopilot albeit it requires additional research which should 
take into account different types of ship nonlinear models (e.g. 
including a steering gear) as well as consider various sailing 
conditions (in this phase of research the disturbances from 
waves has been deliberately omitted). It is worth noting that 
the modern adaptive PID controllers only partially solve the 
problem of ship model parameter changes. They can tune (as 
a rule) their gains only in response to the ship velocity changes 
while the number of factors responsible for the ship dynamic 
variations is much greater (e.g. ship load). The factors, when not 
taken into account by the controller, may lead to the degradation 
of steering process performance. Therefore, each new design 
proposal in this respect should attract attention. 

The presented paper is focused on the properly formulated 
theory rather than on the practical engineering issues and is 
devised as a theoretical introduction to more practice-oriented 
engineering research. 
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ABSTRACT 


Over the past few years considerable advances in sonar technology, spatial positioning capabilities and 
computer processing power have lead to significant improvements in mapping, imaging and technologies 
of seafloor exploration. Recently, modern multibeam echosounder systems (MBES) capable of recording 
backscatter data for the whole water column, not just for the seabed, have become available thus providing 
data allowing for visualization and analysis of objects other than the seabed such as single fish, fish schools 
or pollution. Unlike bathymetric sonars, which only capture the seafloor, multibeam systems produce very 
large amounts of data during surveys. Because of this, storing the data collected during hydrographic or 
scientific cruises becomes a crucial problem. In this context, the paper proposes a new approach for efficient 
reduction and storage of MBES records. The results of a sample implementation of the algorithm being tested 
on several different sets of MBES data are also discussed. 
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INTRODUCTION 


Modern multibeam echosounder systems (MBES) are 
capable of recording backscatter data for the whole water 
column, which allows for analysis and visualization of not only 
the seabed, but objects such as fish schools and single fish as 
well. They produce very large amounts of data during surveys, 
as it has been shown eg. by Buelens et al. (2005). Current 
trends in sonar development involve the use of innovative 
transducer materials as well as application of sophisticated 
processing techniques including focusing algorithms that 
dynamically compensate for the curvature of the wavefront 
in the nearfield and thus allow narrower beam widths (higher 
lateral resolution) at close ranges. Future developments will 
probably focus on “hybrid”, phase-comparison/beam-forming 
sonars, the development of broad-band “chirp” multibeam 
sonars, and perhaps synthetic aperture multibeam sonars 
(Demkowicz & Bikonis, 2006). Not long ago, Buelens et 
al. (2006) considered storage and compression as one of 
the most important computational challenges in processing 
of MBES data. However, since Ferguson’s and Chayes’s 
(1995) proposal of a binary file format for multibeam sonar 
data storage, there has been little development in the field of 
efficient MBES data processing which would allow for fast, 
semi-real-time sharing of the results between diverse groups 
of interest such as fishermen, hydrographers or researchers 


(Chybicki 2008). Apart from works of Wu and Zielinski 
(1997), not much attention has been given to research of 
algorithms for storing and archiving MBES data that would 
allow for efficient browsing, analysis and visualization of 
collected information. 

Even the prospect of having a common file format in which 
the same data sets can be shared between researchers and 
hydrographers offers an interesting new perspective, leading 
to savings in instrumentation and survey costs. Efficient MBES 
data storage algorithms can also make Remotely Operated 
Vehicles (ROV’s) and particularly Autonomous Underwater 
Vehicles (AUV’s) to become more widely adopted as platforms 
for seafloor mapping systems. 

This paper presents the concept of a MBES data reduction 
and storage algorithm, developed by the authors for RESON 
Inc., along with the results of its sample implementation. 


MATERIALS AND METHODS 


MBES data, just like every other type of information, can be 
encoded using well-known methods. However, given the nature 
of MBES, a specialized approach could yield substantially better 
results. During the presented research, the authors designed and 
utilized a modification of Huffman coding (1952) for reduction 
of MBES record sizes according to the rules of information 
theory (1948), as described in the following section. 
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Entropy coding 


Entropy coding is a lossless data compression method based 
on information probability schemes. This area of knowledge 
was first mathematically described by Claude E. Shannon. In 
1948 he defined the basics of information theory in which he 
proposed concepts of measuring the amount of information. He 
assumed that a set of possible events S = {s,,s,, ...,8,} is called 
the alphabet if every element x; is a symbol used to construct the 
communique. He also assumed that the probability of every event 
is known as p(x;) = p;. The set of probabilities is described as 


P= {p,, P», ---P,} assuming that for every I(p;) > 0 and > p;= l. 


In this case, the function: I(x,) = — log,(x;) for k > 0 and k #1 
is called the self-information associated with event x;. 

According to Shannon (1948), self-information is the 
function for measuring the information brought by the particular 
event x;, which is reversely proportional to the probability of 
the event. 

The self-information thesis was extended by the definition 
of entropy which says that the function defined as: 


H(p,, Po, <--> Pa) 


is equal to: 


H(S) = -È pix log, (p,) (1) 


with the set of n independent events S = {s,, S» ..., S} and the 
set of corresponding probabilities P = {p,, P> ..., Pay- Entropy 
is the mean auto-information of the set of particular events 
and, when assuming a binary source, can also be considered 
as the smallest average number of bits required for a symbol 
to be encoded, as shown by Balakrishnan et al. (2007) and 
Titchener (2000). 


Huffman coding 


The first optimized method of compression coding based on 
the information theory was given by David Huffman (1952). 
The term refers to the use of a variable-length code table for 
encoding a source symbol (such as a character in a file) where 
the variable-length code table has been derived in a particular 
way based on the estimated probability of occurrence for each 
possible value of the source symbol. 

Huffman coding uses a specific method for assigning 
the n-bit representations for each entry symbol, known as 
a prefix code. The method consists of replacing commonly 
used symbols with shorter representations and less commonly 
used symbols with longer representations. To achieve this 
goal Huffman coders create a prefix-free binary tree of non- 
overlapping bit-sequences, where the length of each sequence 
is inversely proportional to the likelihood of the symbol to be 
encoded. 

Huffman binary code is biunique, meaning that the 
code is uniquely decodable and the sum of the probabilities 
associated with all symbols is always less than or equal to one. 
If the sum of probabilities is strictly equal to one, the code is 
termed a complete code. In general, a Huffman code need not 
be unique, but it is always one of the codes minimizing the 
complete code length. 

Huffman coding is a particularly good compression method 
for datasets with non-linear distribution, and can be applied 
to many different types of data such as voice, sound, text and 
others. The Huffman algorithm is currently one of the most 
popular and efficient compression methods, and it can work 
in linear time if input probabilities (also known as weights) are 
sorted. It has been shown (eg. by Larmore, 1986) that when the 
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actual symbol frequencies agree with those used to create the 
code, no other mapping of individual source symbols to unique 
strings of bits will produce a smaller average output size. 

Although Huffman coding is optimal for a symbol-by- 
symbol coding with a known input probability distribution, 
its efficiency has since been surpassed by several compression 
algorithms developed later. Huffman coding is also used as 
one of the stages for more complex compression techniques, 
eg. Cosine Transform based compression methods such as 
JPEG (Hashemian, 2003), Wavelet Transform Compression 
methods such as JPEG 2000 (Chang et al., 2006) or Deflate/ 
Inflate algorithms (Khademi & Krishnan, 2006). These coding 
techniques often have better compression capability although 
they are characterized by greater computing complexity and 
larger amount of memory that needs to be allocated during data 
compression and decompression. 


MBES and it’s data structure 


Multibeam sonars are a type of sonar systems widely used 
for surveillance and sensing in seas, oceans and other water 
containers (Mayer et al., 2007). They provide a wide swath 
survey of the ocean floor at each ping by pulsing the bottom 
with a series of soundings normal to the track of the vessel 
and recording the reflected echoes in an orientation parallel 
to the vessel. 

A typical multibeam sonar system consists of several 
elements, including transducer, Digital Signal Processor 
(DSP) and workstation, as shown in Fig. 1. The transducer, 
which uses linear or cylindrical arrays of receivers and 
projectors mounted at appropriate angles to each other 
(typically in so called the “Mills cross” arrangement, where the 
receive array lies in the athwart ship direction and the transmit 
array lies in the fore-aft direction) is responsible for ping 
generation and registration of backscattered acoustic wave 
signal. The registered data is transmitted to the DSP where 
the beamforming (based on sum-delay or Fourier transform 
algorithms), bottom detection (amplitude or phase detection 
algorithms), georeferencing, calibration and stabilization 
processes are performed. For each receive beam, the bottom 
echo from the intersection of the transmit and the receive 
beam footprint is detected. The entire system is managed by 


Ethernet 


Transducer 


Processor 
Fig. 1. The overview of multibeam sonar system structure 


the workstation, which is a PC unit operating under Windows 
or Linux. The workstation also enables storing and archiving 
of processed data. 

The minimum and maximum depth values for data 
acquisition depend on sonar model and vary between a few 
meters (min) to over 1000 m (max). Nominal beamwidth ranges 
from 0.5° to 1.5° x 0.5° to 1.5° and the frequency of performing 
the transmit-receive cycle (ping-rate) is 40 Hz. The commonly 
applied operating frequency of mutlibeam systems is between 
200 kHz and 400 kHz. Most mutlibeam systems offer integrated 
modules capable of acoustical imaging, mapping and data 
storage. The data recorded by these systems is stored in binary 
files defined by vendors. Although there are several formats of 
MBES data storage, their syntactic contents are similar due to 
the common architecture of multibeam systems. 

The echo arrival time and the angle of the receive beam 
provide information for bathymetry, and the backscattering 
strength is used for seabed imaging (Shah & Talukdar, 1999). 
The process of forming beams is called beamforming and has 
been widely described in literature e.g. by Rinhovde (1999) or 
Thorner (1990). In this context, sea-bottom charts and images 
are formed by appropriate processing and combining data from 
many consecutive swaths. 

The result of surveys made with the use of MBES systems 
is usually stored in two basic configurations: 

a) The data is processed and logged by dedicated operating 
software like QUINSy, HYPACK, PDS 2000 or other. 

In this case the results of surveys are kept in dedicated 

databases or database files. 

b) The raw records retrieved from MBES system are logged 
as a binary file in order to perform further analysis. 


In both of these configurations the MBES data is received 
as aset of binary records containing various information 
delivered by the multibeam system, such as bathymetry, water- 
column data, beam geometry data and sonar configuration data. 
Additionally, external sensors data such as CTD probe, GPS, 
Compass, Gyro and others may be also included. The exact 
structure of the dataset depends on the particular multibeam 
sonar model, survey purpose, configuration and connected 
external sensors. Nevertheless, the following characteristics 
are common for every multibeam sonar system: 

e The multibeam system must be connected to a workstation 
in order to collect and process the acquired data and present 
the results; 

e Data is delivered as a set of binary datagrams (records) 
via a connection between the sonar processor and the 
workstation; 

e The MBES produces several types of datagrams, each 
containing different types of information such as GPS data, 
bathymetry data, external sensors data (i.e. heave, peach, 
roll) and water column data; 

e No compression of MBES records is performed either by 
the sonar hardware or the specialized software. 


Although collecting water-column data is not mandatory 
during multibeam hydrographic cruises, it delivers a lot of 
additional information about the monitored marine environment. 
It is also used by diverse groups of researchers as an important 
component of survey results, delivering acoustic information 
about various underwater objects, as shown by Gerlotto et al. 
(1999) and Mayer et al. (2002) and for seafloor characterization 
and classification (Lubniewski, 2010). The water column data 
represents a set of signals acquired from multibeam sonar 
sensors in one sonar ping. The set of beamforming functions 
depends on the sonar specification. The datasets presented in 


the paper were acquired by RESON Seabat 7125 mutlibeam 
system which utilizes standard sum-delay beamforming. Two 
different visualizations of the water-column data acquired from 
this device are shown in Fig. 2. In this sonar the water-column 
is available as a two-dimensional binary array where the value 
of point defined by indices m and n is the 16-bit integer value 
of beam numbered m and sample numbered n (Marcus, 2007). 
The presented research was conveyed with the use of eight 
different datasets acquired by means of the RESON Seabat 
7125 MBES. The overall characteristics of this data are shown 


in Tables 1 and 2. 


50 100 150 200 250 300 


Fig. 2. Sample water-column data acquired 
by RESON 7125 Seabat Multibeam System. 


In order to explain the importance of efficient water- 
column data processing, a more detailed description needs to 
be given. Although the tested system generates 256 beams, 
in some multibeam systems the maximum number of beams 
can reach 3520 per one sonar head. The number of registered 
samples for individual beams depends on the depth and 
probing frequency and can reach up to 32 thousands, as shown 
eg. by Van Buren and Blue (1990), Xiaofeng and Wenjun 
(1998) or Foote et al. (2003). Since each value is stored as 
a 16-bit float or integer, the size of water column data can 
easily exceed 95% of all data collected by a multibeam system, 
as shown in Table 1. This is particularly apparent when the 
data is stored in linear scale, as it is the case in the datasets 
utilized in the presented work. 

In literature the value of water-column datagram samples 
is calculated according to the sonar equation: 


EL = SL + TS — 2TL 
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In other words: Received echo level (EL) is equal to 
transmitted source level (SL) plus target strength (TS) minus 
two-way transmission loss (2:TL) expressed in dB. The two- 
way transmission loss is: 


2TL=2 aR + 40 log R (re R=1m) 


where a is absorption coefficient (dB/m) and R is slant range 
(m). 

In a water column datagram, the received echo level is 
adjusted for system dependent parameters in addition to being 
amplified with a range dependent TVG function: 


Awe = EL— SL- 10 log Qiy ‘Qax + 
+2aR+40* logR+C 


where Q-y and Qgy are along-track (transmit) and across-track 
(receive) beam widths respectively (expressed in radians). 
These values are dependent on actual sonar model and operating 
parameters. For instance, in the RESON 7125 MBES which 
was used as a source of data during the presented work, the 
parameter Q,, was equal to 0.5° and Q,, equaled 1°. The 
parameter C represents gain selection of MBES transceiver used 
to adjust the signal level to system dependent parameters. 

In the presented formulas the echo levels calculations are 
expressed in dB. Also, the transducers accuracy used in the 
RESON 7125 MBES system is lower than 0.5 dB, as shown eg. 
by Van Buren and Blue (1990), Xiaofeng and Wenjun (1998), 
or Foote et al. (2003). This is the reason why some MBES 
manufacturers (eg. Kongsberg) apply conversion from linear 
scale to dB scale with 0.5dB re 1uP precision. That has been 
also applied during the presented research. 

In this context, the objectives of this work are to achieve 
improved data reduction ratio and offer reduced time of 
processing in comparison to standard archiving tools. 


(2) 


The proposed algorithm 


The developed algorithm is a modification of static Huffman 

codes dedicated for multibeam sonar data. It aims to: 

e enable efficient coding/decoding of large binary files 
containing MBES data; 

e lay foundations for a standard file format for compressed 
MBES data; 

e deliver an efficient archiving tool for online processing 
of particular MBES records, fulfilling the following 
requirement: 

- compression time of every record produced by MBES 
system must be shorter than the time of creating it by 
the sonar system; 

- program memory use must be unrelated to the size of 
compressed files, and be possibly small. 


The algorithm aims to improve MBES data processing speed 
in comparison to standard Huffman coders by assuming that the 
histograms of records and probabilities for particular symbols 
(byte values) to be encoded are similar within individual 
datagram types. Using this knowledge, the algorithm composes 
several Huffman trees based only on the first datagram of 
particular type delivered by the MBES system. 

According to (1), the entropy of a dataset (datagram) is the 
minimum average number of bits that can represent an encoded 
symbol using binary prefix sad which can be expressed as: 


-S plog, (pi) 


where p; is the probability AAR to symbol i that 
can also be defined as n/N where n; is the number of bytes 


E pattern = 
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containing value i in first pattern datagram and N is the total 
number of bytes. The latter is usually constant for every 
datagram type, as shown by Marcus (2007). According to 
Huffman theory, log,(p;) represents the number of bits required 
to encode symbol i. In this case, E atem is the entropy of the 
pattern datagram, which defines the Huffman tree. Once a tree 
is created for the analyzed type of datagram, the calculated 
Huffman structure is applied for encoding the rest of the 
records of particular type. In this case the entropy, calculated 
for a consecutive record is defined as: 


= -$ (pi + Ap; log» (pi) (3) 

i=0 
where Ap, is the difference of probability between the 
probability symbol i in the pattern datagram and the current 


datagram, and Yan =0 


To make daa uniquely decodable, additional information 
about the compressed records, such as the type of compressed 
datagram and its original and compressed sizes, must be saved 
to the result file. 

Fig. 3 depicts the entropy of the water-column data 
calculated by (3) on five sample datasets. As shown in the 
diagram, the entropy difference of at least 30 consecutive 
datagrams is approximately constant. Thus, assuming that the 
entropy is the lower bound for Huffman coding algorithm, 
compression ratio shouldn’t vary much between consecutive 
swaths. It is also worth mentioning that in some cases an 
increase of entropy was observed after several pings, which 
is aresult of on-device gain adjustment when forming beam 
geometry data. This will have a considerable impact on the 
presented results. 

It may be noticed that the number of datagram types 
produced by the MBES system is constant and, according to 
technical specification, does not exceed 20-25 in most cases. 
This issue is described in detail by Marcus (2007) and Van 
Buren & Blue (1990). The entropy characteristics shown in 
Fig. 3 allow the algorithm to assume that the consecutive 
records of the same type generated by the MBES system are 
similar in content. Thus, all consecutive datagrams of the same 
type are encoded using the same Huffman tree which was 
generated for the first datagram. Because the Huffman tree 
is generated only once for each datagram type, the algorithm 
is granted a considerable increase in performance without 
negative influence on compression ratio. This mechanism is 
depicted in Fig. 4. 

As it has been said in the preceding section, in some MBES 
systems (e.g. RESON 7k series) the water-column data is stored 
in linear scale, which causes unnecessary waste of storage 
space. In order to further optimize the size of the resulting data, 
its precision has been made a user-controlled parameter which 
can be set in the range of 0.1-1 dB. For best data reduction 
efficiency without losing information relevant for visualization 
or analysis of the water-column data, the parameter must be set 
in such a way as to match the accuracy of the MBES system. 
During the tests, this parameter was set to 0.5 dB, which is 
more than the native precision of the RESON 7125 MBES, 
as shown eg. by Van Buren and Blue (1990), Xiaofeng and 
Wenjun (1998), or Foote et al. (2003). 

An additional advantage of this approach is the fact that 
the proposed compressed file format enables access and 
decompression of particular records without exigency of 
decoding the whole dataset, while maintaining the efficiency 
of decompression. In this case, the only elements that must be 
read from the compressed file are the appropriate Huffman tree 
and the position of the compressed data. Reading the structure 
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Fig. 3. Entropy of consecutive water-column records according to (3) 
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Fig. 4. Block diagram of the proposed online compression algorithm 


of the compressed file and retrieving the information such 
as the total number of datagrams, original size of particular 
datagrams and their types is also possible without the need of 
decoding the whole compressed dataset. That is particularly 
useful when processing files bigger than 500 MB. In this case 
the decompression process would last considerably longer when 
using standard archiving tools. 


RESULTS 


A sample application of the proposed algorithm for 
multibeam data processing was implemented using C++ 
programming language. The software processes binary files 
containing 7k series multibeam sonar records. This section 


contains the results of testing the software as well as the 
information on the testing methodology and input datasets. 


Testing methodology and input data 
characteristics 


The basic acceptance criteria for the proposed method 
are: time of compression of particular datagram types, time of 
processing of whole dataset, and compression ratio. 

Tables 1 and 2 show the information on the datasets 
used for testing the compression algorithms. The files were 
acquired during scientific or research surveys in an area of 
Santa Barbara, California, USA, which is characterized by 
relatively flat and homogenous bathymetry. The data was 
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Tab. 1. Description of sample datasets acquired during research in the area of Santa Barbara, California, USA 


Percentage of Operating 
Dataset name water-column data frequency 
in entire dataset [kHz] 


20060719 _204657__ 7125 (400kHz).s7k 85.72% 


Sample rate Result 
[kHz] file size[ MB] 


20051014 185729 7125 (400kHz).s7k 83.09% 


20070720_170910.s7k 82.77% 


20070720_171128.s7k 82.82% 


20060201 190255 7125 (400kHz).s7k 83.07% 


20070315_184108.s7k 97.96% 


20070720_170721.s7k 82.84% 


20070831 185543.s7k 87.87% 


Tab. 2. Acoustic parameters of water-column data processed during the tests 


20060719 _204657__7125 | 20051014 185729 __7125 
20070720_170910.s7k 20070720_171128.s7k 
(400kHz).s7k (400kHz).s7k 


ZEN BAN TAS 


Number of pings: 146 Number of pings: 274 Number of pings: 1013 Number of pings: 429 
Number of samples = 1595 | Number of samples = 1829 | Number of samples = 401 | Number of samples = 401 
Number of beams = 512 Number of beams = 256 Number of beams = 301 Number of beams = 301 
Range selection [m] =35 | Range selection[m]=40 | Range selection [m] = 50 | Range selection [m] = 50 
Power selection [dB] = 218 | Power selection [dB] = 207 | Power selection [dB] = 185 | Power selection [dB] = 200 
Gain selection [dB] = 30 | Gain selection [dB] = 43 | Gain selection [dB] = 83 | Gain selection [dB] = 25 


20060201_190255__7125 
20070315_184108.s7k 20070720_170721.s7k 20070831_185543.s7k 
(400kHz).s7k 


a baw 


Number of pings: 278 Number of pings: 9 Number of pings: 843 Number of pings: 137 


Number of samples = 1839 | Number of samples = 1595 | Number of samples = 401 | Number of samples = 806 
Number of beams = 256 Number of beams = 512 Number of beams 301 Number of beams = 301 
Range selection [m]=40 | Range selection [m] = 300 | Range selection [m] = 50 | Range selection [m] = 100 
Power selection [dB] = 220 | Power selection [dB] = 220 | Power selection [dB] = 175 | Power selection [dB] = 220 
Gain selection [dB] = 30 | Gain selection [dB] = 19 | Gain selection [dB] = 40 | Gain selection [dB] = 10 
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acquired using maximum pinging frequency of 40 Hz. It is 
important to note that the surveys were made under various 
conditions, with varying place of measurement, depth, bottom 
type, sonar frequency, probing frequency, etc. This is to 
clearly show that the research was based on a broad range 
of datasets, which cover a lot of the possible MBES survey 
configurations. 

Table | clearly shows that appropriate methods of processing 
water-column data are the key to reducing the size of MBES 
records. This means that an efficient and fast algorithm for 
processing water-column data will guarantee optimal results 
of processing the entire file. Table 2 presents the information 
on acoustic parameters and sample visualizations for every 
water-column dataset processed during the tests. The images 
display raw frames of data without any form of post-processing 
such as fanshape correction. 


Measured efficiency of the algorithm 


The presented results were obtained on a standard desktop 
PC equipped with an Athlon X2 3600 + CPU paired with 
1 GB of RAM, running Windows XP SP3. For each dataset, 
Table 3 lists the time of compression for water-column records 
compared with the time in which the records were generated by 
the sonar system. Additionally the compression ratio acquired 
with the use of the proposed method and that acquired by the 
standard ZIP compression tool are also presented. The time of 
compression of the remaining elements ofa record (bathymetry, 
position etc.) is usually about 1-2 ms and is always shorter than 
the time of record generation. 

The proposed method can be also applied as an archiving tool 
for MBES data. Currently, the most commonly used tools for 


archive volume reduction are the well known file compression 
tools such as ZIP, RAR or the newly developed 7-ZIP. These 
tools do not utilize the information about the character of MBES 
data and consider the input simply as a series of bytes. Table 
4 presents a comparison of the proposed archiving tool with 
other standard data compression applications based on the total 
time of archiving the whole dataset. 


DISCUSSION 


As shown in Table 4, the proposed method offers superior 
archiving speed in all test cases, being on average 68% faster 
than the next best ZIP algorithm, while offering up to 111% 
smaller output file size. In most cases the proposed tool 
also offers greatest reduction of the input file size, with two 
exceptions (test cases no. 6 and 7) where 7-ZIP performs 15% 
better, but at a cost of up to 829% longer processing time. There 
is also one special case where 7-ZIP offers 91% better results. 
This is due to statistical change of datagram contents during the 
survey where the pattern datagram doesn’t describe the content 
of records registered at later stages of the survey. In this case 
additional analysis of the data is required to update Huffman 
trees describing the statistical information about the content 
of the records (as seen previously in Fig. 3). Even though the 
proposed algorithm fails to produce the smallest file size in this 
test case, it is 384% faster than 7-ZIP. 

To better understand the origins of the obtained results, 
a more in-depth analysis of the data used in test cases no. 6 and 
7 is required. Figures 5 and 6 show the differences between 
the contents of consecutive water-column datagrams in the 
problematic surveys, in the form of images representing the 
actual data and corresponding histograms. 


Tab. 3. Compression ratio yielded with proposed methods of particular water-column datagrams compared 
to time of generation of the records by RESON 7125 MBES system 


Dataset name 


20060719 204657 7125 (400kHz).s7k 80 


Average water-column 
record compression 
time [ms] 


Time of 
generation of 
water-column 

record [ms] 


150-200 


Compression 
ratio 


ZIP compression 
ratio 


0.28-0.32 


20070315 _184108.s7k 350-400 


450-500 0.31-0.33 


20070720_170721.s7k 15-20 


100-120 0.19-0.21 


20070831 185543.s7k 20-25 


300-350 0.32-0.35 


20051014 185729 7125 (400kHz).s7k 40 


70-80 0.26-0.27 


20070720_170910.s7k 12 


120-130 0.32-0.33 


20070720_171128.s7k 


120-130 0.23-0.28 


20060201 190255 7125 (400kHz).s7k 


70-80 0.30-0.31 


Tab. 4. Compression ratios and the time of archiving of selected compression tools in comparison to the proposed method 


Result 


Dataset name file size[ MB] 


20060719 204657 7125 (400kHz).s7k 


Proposed 
methods 
ratio/time [s] 


0.31/65 


ZIP ratio/ 
time [s] 


-ZIP ratio/ 
time[s] 


RAR ratio/ 
time [s] 


0.65 / 155 0.50 / 295 0.54 / 220 


20070315_184108.s7k 


0.79 / 35 0.63 / 73 0.69 / 42 0.31/12 


20070720_170721.s7k 


0.51/255 |0.23/11 min.| 0.44/380 0.44 / 172 


20070831_185543.s7k 


0.48 / 80 0.40/195 0.41 / 110 0.19/29 


20051014 185729 7125 (400kHz).s7k 


0.63/178 0.48 / 380 0.57 / 182 0.28/90 


20070720_170910.s7k 


0.53/305 |0.40/13 min.| 0.46/ 430 0.49 / 234 


20070720_171128.s7k 


0.52/125 |0.40/13 min.| 0.46/170 0.46/94 


1 
2 
3 
4 
5 
6 
7 
8 


20060201 190255 _7125 (400kHz).s7k 


0.66 / 170 0.54 / 340 0.61 / 225 0.31/82 
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x 5 7 80 


(3) 


Fig. 5. Water column data and corresponding histogram for consecutive records in 20070720_ 170910. 
Left — ping no. 1, middle — ping no. 150, right — ping no. 300 


EJ 


(1) 


Fig. 6. Water column data and corresponding histograms of records in 20070720_171128. Left — ping no. 1, middle — ping no. 100, right - ping no. 220 


According to (2), such parameters as TVG, gain, and power 
of the sonar, as well as others that can be automatically or 
manually adjusted during the survey, influence the statistical 
parameters ofthe data contents and thus can make compression 
less effective. The change of gain as a function of ping number in 
the two problematic datasets is shown in Fig. 7. a modification 
of the algorithm that would automatically update the Huffman 
trees when the compression of consecutive records becomes 
less effective is an important issue to consider, especially for 
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dataset structures similar to those of 20070720 170910 and 
20070720_ 171128. 

What is also important is the fact that the time of processing 
particular datagrams is shorter than the time interval between 
delivering consecutive records by the MBES system. Therefore, 
the proposed method can be used as on-the-fly record size 
reduction tool in routine research or commercial surveys, 
or efficient archiving tool for bigger multibeam sonar data 
warehouses. 
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Fig. 7. Gain selection presented as a function of ping 
no. in 20070720_170910 (left) and in 20070720_ 171128 (right) 


An additional advantage of the presented algorithm is that 
it is dedicated to processing MBES data, and thus it allows 
unique identification of the type of processed information. As 
a consequence, the proposed compressed file format enables 
selective access to data such as individual datagrams and their 
types without the need of decoding the whole compressed 
dataset. This ascertains considerable time savings in comparison 
to standard archiving tools, in particular for larger datasets. 

The versatility of the algorithm allows its application to on- 
the-fly MBES data processing. The presented performance results 
were obtained on a standard PC-class computer. Implementing 
and running the algorithm on a dedicated processor such as 
a CUDA-enabled GPU provides a substantially different level 
of performance, opening a new range of possible applications of 
the algorithm. Preliminary research shows that utilizing a GPU 
allows the algorithm to achieve semi-real-time compression of 
water column data. This in turn allows the compressed frames 
to be dynamically visualized in three dimensions via a Web- 
based Geographic Information System. a sample visualization 
of animated MBES swath in 3D is shown in Fig. 8. 

Because the contents of every MBES frame have been 
already processed, they may be easily analyzed according to 
the paradigms of Geovisual Analytics. Application of online 
analysis algorithms enables extraction of water-column 
data such as pelagic fish schools on-the-fly as well as their 
visualization in 3D context alongside digital terrain data 


Fig. 8. A three-dimensional visualization of processed sonar data 


such as bathymetry (Kulawiak 2010). Fig. 9 shows a three- 
dimensional visual analysis of processed water column data 
in the context of a standard two-dimensional map via WebGL. 
For this purpose, GPS location data from the vessel is used 
to localize the collected MBES information in geographical 
context in the form of a marker on a two-dimensional map 
constructed with the use of HTMLS. The contents of the data 
are then visualized inside the marker popup cloud on a three- 
dimensional canvas. The water-column data is streamed from 
the server in compressed form, which enables its dynamic 
visualization on the client. The data presentation has been set 
up to exclusively depict the ocean bottom as well as objects 
floating below the water surface. Changes in the shape of the 
bottom have been reflected by different colors of the three- 
dimensional model. The backscatter objects in water column 
data are colored in white. 


Fig. 9. Three-dimensional geovisual analysis 
of water-column backscatter data 


The proposed modification of Huffman coder, optimized for 
storing and compressing MBES records, enables its application 
in devices which process MBES data in real-time, while its 
efficiency rivals that of well-established compression tools 
like ZIP and RAR. However, the presented study of multibeam 
sonar data reduction and storage methods has been conveyed in 
a field that has not yet been properly explored. This will likely 
change as the amounts of data collected during research, as well 
as commercial surveys, will continuously increase. Properly 
designed MBES data processing algorithms may result in 
opening up new worlds of interactive exploration to a multitude 
of users. Despite their preliminary nature, the presented results 
show great promise of fast and efficient MBES data reduction, 
storage and retrieval. Some important factors influencing the 
compression ratio such as TVG, gain, and bottom type are yet to 
be considered. The presented MBES data reduction and storage 
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techniques can become a powerful tool for users of MBES 
systems, used as an efficient archiving tool for researchers, 
hydrographers, ecologists, fishers and commercial groups of 
interest, such as dredging industry or pipeline inspection and 
chart production companies. 
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ABSTRACT 
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£ 


This article presents a wider perspective on possibilities of Motorways of the Sea 
development in the Adriatic Sea, with a special emphasis on the eastern coast of the Adriatic 
Sea. A complete overview of Motorways of the Sea development has been presented, and 
key elements influencing the development of Motorways of the Sea were analysed. It has 
been ascertained that there exist a lot of bottlenecks in the connection to the transport 
infrastructure in South East Europe. Superannuated port infrastructure, unsuitable inland 
connections, and nonexistent IT tools for electronic data exchange between all participators 


in the logistics chain hinder dynamic development of Motorways of the Sea. All these elements were 

thoroughly analysed, and a proposal for a macro transport strategy suitable for South East Europe has 

been exposed. It has been accentuated that there is a future for Motorways of the Sea introduction but 

all described bottlenecks have to be further analysed and removed in a short period of time, to stimulate 
private sector for financial investments. 


Key words: Motorways of the Sea; port infrastructure; intermodality; 
inland connections; South East Europe; Adriatic Sea 


INTRODUCTION 


Today, intermodality and the use of sea transport route is 
an important topic of European transport policy. Consequently, 
these fields are subjects of detailed analysis and are included 
in the European development plans. Establishing intermodality 
and developing Motorways of the Sea to obtain sustainable 
passenger and cargo mobility are essential goals for the 
European Commission. Such a policy is not new, as it was 
developed in early nineties of last century, when a green paper 
on the impact of transport on the environment was released by 
the European Commission. In addition, a white paper on the 
common transport policy was presented in 1992, where the 
industry was called to support co-modality between different 
transport modes, especially between sea and rail transport, and 
to establish efficient transport services using intermodal nodes 
within the EU network. 

The idea and the concept of MoS were firstly developed in 
June 1992 when Viamare S.p.A. introduced the first initiative 
to shift road transport to the sea. It was a sea line between 
Genoa and Termini Imersi in Sicily that was the starting point 
of today’s common European strategy called Motorways of 
the Sea (MoS). According to Paixao [1], Bagchus and Kuipers 
[2] presented the concept “autostrade del mare” for the cargo 
shift from the road to the sea in the Portugal and Netherlands 
Corridor. This MoS concept started a dynamic development in 


Europe, as Grimaldi introduced the service between Genoa and 
Palermo, and the International Association of Turkish hauliers 
introduced a RO-RO (roll on-roll off) service between Turkey 
and Italy. Different services were established in the next years, 
first a RoPax (roll on-roll off and passenger) service between 
Ancona and Patras, followed by different services in the Baltic 
and North Sea region. 

Increased industrial interests stimulated the European 
Commission to highlight the MoS as a European concept to 
decrease the use of road transport, and consequently decrease 
the air pollution and traffic congestion in the European territory. 
Psaraftis [3] says that MoS will concentrate flows of freight on 
specific sea routes with the aim to establish new viable, regular 
and frequent maritime links between Member States, reduce 
road congestion, and improve the access to peripheral EU 
countries. For this reason the MoS concept has been promoted 
as an intra-community strategy and a cross-border project 
between all EU states. Special emphasis was put on ports 
with their infrastructure and hinterland connections, especially 
the railway network connections. European Commission 
recognised that MoS should be developed under a long-term 
transport strategy. Thus, different studies and project were 
launched and financially supported. Baird [4] recognises that 
from an academic standpoint these projects also contributed 
towards the development of an analytical framework to 
determine the feasibility of MoS. 
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Some EU countries already developed MoS services for 
regular sailing from their ports. A big challenge is foreseen for 
the region of South East Europe, especially for the Adriatic Sea 
region. This was also recognised by previous research studies 
done by Haralambous [5] and Tilling [6]. Consequently, new 
research activities are necessary to stimulate the development 
of MoS in South East Europe and in the Adriatic Sea. 


THE MOS CONCEPT 
Development of MoS in EU transport network 


The MoS were developed on a European Short Sea 
Shipping platform, as they are involving different operators 
and logistics players in a unique transport chain. The MoS 
concept was placed by EU as a strategic transport policy for 
the second time in 2001. Previously the European Commission 
worked on the ‘European Marine Motorways’ project under 
the Transport Research and Technological Development 
Program of the 4th Framework. Within this program, 
commercial viability of conventional and high-speed RO-RO 
services was studied in order to find an alternative to freight 
road transport. 

With a white paper on transport policy presented by 
European Commission in 2001 some important principles 
regarding the implementation of MoS were incorporated into 
a strategic European document. The use of European inland 
waterways to support more balanced usage of transport modes 
and to develop intermodality placed a special focus on port 
organization and infrastructure in use. In addition, European 
Commission suggested that the industry prioritises ports with 
good hinterland rail connections and supports the creation of 
important intermodal points in the Mediterranean, Atlantic and 
North Sea coasts. 

On this base, the study to define four MoS regions was 
carried out by the High Level Group, chaired by Karel van 
Miert and supported by European Commission [7]. According 
to the study presented in June 2003, the following four 
regions were proposed and later on accepted by the European 
Commission: 

- Motorway of the Baltic Sea (linking the Baltic Sea with 
Central and Western Europe, including the route through 
the North Sea/Baltic Sea Canal); 

- Motorway of the Sea of Western Europe (leading from 
Portugal and Spain via the Atlantic Arc to the North Sea 
and the Irish Sea); 

- Motorway of the Sea of South East Europe (connecting 
the Adriatic Sea to the Ionian Sea and the Eastern 
Mediterranean, including Cyprus); 

- Motorway of the Sea of South West Europe (western 
Mediterranean, connecting Spain, France, Italy, Malta and 
linking South East Europe and also the Black Sea). 


According to the experience from the Short Sea Shipping 
development, a lot of bottlenecks were identified. These are: 
complicated customs procedures, a lot of paper work, different 
operational and commercial documents in use, administrative 
procedures, expensive services’ implementation, mostly 
insufficient critical mass of cargo to be shifted in a short period, 
etc. During the last decade a lot of bottlenecks were eliminated, 
resulting in better starting conditions to establish and promote 
MoS. The European Commission played a crucial role in these 
actions, especially in reducing administrative barriers and in 
simplifying customs formalities on the entire transport route. 

The major unsettled problems are still connected to the 
ports, their infrastructures and the railway networks, which 
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are especially in eastern and southern Europe underdeveloped. 
Such infrastructure hardly supports cargo’s expectation for 
fast and cost effective services. In addition, the economy 
in the region is underdeveloped. This was ascertained also 
by Midelfart-Knarvik et al. [8], as they say that the MoS of 
the Baltic Sea and the MoS of Western Europe offer greater 
potentials in MoS development; as they can be supported 
by higher industrial distribution at a European level and by 
better trade relations. Moreover, Paixao [1] says that not all 
ports will be MoS gateways, only 87.2% of port authorities 
surveyed in a research of EU ports are predisposed to become 
MoS interfaces. 


Key elements influencing development of MoS 


It is of crucial importance to identify key elements that 
influence the development of MoS, because they have a deep 
impact on strategic decisions of shippers and intermodal 
operators. Based on the analyse of the European Marine 
Motorways project (EMMA) under the 4 Framework Program, 
Baird [9] exposes six main elements influencing adoption of 
MoS in EU transport framework. These elements are: 

- Price — MoS transport price versus actual road transport; 

- Service schedule — minimum service schedule of MoS 
should be a daily service to compete with actual road 
service; 

- Reliability — regular departures of vessels, with no weather 
influences; 

- Transit time — MoS should secure 24 to 48 hours door-to- 
door services to compete with actual road service; 

- Efficiency in port — factors as handling speed, cargo security, 
24-hour working in the port should be secured to maintain 
low-cost transport and service schedules; 

- On-board facilities — a range of services for drivers should 
be secured to accompany the cargo and to be included in 
the total transport costs. 


The price of transport services provided by MoS is 
a sensitive factor, because shippers or cargo owners might 
decide to shift their cargo from road transport to a MoS service, 
just in case the “all in” price is competitive. According to the 
United Kingdom Marine Motorways study performed by Napier 
University and Partners [10], a price of a RO-RO door-to-door 
service consists of the following three groups of costs: 
- 50% for the sea RO-RO transport; 
- 42% for road transport from door to port services and vice 
versa; 
- 8% for cargo handling. 


From the MoS price perspective it has to be emphasised that 
the price of inland transport is still very crucial in total costs for 
MoS services. Even if we shift cargo from road to the sea, over 
40% of all costs will be generated by inland services. Almost 
one tenth of all costs are produced by cargo handling. This 
poses very important pressure to the ports, in order to increase 
efficiency in every sphere of their organisation. Maintaining 
service schedule, service reliability and total door-to-door 
transit time are key elements to be secured by MoS. Shinghal 
and Fowkes [11] say that high frequency is a central attribute 
when customers are determining mode choice. 

In addition to the above six key elements, three important 
factors should be also mentioned. The first one is the degree of 
the developed hinterland infrastructure. The second one is the 
location of port terminals or hinterland intermodal terminals, 
and the third one is an IT connection between all participators 
in the logistics chains. 


The efficiency of ports, port services and especially port- 
hinterland connections are important elements for MoS [12]. 
The issue of developed or modernised hinterland infrastructure 
is for sure a big challenge for the underdeveloped regions and 
economies. Even if a super fast RO-RO service is introduced 
between the two ports, the problem of superannuated inland 
infrastructure can eliminate the advantage of fast sea service. 
The location of ports or hinterland intermodal terminals 
influences the transfer of cargo from different transport means. 
Ports or intermodal terminals must be positioned out of urban 
centres and massive traffic, to guarantee congestion-free access 
between sea and inland transport. Special transport corridors 
should be introduced to speed up the transfer of intermodal 
transport units out of the port or terminal systems. This must 
be supported by uniformed IT programs and connections, and 
by simplified documentation procedures. 

Based on all described key elements influencing the 
development of MoS services in the European transport network 
it can be stressed out that intermodality plays a crucial role in 
MoS development. The MoS can be more easily introduced 
in the areas where intermodal transport is already present and 
frequently used. Paulauskas and Bentzen [13] accentuate that 
a few parallel and independent transport routes that exclude any 
negative technical, political or economical influence are needed 
to develop MoS. Moreover, IT platforms and IT tools between 
port authorities, port operators, financial institutions, cargo 
owners, logistics operators etc. should be uniformed to simplify 
logistics processes in the door-to-door cargo delivery. 

For sure the financial perspective of supporting MoS 
development has to be considered too, as MoS implementation 
is directly connected to the financial support of the national 
economy of an interested state, or European founds. Of course, 
financial investments from the private sector are of significant 
importance, too. Beside this issue, the primary focus is on the 
infrastructure, which has to be previously secured. Only in 
this way will shippers and operators support the modal shift 
offered by MoS. 


TRANSPORT CHARACTERISTICS OF 
EAST MEDITERRAN AND ADRIATIC SEA 


The region of South East Europe is becoming very 
interesting research matter from the transport and logistics 
perspective. This is mainly connected with the analysed 
possibilities of further economical development and investment 
possibilities performed by global manufacturing and logistics 
companies. The region of Balkans is therefore presenting an 
opportunity to partly cover manufacturing necessities for the 
western economies, as in the era of global crisis transport costs 
have become an important issue of total price reduction. Far 
East production established in last decades is moving closer 
to the final consumer markets. Eastern Europe was recognized 
as an important alternative. These possibilities were already 
recognized by Fiat, Bosch, Continental, IKEA etc. 

South East Europe is attracting developed economies, but 
investments in massive manufacturing need a support from the 
transport and logistics sector in order to ensure “just in time” 
and “just in sequence” production and distribution of products 
to the final consumers. Consequently, the development of 
transport infrastructure and new transport concepts will play 
a significant role in future development of the region. 

Transport concepts like MoS can significantly support 
increased cargo flows of raw materials in import and finished 
product in export from the Balkans. The MoS will also fulfil 
higher transport quality expectations from the investors 
and logistics operators. Analysing the existing transport 


infrastructure and procedural bottlenecks can speed up the 
economy and industry to prioritise the settlement of the 
mentioned tasks. 


Analysing and removing existing bottlenecks 


Bottlenecks influencing the development of MoS in South 
East Europe were analysed in different research papers [5, 6, 
14, 15]. With the S.W.O.T analyse Haralambous [5] exposed six 
main weaknesses of MoS in South East Europe. These are: 

- Inadequate port infrastructure to support intermodality in 
all states, 

- Problematic port-hinterland connections, mainly with 
railway routes, 

- Insufficient railway infrastructure in all parts, delay in 
the realisation of the TEN railway projects in South East 
Europe, 

- Over-aged vessel fleet shipping in the area, unable to support 
fast modal shift, 

- Lack of synergy among states and economies, 

- Absence of large Trans-national logistic operators that could 
play the leading role for uniting interested actors. 


In addition, Haralambous [5] says that threats, such as 
congestion in core ports and lines, dependences of regional 
economies, mostly on road transport, and the risk of failure of 
viable schemes after the end of funding, due to high investment 
and operational costs, are influencing the development of MoS. 
Consequently, it can be stated that transport infrastructure 
and services directly related to it are presenting the main 
bottleneck. 

Tilling [6] also sees the main bottleneck in the underdeveloped 
infrastructure and complete transport sector, and exposes that 
drastic reform in the transport industry is still advocated as 
the only chance of the survival and economic growth of the 
South East European countries, as the capacity of countries to 
maintain the existing infrastructure and operational services is 
practically non-existent. 

Beside the superannuate infrastructure, an important 
bottleneck is foreseen in the field of cargo information exchange 
between ports, shipping lines and other logistics operators. 
The lack of used IT tools in the entire logistics chain definitely 
hinders the MoS development. This field should be thoroughly 
analysed, and a detailed action plan appropriate for South East 
Europe should be prepared shortly. 

As far as other bottlenecks are concerned, a severe 
doubt appears inevitable, as some economies and important 
companies are significantly depending on road transport and 
they might lose significant profit from shifting from road to sea 
or MoS. Consequently, strong resistance to change traditional 
transport routes and modes is being predicted. For this reason, 
it is important that the economy and industry recognise long 
term benefits of MoS implementation. 


Main possible routes of MoS 


Possible routes of MoS are developed according to market 
needs. Capacity utilization is an important element that should 
be considered. Beside this element the industry should consider 
also following elements: trade imbalance and variations in 
demand, competitive situation, state of the market, and types 
of customers and cargo [16]. 

Taking into account possible route characteristics different 
studies proposed arange of MoS corridors for South East 
Mediterranean and the Adriatic Sea. A study carried out 
for the “East Med MoS Master Plan” [17] funded by the 
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European Commission — DG TREN proposed nine main MoS 
corridors. Five corridors are important for the cargo flow in the 
Adriatic Sea. These corridors are: Igoumenitsa-Koper, Venice- 
(Igoumenitsa)-Patra-(Korinthos), Igoumenitsa-Ancona-Koper, 
Venice-Koper-Ploce and Malta-Venice. 

Just three MoS corridors include ports on the eastern coast 
of the Adriatic Sea. In all three corridors Koper port plays an 
important role. Meanwhile, from a range of other ports just port 
of Ploce is proposed for the line with Venice port. 

The Corridor Igoumenitsa-Koper should be an alternative 
to the strong road transport between south Europe and central 
Europe. Markets as Hungary, Austria and other central European 
countries should shift an important part of a daily truck line 
to the MoS. A yearly traffic of 400 000 to 600 000 tons has 
been foreseen till the year 2015. This MoS corridor organized 
by RO-RO service should be an alternative to the existing 
X. Pan-European Corridor and should use V. Pan-European 
Corridor to link Koper port with the hinterland destinations 
by rail, as already almost 70% of all cargo is transported from 
the port by rail. 

The described MoS corridor will face troubles to attract 
significant cargo flow, because the length of the journey 
from Greece to Slovenia via X. Pan-European Corridor is 
approximately 1 450 km. It can be foreseen that only one third 
of the actual road transport might be transferred to MoS service. 
Consequently, the total annually shifted traffic could reach up 
to 600 million ton/km. 

The second proposed corridor Igoumenitsa-Ancona-Koper 
is almost the same corridor as the previous one, but important 
cargo flow from central part of Italy should be collected by 
RO-RO service for central Europe, the Balkans, and for the 
western part of Greece. The study “East Med MoS Master 
Plan” shows that from Ancona to Koper and vice versa over 
2 million tons of goods should be shifted to MoS by 2015. 
Main countries interested in such service should be Slovenia, 
Croatia and Hungary, as they have strong import and export 
goods flow with the central part of Italy. 

The MoS corridor Venice-Koper-Ploce should be mainly 
established to connect the northern and central part of Italy to 
the Balkans, especially markets of Croatia and Serbia. With this 
RO-RO service around 800 000 to 1.3 million tons of goods 
might be shifted from roads to the sea. This corridor might 
also experience troubles with shifting significant cargo flow 
and passengers to MoS service, because the distance by road 
between northern Italy and Croatia is just 500 to 600 km, and 
these markets are connected mainly by highways. 

Beside the proposed MoS corridors, an important corridor 
should be additionally supported in the future. It is the RO-RO 
service between northern Adriatic ports (Trieste or Koper) 
with Albania and its port in Durres. Durres port is situated on 
important VIII. Pan-European Corridor and directly connects 
markets of FYR Macedonia, Bulgaria, Romania, Serbia and 
partly Greece. Although Albania is not an EU country the 
service should be seen as an important European intermodal 
link, which should offer and gain important benefits from 
modal shift. 


THE CASE STUDY OF THE ADRIATIC SEA 
REGION 


The case study performed on the Adriatic Sea region 
concentrates mainly on the eastern coast of the Adriatic Sea. 
A lot of bottlenecks have been foreseen by different authors, 
therefore a complete analyse of all main key elements 
influencing the MoS development is necessary. Consequently, 
all main ports were analysed and compared, including the 
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actually used infrastructure and IT tools or programs used by 
terminal operators. In addition, the inland infrastructure and 
hinterland terminals were analysed to formulate a strategy of 
faster MoS development in the Eastern Mediterranean and in 
the Adriatic Sea. Moreover, it becomes very important to set up 
a macro strategy for MoS development for a longer period, in 
order to attract financial support from the private sector. 


Port infrastructure and superstructure 


Port systems are one of the important elements in MoS 
development. They should cut down all waste operations, 
extra costs on door-to-door transportation, and perform lean 
modal shift. Problems of agile port operations have been 
deeply analysed and discussed by Vis and de Koster [18], 
Steenken et al. [19], and Notteboom et al. [20]. Paixao and 
Marlow [21] propose that ports must adopt a completely new 
logistic approach and agility to cope with new trends on the 
global market. 

It has been realised that there exist big differences between 
the ports and their terminals in sense of technical equipment, 
handling capacities and infrastructure. But it is valid for all 
these marine systems that are faced with more and more 
intermodal units to be handled in short time and at low cost. 
It is essential for operators to reduce unproductive time and to 
offer effective processes in order to be in line with constantly 
increasing competition among ports or terminals. 

According to the research of all main ports on the eastern 
coast of the Adriatic Sea, the infrastructure and superstructure 
is very poor. All maritime ports, except Koper and Piraeus, are 
underdeveloped on the sea side as well as on the land side. The 
dedicated intermodal infrastructure practically does not exist. 
In Tab. 1 basic infrastructure data for container and RO-RO 
terminals were collected. The data for Piraeus port were not 
included as this port is not representative due to well developed 
infrastructure and important role for containerised and RO-RO 
traffic in the Mediterranean. 

All five analysed ports together have just 10 specialised 
RO-RO ramps. In most cases, the quay or berths for RO-RO 
and container vessel are the same. Among these ports only 
Koper has suitable infrastructure to support efficient MoS 
services. Other ports can also be a gateway for MoS, but their 
infrastructure and superstructure cannot assure efficient modal 
shift which MoS should offer to the cargo owners. Beside the 
infrastructure, the superstructure of the maritime ports is also 
superannuated and actually cannot support higher volumes of 
intermodal cargo. Ports as Ploce, Bar and Durres do not have 
specialised container berth cranes. Almost the same situation 
is on the yard, where container gantry cranes are not in use. All 
handlings of intermodal units are performed by reachstakers 
or forklifts. 


Tab. 1. Basic port infrastructure data for container and RO-RO terminals 


Quay length | Draught | RO-RO 
[m] [m] 


Yard capacity 


ramps [m] 


65 000 
30 000 
1 135 000 
38 000 
135 500 


Moreover, it has to be jeopardising that in the South East 
Europe region there are no functioning modern hinterland 
intermodal terminals. This has been acknowledged also by 


Šakalys and Pal¥atis [22], as they came to the conclusion that 
a very small quantity of inland terminals exist in Southern 
and Eastern Europe, and they operate with limited handling 
equipment and on limited land area. The only operating 
railway stations that can accommodate containers and other 
intermodal units are situated in capital cities such as Belgrade, 
Zagreb, Ljubljana, Sarajevo, etc. These terminals do not use 
modern handling technologies, the degree of automation is very 
low, and static and dynamic capacities cannot secure optimal 
handling processes. 

Obviously the actual port infrastructure and superstructure 
on the eastern coast of the Adriatic Sea is not in a position to 
fulfil all key elements of MoS development. On the one hand, 
port systems cannot be cost effective, and on the other, they 
cannot maintain requested service schedules and so important 
transit times. The reliability of such MoS would be under a big 
question. 


Inland connections and infrastructure 


Intermodal transport and appropriate infrastructure is an 
important element of MoS development. It can be expressed as 
a basic platform to develop MoS, and it is directly linked with 
the health of the regional economy. Stronger economy uses 
intermodal transport advantages in order to achieve stronger 
long-term development. Consequently, such economies 
are directly investing in intermodal infrastructure. The 
comparison between Western European infrastructure with 
South Eastern one shows that the last one is underdeveloped 
and superannuates. 

The difference is not as big in road infrastructure as it is 
in railways. States of South East Europe invested in highways 


development between the capital cities. Highway connections 
with the maritime ports are practically nonexistent, except Koper 
port and to some extent the ports of Rijeka and Durres. 

The MoS for the South East Europe should be developed 
with connection to the rail transport. The railway should be 
used as the key transport solution between maritime terminals 
and hinterland hub terminals. Consequently, the road transport 
should be used only on shorter transport routes up to 150 km 
and on direct to door deliveries. Unfortunately, the use of the 
rail transport in South East Europe drastically decreased in the 
last decade, and at the same time the railway infrastructure was 
not modernised accordingly. As shown in Tab. 2, the total length 
of railway lines in South East Europe is 13 276 km and was 
reduced by over 200 km in the analysed countries. Meanwhile, 
in the same period Western Europe increased them, especially 
in Spain and Germany. 

The length of the road network in South East Europe is 
about 182 600 km. Just 3 910 km or 2.14% of all road networks 
are highways, where just Greece and Croatia are above this 
low average. Consequently, it can be stressed out that such an 
infrastructure does not allow high transport speed and cannot 
secure high security standards. 

Beside the fact that the road network is not developed and 
does not permit high transport speed, it must be accentuated 
that the rail network is even more superannuated. As analysed 
in Tab. 3, on an average the rail network represents just 6.8% 
of total inland transport network. Furthermore, the railway 
transport speed is in average below 40 km/h, in some sections 
even below 20 km/h. 

The analysis can be summarized that inland connections 
are not supporting modern logistics concepts. This can have 
positive and negative impacts on the development of MoS. 


Tab. 2. Total length of railway lines (km) [23] 


| | 2001 | 2002 | 2003 | 2004 | 2005 | 2006 


Belgium 


2008/2001 


Czech Rep. 


Germany 


Spain 
Sweden 


Netherlands 


Slovenia 


Croatia 


Macedonia 


Bulgaria 


Romania 


Greece 


Highways 
[km] 


Bosnia& Her. 


Road/rail 
[%] 


2.48 


Roads/Highways 
[%]) 


0.94 


Croatia 


Greece 


Macedonia 


Montenegro 


Serbia 


Slovenia 
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On the one hand, superannuated inland infrastructure calls for 
Mos, in order to shift cargo from inland transport to potentially 
better and faster sea transport. On the other hand, inland 
transport infrastructure is in worse conditions close to the ports, 
as some of them are still not connected with the highways. 
Consequently, road door-to-door services are preferred. 


IT systems and documentation procedures 


Modern information technologies and systems are very 
important tools to enhance the integration of intermodal 
operations and to speed up the modal shift at intermodal 
points. Efficient data exchanges between all logistics operators 
and cargo owners are especially important to establish just 
one IT data network, with uniformed quantity and quality of 
information flow. European Union and European Commission 
recognised this important goal a decade or two ago. The main 
aim is to support research which addresses practical problems 
of implementing the latest information technologies in 
intermodal logistic organisation, and help to break barriers by 
introducing cutting-edge technologies to the system, standards 
and services [24]. 

Different EU projects analysed actual needs of the 
industry to work in a real-time logistics information network. 
Important projects were INTRARTIP, MARNET, ITESIC, 
INFOLOG, etc, which provided proposals for full information 
network architecture, consisting of a uniformed information 
platform and standardised solutions to integrate Electronic 
Data Interchange (EDI) in the entire logistics chain. Services 
like MoS need up-to-date booking services, tracking and 
tracing management, paperless working processes, simplified 
administrative services, etc. 

Unfortunately, ports and hinterland terminals do not use 
uniformed IT systems, which should be easily connected 
through standardised information platform. Among the 
analysed ports and terminals only ports of Koper and Piraeus 
have integrated an EDI system that to some extent simplifies 
information exchange between shipping lines, ports and the 
inland transport operator. Ports of Rijeka, Ploce, Bar and 
Durres do not use modern IT platforms to simplify data 
exchange between logistics providers. Moreover, in Durres 
port all the procedures and direct orders to the port are still 
managed through hard copy documents, which can be hand 
over manually or as scanned documents by e-mail system. 
The same procedure is in practice also with other authorities 
as customs, phytosanitary inspections, police etc. 

IT communication is one of the key elements influencing 
MoS development, and the present situation in the ports and 
hinterland terminals in South East Europe is not satisfactory 
for fast MoS development. A two-pillar action plan should be 
made in a short time to remove the barriers of interoperability 
between different information systems. Firstly, ports must 
integrate an EDI practice in their processes, even if they 
are inhomogeneous with other operators. And secondly, 
a centralised communication node should be developed. Such 
a node should connect and integrate shipping lines, their agents, 
logistics providers, forwarders and cargo owners into a main 
‘umbrella’ module. Consequently, an IT platform should be 
successfully implemented, improving cargo information 
exchange and securing on-line cargo tracing and tracking. Thus, 
the needs and expectations from MoS should be provided. 


A need for macro transport strategy 


Based on the described findings, it is necessary that the 
region of South East Europe makes a common strategy on 
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developing intermodality in the region and on how MoS 
should be developed in a long term. According to the research, 
a short-term strategy would not bring necessary actions to 
actively develop MoS in the Balkans. In this situation two 
approaches are possible — a bottom-up or a top-down approach. 
Rodrigue et al. [25] says that with a top-down approach direct 
and fast actions can be implemented, but these actions should 
be executed by the governments or Port authorities. These 
actions should not be focused on road charges only as they 
may not achieve the desired shift from road to sea operations, 
particularly in countries where the economic and social impacts 
of road freight movement outweigh its external impacts and 
costs [26]. 

Parantainen and Merilainen [27] say that a bottom-up 
approach is the best way of developing concrete projects. 
Projects should be based on realistic estimates of transport 
volumes and create permanent and economically feasible 
transport chains. This is also valid for MoS, because the 
development of MoS should be done on a project base. It can 
be anticipated that with the bottom-up approach the industry 
can efficiently enter the sector of intermodality and modal shift. 
With such approach, the industry must present its expectations 
and needs to governmental institutions and dynamically co- 
create a macro transport development strategy in the region. 
With this it is possible to influence governmental investments 
in transport infrastructure and future legislation in the transport 
sector. The industry must press the Port authorities to introduce 
EDI platforms and to reduce documentation processes in 
business cooperation with the ports and hinterland terminals. 
Clark et al. [28] expose that legal restrictions can negatively 
affect port performance and the transport system as a whole. 
Thus, the industry must be active also in this very important 
field of transport regulation. 

Given the analysed infrastructure, superstructure, and 
current data exchange procedures in South East Europe it 
is necessary to develop a long-term transport policy for this 
European region. With a macro transport strategy an innovative 
environment should be established. According to Haralambous 
[5], the intermodal chain must provide free flow without any 
bottlenecks, physical or operational, and has to be a credible 
choice. The concentration on cargo flow is recommended to 
ensure viability, but also the promotion of ports from Maritime 
transport servers to intermodal nodal points that provide fast 
and low cost services, and optimum connections to land high 
level infrastructure are prerequisites for the success. 

Therefore, the MoS development mainly has to focus on 
infrastructure modernization to stimulate intermodality and, at 
the same time, the key elements and goals of modern, lean and 
green logistics should be considered accordingly. 


CONCLUSION 


The MoS initiative is an important goal supported and 
coordinated by the European Commission. The main goal is 
to reduce congestion on European roads, to shift cargo and 
passengers from land transport to the sea, to reduce external 
costs from the extent use of road transport, and to achieve 
long-term economic growth in the Community. Based on 
the above, the MoS concept has been promoted as an intra- 
community strategy, with special emphasis on ports, their 
infrastructure and hinterland connections, especially railway 
network connections. 

The MoS are very important “tools” also for the region 
of South East Europe. Global companies recognised Eastern 
Europe as an important alternative for Far East production. 
But investments in massive manufacturing need a support 


from the transport sector, in order to secure “just in time” and 
“just in sequence” production and distribution. The MoS can 
significantly support increased cargo flows of raw materials 
in import and finished product in export from the Balkans, 
but also need modernised port infrastructure as well as inland 
connections. 

The situation in South East Europe, especially in the 
Balkans, does not offer a short-term solution for massive MoS 
implementation. Firstly, the economy is very poor, especially 
in the era of global crisis. Secondly, the transport infrastructure 
and superstructure is superannuated, and without actual funds 
it is impossible to realise important investments, neither by the 
states nor private sector. 

The only adequate intermodal points are the ports of Koper 
and Piraeus, which invested in infrastructure and superstructure 
during the last decades. These ports can be an important gateway 
for MoS. Although they are well equipped and organised, the 
main bottleneck remains the inland connection by rail. Port 
of Koper is suffering from congestion on the rail transport, as 
over 70% of goods are transported from the port by rail. Other 
ports as Rijeka, Ploce, Bar and Durres do not have adequate 
infrastructure to importantly support MoS. Their main problem 
is also that they do not use uniformed IT systems, which should 
connect shipping lines, shipping agents, forwarders and other 
operators in a unique IT chain. 

Undoubtedly, the MoS have a future in South East Europe 
and in the Adriatic Sea. But it is necessary to develop a long- 
term transport policy for this European region to further analyse 
all described bottlenecks, in order to remove them shortly. In this 
way, a private sector should be importantly stimulated for direct 
financial investments. Consequently, a bottom-up approach 
of MoS development should be established. According to the 
past experience with MoS implementation, this is the only 
appropriate way to have long active MoS services. 
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The Ship Handling Research and Training Centre at Ilawa is owned by the Foundation for Safety of Navigation 
and Environment Protection, which is a joint venture between the Gdynia Maritime University, the Gdansk University 
of Technology and the City of Ilawa. 


Two main fields of activity of the Foundation are: 


Ə Training on ship handling. Since 1980 more than 2500 ship masters and pilots from 35 countries were trained at 
Itawa Centre. The Foundation for Safety of Navigation and Environment Protection, being non-profit organisation 
is reinvesting all spare funds in new facilities and each year to the existing facilities new models and new training 
areas were added. Existing training models each year are also modernised, that's why at present the Centre represents 
a modern facility perfectly capable to perform training on ship handling of shipmasters, pilots and tug masters. 


Research on ship's manoeuvrability. Many experimental and theoretical research programmes covering different 
problems of manoeuvrability (including human effect, harbour and waterway design) are successfully realised at 
the Centre. 


The Foundation possesses ISO 9001 quality certificate. 
Why training on ship handling? 


The safe handling of ships depends on many factors - on ship's manoeuvring characteristics, human factor (operator 
experience and skill, his behaviour in stressed situation, etc.), actual environmental conditions, and degree of water 
area restriction. 


Results of analysis of CRG (collisions, rammings and groundings) casualties show that in one third of all the 
human error is involved, and the same amount of CRG casualties is attributed to the poor controllability of ships. 
Training on ship handling is largely recommended by IMO as one of the most effective method for improving the 
safety at sea. The goal of the above training is to gain theoretical and practical knowledge on ship handling in a wide 
number of different situations met in practice at sea. 


For further information please contact: 
The Foundation for Safety of Navigation and Environment Protection 


Head office: Ship Handling Centre: 
36, Chrzanowskiego Street 14-200 ILAWA-KAMIONKA, POLAND 
80-278 GDANSK, POLAND tel./fax: +48 (0) 89 648 74 90 
tel./fax: +48 (0) 58 341 59 19 e-mail: office@ilawashiphandling.com.pl 

e-mail: office@portilawa.com 
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